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Abstract. 


A  spectral  analysis  and  normal  mode  expansions  are  developed  for 
the  acoustic  propagator 


Au  -  -c2(y)  p(y)  7  •  (p"  (y)  7u) 

of  a  stratified  fluid  with  sound  speed  c(y)  and  density  p(y)  at  depth  y. 
For  an  infinite  fluid  It  is  assumed  that  the  (in  general  discontinuous) 
functions  c(y),  p(y)  are  uniformly  positive  and  bounded  and  satisfy 

f±oo  *±00 

*  |c(y)  -  c(±<»)|  dy  <  «,  ±  |  P (y)  —  p(±°°)  |  dy  <  °°. 

o  '  o 

Semi-infinite  and  finite  fluid  layers  are  also  treated.  The  work 
provides  a  basis  for  the  analysis  of  transient  and  steady-state  sound 
fields  in  such  fluids. 


I 


§1.  Introduction. 

^  This  paper  presents  a  spectral  analysis  of  the  acoustic  fields 
in  stationary  plane  stratified  fluids  whose  densities  and  sound  speeds 
are  functions  of  the  depth.  The  analysis  is  based  on  families  of  normal 
mode  fields  that  have  simple  physical  interpretations. 

The  acoustic  field  in  such  a  fluid  may  be  described  by  an  acous¬ 
tic  potential  or  by  the  excess  pressure.  Each  of  these  is  a  real  valued 
function  u(t,x,y)  that  satisfies  the  wave  equation  [4,  13] 

(1.1)  0  -  c2(y)  p(y)  7  •  (p-1 (y)  Vu)  -  0 


where  t  is  a  time  coordinate,  x  *  (Xj,x2)  are  rectangular  coordinates  in 
a  horizontal  plane,  y  is  a  vertical  depth  coordinate  and 
V  ■  (3/3x! ,3/3x2,3/3y) .  c(y)  and  p(y)  are  the  variable  sound  speed  and 
density,  respectively,  and  p-1(y)  ■  l/p(y). 

The  paper  presents  a  spectral  analysis  of  the  acoustic  propagator 

(1.2)  Au  -  -c2 (y)  p (y)  V  •  (p~l (y)  Vu) 

for  the  cases  of  an  unlimited  fluid  ((Xj,x2,y)  €  R3) ,  a  semi- infinite 
layer  ((xt,x2)  e  R2,  0  <  y  <  “)  and  a  finite  layer  ((xltx2)  €  R2 , 

0  <  y  <  h) .  Only  the  first  case  is  presented  in  detail.  The  modifica¬ 
tions  required  in  the  second  and  third  cases  are  described  in  §9  at  the 
end  of  the  paper. 

In  the  case  of  an  unlimited  fluid  p(y)  and  c(y)  will  be  assumed 
to  be  Lebesgue  measurable  functions  that  satisfy 

1 


2 


(1.3)  0  <  Pm  <  p(y)  <  pjj  <  ®,  0  <  cm  <  c(y)  <  ^  <  », 
for  all  y€  R,  and 

#+00  #+00 

(1.4)  ±1  |p(y)  -  p(±“)|  dy  <  ±1  |c(y)  -  c(±»)|  dy  <  «, 

0  *  0 

where  pm,  pjj,  cm,  Cjj,  p(±°°)  and  c(±“)  are  constants.  It  is  clear  that 
(1.3),  (1.4)  imply 

(1.5)  Pm  <  P(±«)  <  PM,  Cm  <  C(±»)  <  Cjj. 

The  spectral  analysis  of  A  will  be  based  on  the  observation  that 
it  is  formally  selfadjoint  with  respect  to  the  scalar  product  defined  by 

(1.6)  (u.v)  -  [  u(x,y)  v(x,y)  c-2 (y)  p-1(y)  dxdy  , 

JR3 

where  u  is  the  complex  conjugate  of  u.  This  suggests  the  introduction 
of  the  Hilbert  space 

(1.7)  JC  ■  L2(R3,c_2(y)  p_l(y)  dxdy) 

with  scalar  product  (1.6).  Note  that  (1.3)  Implies  that  K  is  equivalent 
as  a  normed  space  to  the  usual  Lebesgue  space  L2(R3)>  although  they  are 
distinct  as  Hilbert  spaces. 

A  selfadjoint  realization  of  A  in  3C  is  obtained  by  defining  the 
domain  of  A  in  K  to  be 

(1.8)  D(A)  -  Ll(R3)  n  {u  |  V  •  (p"l(y)  Vu)  S  L2(R3)} 


where 
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(1.9)  Lj (R3)  -  L2(R3)  n  {u  |  7u  €  L2(R3)} 

is  the  usual  first  Sobolev  space  [1].  All  the  differential  operations 
in  (1.8),  (1.9)  are  to  be  understood  in  the  sense  of  the  theory  of 
distributions.  The  linear  operator  in  X  defined  by  (1.2),  (1.8) 
satisfies 

(1.10)  A  -  A*  >  0 

iff 

where  A  is  the  adjoint  of  A  with  respect  to  the  scalar  product  (1.6). 

A  proof  of  (1.10)  may  be  given  by  the  method  employed  in  [17]. 
Alternatively,  (1.10)  may  be  verified  by  showing  that  A  is  the  operator 
in  X  associated  with  the  sesquilinear  form  A  in  X  defined  by 

(1.11)  0(A)  -  L|(R3)  C  Jf 
and 


(1.12) 


A(u,v) 


L 


Vu  •  Vv  p~  (y)  dxdy; 


see  [8,  p.  322].  A 

The  principal  result  of  this  report  is  a  construction  of  the 
spectral  family  of  A  by  means  of  a  normal  mode  expansion.  The  spectrum 
of  A  is  continuous  and  contains  no  imbedded  eigenvalues.  This  fact, 
which  is  verified  below,  implies  that  the  normal  mode  functions  of  A 
must  be  generalized  eigenfunctions;  that  is,  solutions  $  of  the 
differential  equation 


(1.13) 


-c2(y)  p(y)  7  •  (p~*(y)  V^(x,y))  -  \ip(x,y) 
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chat  are  not  in  3C.  Solution  of  (1.13)  by  separation  of  variables  leads 
to  solutions  of  the  form 

(1.14)  <|>(x,y)  -  eip  x  iKy).  P'“  (PX,P2)  €  R2, 

where  p  •  x  =*  PjXx  +  p2x2  an<^  *Ky)  is  a  solution  of  the  equation 


(1.15) 


-c2(y)  p(y)  ~  |p-1(y)  -  |p|2^j  ■ 


with  |p | 2  •  p2  +  p2. 


The  operator  A 


defined  by 


(1.16) 


Ay  * 


will  be  called  the  reduced  acoustic  propagator.  For  every  y  >  0,  A^  has 
a  selfadjoint  realization  in  the  Hilbert  space 

(1.17)  3C(R)  -  L2 (R, c-2 (y)  p-1  (y)  dy) . 

The  domain  of  A^  is  the  set 

(1.18)  D(Ay)  -  L*(R)  n  |  j-  [p“(y)  Mj  e  L2(R)}. 

The  properties 

(1.19)  -  A*  >  c2  y2 

can  be  verified  by  showing  that  A^  is  the  operator  in  JC(R)  associated 
with  the  sesquilinear  form  in  JC(R)  defined  by 


.IL.IBII  "'■■1L- 
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(1.20) 

D(V  ■  L2(R)  C  3f(R) 

(1.21) 

Ay(<f>.^)  *  j  dy  +  y2  ^  p~1(y)  dy 

The  spectral  analysis  of  A  will  be  derived  from  that  of  A^.  The 
main  steps  of  the  analysis  are  the  following.  First,  hypotheses  (1.3) 
and  (1.4)  are  used  to  construct  solutions  of  A^  <j>  *  A<j>  that  have 
prescribed  asymptotic  behaviors  for  y  ■*  ±«°.  Second,  these  solutions  are 
used  to  construct  an  eigenfunction  expansion  for  A^.  The  construction  is 
based  on  the  Weyl-Kodaira  theory  of  singular  Sturm-Liouville  operators. 
Finally,  the  expansion  for  A^  and  Fourier  analysis  in  the  variables  xt, 
x2  are  used  to  construct  a  spectral  representation  for  A.  This  method 
has  been  applied  to  the  special  cases  of  the  Pekeris  and  Epstein  profiles 
[5,  17]  where  explicit  representations  of  the  solutions  of  A^  $  »  A<j>  by 
means  of  elementary  functions  are  available.  Thus  the  main  technical 
advance  in  the  present  work  is  the  construction,  for  the  class  of 
density  and  sound  speed  profiles  defined  by  (1.3)  and  (1.4)  of  solutions 
of  Ay  <j>  ■  A<j>  that  have  prescribed  asymptotic  behaviors  for  y  -*•  ±»  and 
sufficient  regularity  in  the  parameters  A  and  y  to  permit  application  of 
the  methods  of  [5,  17]. 

The  remainder  of  the  introduction  contains  a  description  of  the 
eigenfunctions  and  generalized  eigenfunctions  of  A^,  the  corresponding 
normal  mode  functions  for  A  and  their  physical  interpretations. 

The  limiting  form  for  y  -*■  ±«  0f  the  equation  A^  <J>  -  A<j>  is  the 

equation 


ivf.r 


6 


(1.22) 


2 

+  (Xc_2(±»)  -  y2)  <j>  -  0 


whose  general  solution  when  X  #  c2(±°°)y2  can  be  written 

f  *  -  Cl  eq±(ti’A)y  +  c2  or 


(1.23) 


where 


(1.24) 


and 


(1.25) 


siq±(U,A)y  +  c  e-iq+(U»X)y 


Cj  e  -  -r  u2 


q;(y,X)  -  (y2  -  Xc~2(±®))l/2  >  0 
q+(y,X)  -  ~iq+ (y  ,X) 


>  for  X  <  c2(±»)y2, 


q^(y,A)  -  i  q±(y,X) 

q+(y,X)  -  (X  c-2  (±«)  -  y2)I/z  >  o 


for  X  >  c2(±°°)y2 


In  particular,  the  solutions  are  oscillatory  when  X  >  c2(±°°)y2  and 
non-oscillatory  when  X  <  c2(±»)y2.  The  first  step  in  the  analysis  of 
will  be  to  construct  special  solutions  of  A^  $  •  X$  that  have  the 
asymptotic  forms 


(1.26) 


4>i(y»y»X)  *  e<lt'^’A^y  [1  +  o(l)],  y  ■*  +°° 

4>2(y.y»A)  *  e~q+^lJ,x^y[i  +  o(i)],  y  ■*  +» 

<J3(y,y,X)  -  eq-(U,X)y  +  0(1)]1  y  •+ 

0„(y.y»X)  -  e“q-^u,x^y[l  +  0(1)],  y  -*■  — 


X  i  c2(«)y2 , 

X  +  c2(-°°)y2. 
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It  follows  from  an  asymptotic  calculation  of  the  Wronskians  that  and 
4>2  are  a  solution  basis  for  A  if  ■  X<j>  when  X  #  c2(°°)p2,  while  and  <J>4 
are  a  basis  when  X  ^  c2(-")p2. 

The  nature  of  the  spectrum  and  eigenfunctions  of  can  be 
inferred  from  (1.26).  It  will  be  assumed  for  definiteness  that 


(1.27) 


c  (®)  <  c  (-«)  . 


It  follows  that  if  X  <  c2(»)y2  then  A^ 
if  (p2  and  <|>3  are  linearly  dependant. 

(1.28)  F(u,X)  -  p" 


=  X<J>  has  bounded  solutions  only 
Thus 


1  W(<J>2,c|>3)  =  0 


is  an  equation  for  the  eigenvalues  of  A^,  where  W  denotes  the  Wronskian 
and  p  1  W  is  independent  of  y.  The  corresponding  solutions 

(1.29)  ipk(y,h)  -  ak(u)  <J>2(y, U,Xfc(y))  =  a£(p)  (J)3  (y ,p,Xk(y)) , 


where  X  *  Xk(p)  is  a  root  of  (1.28),  are  square  integrable  on  R  and  hence 
are  eigenfunctions  of  A^.  Moreover,  A^  can  have  no  point  eigenvalues 
X  >  c2(°°)u2,  by  (1.26).  Thus  a0  (A^)  ,  the  point  spectrum  of  A^,  lies  in 
the  interval  [c2y2,c2(®)u2] •  Criteria  for  ^(A^)  to  be  empty,  finite  or 
countably  infinite  are  given  in  §3. 


It  will  be  shown  that  the  continuous  spectrum  of  A  is 


[c2  (»)y2 ,<»)  and  corresponding  generalized  eigenfunctions  will  be  deter¬ 
mined  from  (1.26).  For  c2(°°)y2  <  X  <  c2(-°°)y2  there  is  a  single  family 
of  generalized  eigenfunctions  of  the  form 


(1.30) 


<l>0(y.U,X)  -  a0(u,X)  4>j  (y,h,X) . 
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For  A  >  c2(-°°)y2  there  are  two  families  defined  by 


(1.31) 


^(y.U.A)  *  a+(y,A)  <J>4(y,u,X) 
<i»_(y.M.X)  -  a_(y,A)  <frx(y,y,A) 


It  will  be  shown  that  these  functions  have  the  following  asymptotic 
forms. 


(1.32) 


\(y.u)  ~ 


-q! (u,x.  (y))y 

cj(u)  e  K  ,  y  -  4-, 

q^(y,Ak(y))y 
ck(u)  e  ,  y 


(1.33)  i|>0(y,y,A)  -  c#(y,A) 


e-iq+(y,A)y  +  eiq+(lJ’X)y,  y  -K*>, 

T0(y,A)  eqI(lJ’X>y  ,  y  -  — , 


(1.34)  i|>+(y,y,A)  ~  c+(y,A) 


e-iq+(y,A)y  +  R+(PjA)  eiq+(y,A)y> 

y  -►  4-00^ 

l+(u M  ,  y  * 


(1.35)  <Ji_(y,p,A)  -  c_(y,A) 


T.(u.A)  ,  ,  - 

,H-(u.Wy  +  e-iq.(u.X)yi  y  . 


Here  ak(y) ,  a£(y),  a0(p,A),  a±(p,A),  ck(y),  c0(y,A),  c+(y,A),  R0(y,A), 
R+(U.A),  TQ(y,A)  and  T+(y,A)  are  functions  of  y  and  A  that  will  be 
calculated  below. 

Families  of  normal  mode  functions  for  A  may  be  constructed  from 
those  for  A^j  by  the  rule  (1.14).  The  following  notation  will  be  used. 


La’iuTigi'iA.  - -  ■  1  ..  " — .  na 
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(1-36)  <|>+(x,y,p,A)  =  (2tt)_1  eip'x  <|>+(y,  |p|  ,A) ,  (p,A)  e  £2, 


(1.37)  <Mx»y»P»*)  *  (2tt)-1  eip’x  ij>  (y,  |p|,A),  (p,A)  e  £20, 


(1.38) 


Vk(x,y,p)  *  (2tt)_1  eip'x  4>k(y.  I P I )  *  ^  k  >  1, 


The  parameter  domains  (2,  £2Q,  £2^  are  defined  by 

(1.39)  £2  *  {(p,  A)  e  R3  [  c2(-°°)  |p  | 2  <  A}, 


(1.40)  £2fl  =  {(p.  A)  €  R3  j  c 2 (°°)  |p  | 2  <  A  <  c2(-«>)  | p  |  2 } 


(1.41)  £2k  »  {p  6  R2  |  |p |  6  0k},  k  >  1, 


where  0,  is  the  set  of  y  >  0  for  which  A  has  a  kth  eigenvalue;  see  §6. 

The  three  families  have  different  wave-theoretic  interpretations 
that  are  characterized  by  their  asymptotic  behaviors.  Thus  for 
(p,A)  €  £2  one  has 


(1.42)  <J>+(x,y,p,A)  ~ 


c+(|p|,A) 


(1.43)  ij/_(x,y,p, A)  ~ 


2tt 


c  ( |p | , A) 


ei(p-x-q+y)  +  R+  ei(P.x+q+y)j  y  _  ^ 


T+e 


i(p*x-q_y) 


2tt 


T  eKp*x+q+y) 
gi(p*x+q_y)  +  R  ei(p*x-q_y) ^ 


,  y  +°°> 


where  q±  -  q±(|p|,A),  R+  -  R+(|p|,A),  etc.  Hence  i|/+(jc,y»P»A)  behaves 
for  y  ■*  +°°  like  an  incident  plane  wave  with  propagation  vector 


2£ 
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k^  *  (p,-q_p  plus  a  specularly  reflected  wave  with  propagation  vector 
*  (p»q+)>  while  for  y  -*■  it  behaves  like  a  pure  transmitted  plane 
wave  with  propagation  vector  lct  «  (p,-q_).  The  incident  and  transmitted 
plane  waves  can  be  shown  to  satisfy  Snell's  law  n(°°)  sin  0(°°) 

=  n(-°°)  sin  9(-°°)  where  8(»)  and  6(-°°)  are  the  angles  between  the  y-axis 
and  k^  and  kt>  respectively,  and  n(±°°)  *  c  1  (±°°) .  ^_(x,y,p,A)  has  a 

similar  interpretation. 

For  (p,A)  €  G0  one  has 


co<Ip! .A) 

(1.44)  >Mx,y,p,A) - 

2tt 


ei(p*x-q+y)  +  R  gKp-x+q+y)^  y  _  ^ 
T0  eip-x  £q:y  ,  y  -  -®. 


Hence  for  y  -*•  +»  <jj0  (x,y,p,A)  behaves  like  an  incident  plane  wave  plus  a 
specularly  reflected  wave  while  for  y  ■*•  -»  it  is  exponentially  damped. 
This  is  analogous  to  the  phenomenon  of  total  reflection  of  a  plane  wave 
in  a  homogeneous  medium  of  refractive  index  n(°°)  *  c  1  (<*>)  at  an  inter¬ 
face  with  a  medium  of  index  n(-«)  *  c-1  (-00)  <  n(“>) .  Indeed,  the 
condition  A  <  c2(-«°)  |p|2  is  equivalent  to  the  condition  for  total 
reflection:  n(°°)  sin  d(°°)  >n(-»). 

For  p  6  k  >  1,  one  has 


(1.45) 


\(x,y,p)  -  j 


Ck(|p|)  ip.X  ,q.y 

-  e  v  e  ,  y  -*■  +°°, 

2ir 


ck(lp|)  ip.x  q^y 
-  e  r  en  J  ,  y 

ft T 


Hence  the  functions  ^(x^p)  can  be  interpreted  as  guided  waves  that  are 
trapped  by  total  reflection  in  the  acoustic  duct  where  c(y)  <  c(±»). 


They  propagate  in  the  direction  k  -  (p,0)  parallel  to  the  duct  and 
decrease  exponentially  with  distance  from  it. 


The  coefficients  R+,  RQ  and  T+,  T0  in  (1.42),  (1.43),  (1.44)  may 
be  interpreted  as  reflection  and  transmission  coefficients,  respectively, 
for  the  scattering  of  plane  waves  by  the  stratified  fluid.  They  will  be 
shown  to  satisfy  the  conservation  laws 

(1.46)  q±  | R± | 2  +  q-  |T±| 2  -  q±,  |R#|  -  1. 

The  completeness  of  the  set  ,i^2 , . . . }  of  normal  mode  functions 

is  proved  in  §8  below. 

The  three  families  i|/+,  \p_  and  ip0  represent,  collectively,  the 
response  of  the  stratified  fluid  to  incident  plane  waves  exp  (i(p*x-qy)}, 
(p,q)  €  R3.  To  see  this  consider  the  mappings 

' 

(p.q)  ■  x+(p,A)  *  (p »q+( |p  | » A)) ,  (p,A)  e  n, 

(1.47)  (p,q)  -  Xfl(p,A)  -  (p,q+( |p| ,A)) ,  (p,A)  6  J20 , 

(p.q)  *  X_(p,A)  =  (p,-q_ ( |p | , A) ) ,  (p,A)  €  ft. 

X+  is  an  analytic  transformation  of  (2  onto  the  cone 

(1.48)  C+  -  { (p ,q)  |  q  >  a  |p|} 
where 

(1.49)  a  -  ((c(-»)/c(»))2  -  l)1/2>  0. 

Similarly,  X0  is  an  analytic  transformation  of  onto  the  cone 


(1.50) 


C0  ■  {(p,q)  |  0  <  q  <  a  [ p  | } 
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and  X  is  an  analytic  transformation  of  ft  onto  the  cone 

(1.51)  C_  -  { (p »q)  |  q  <  0}. 

Thus,  the  asymptotic  forms  of  4i+  and  \pg  for  y  -*■  ±»  show  that  t|>+(x,y,p,A) 
with  (p,A)  €  ft  is  the  response  of  the  fluid  to  a  plane  wave 
exp  {i(p*x  -  qy)>  with  (p,q)  £  C+,  ip0(x,y,p,\)  is  the  response  to  a  plane 
wave  with  (p,q)  €  C0  and  i<>_(x,y,p,A)  is  the  response  to  a  plane  wave 
with  (p,q)  6  C_.  Note  that 

(1.52)  R3  »  C+  U  C0  U  C_  U  N 

where  N  is  a  Lebesgue  null  set. 

The  interpretation  of  ijj+,  ij>_  and  \ 1>0  given  above  suggests  the 
introduction  of  a  composite  eigenfunction 

(1.53)  <J>+(x,y,p,q)  -  (2ir)”1  eip’x  4>+(y,p,q),  (p,q)  6  C+ u  C0  u  C_, 
where 

(2q)1/2c(»)  4>+(y,|pl»X),  (p,A)  -X+1(p,q),  (p,q)€C+, 

(1.54)  <{)+(y,p,q)  -■  (2q)  1/2  c(o»)  *J/q  (y,  |  p  [ ,  A)  ,  (P»A)  “Xj^p.q),  (p,q)€C„, 

( 2 1  q | ) 1/2  c(-“)  \i>_(y ,  1  p|  ,X)  ,  (p,A)  -  X_1(p,q), 

(p»q)^  c_. 

The  normalizing  factors  (2q)1^2  c(«)  and  (2jqj)1/2  c(-<»)  are  the  square  roots 
of  the  Jacobians  of  x"1 ,  X"1  and  X"1 .  The  function  4>+(x,y,p,q)  is  a 
solution  of  the  differential  equation 

(1.55)  A<|>+(*,p,q)  -  A(p»q)  $+(*,P,q) 


where 
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(1.56) 


*(p,q) 


c2(°»)(|p|2  +  q2),  (p.q)  €  C+  U  C0, 
c  2  (-00)  <  |  p  [  2  +  q2),  (p,q)  6  C_. 

Ics  asymptotic  behavior  is  described  by 

f  ei(P-»-W)  +  K  £  C+, 

tl(p-x-,y)  +  .KP-WW),  (p,q)  €  C(> 
I _  ,l(p-x^M.(|p|,X)y)  ,  (Pt,)6c., 


(1.57)  4>+(x,y,p,q)  ~  c(p,q) 


for  y  -*■  +a>  and 


fT  i(p**-q-(|p|  .*)y) 

+  e 


(1.58)  $+(x,y,p,q)  ~  c(p,q) 


(p.q)  e  c+, 
(p.q)  e  c0, 
l  ei(P*x-^)  +  R.  ei(p'X+qy),  (p.q)  €  c. 


T  eip‘x  eq:('pl’X)y, 


for  y  -*■  -®.  In  §8  it  is  shown  that  one  may  take 


(1.59) 


c(p,q)  - 


(2ir)~3/2  c(»)  p1^2  (®),  (p.q)  6  C+  u  C0, 
(2it)“3/2  c(-®)  p1/2  (-«),  (p.q)  €  C_, 


and  the  completeness  of  the  set  { , 4^ i  >'l>2 » •  •  • )  of  normal  mode  functions 
is  derived  from  that  of  »4'2 » •  •  •)  • 

Another  family  of  normal  mode  functions  for  A  is  defined  by 

(1.60)  <J>_(x,y,p,q)  ■  4>+(x,y,-p,q) ,  (p.q)  €  C+ U  C0  U  C_. 

It  is  clear  that  A  -  X(p,q)  and 


where 


(1.62)  <>_(y.p,q)  ■  <J>+(y.p.q). 

The  asymptotic  behavior  of  for  y  ■*  ±°»  may  be  derived  from  (1.57), 
(1.58),  (1.59).  It  is  given  by 

(1.63)  4>_(x,y,p,q)  ~  c(p,q) 

for  y  -*•  +“  and 

(1.64)  (J»_(x,y,p,q)  -  c(p,q) 

for  y  •*>  These  relations  clearly  imply  that  $_(x>y,P,q)  is  not 
simply  a  multiple  of  $+(x,y,p,q).  By  contrast  the  guided  mode  functions 
have  the  symmetry  property 

(1.65)  <l>k(x,y»p)  ■  Yk(x.y»~p)  ,  k  >  1, 

because  they  are  real-valued  and  depend  on  p  only  through  jpj. 

The  completeness  of  the  family  ,if»2 , . . . }  is  derived  from 

that  of  { ,^2  f  •  •  in  §8-  The  existence  of  the  two  families  and 
is  a  consequence  of  the  invariance  of  Che  wave  equation  (1.1)  under 


f  .«■>•*♦«>  *  E+  (p>q)  e  c+> 

ei(p-«-w)  +fj  (p,„  e  c,, 

T.  et<p-«-<u.(lp|,Wy)  >(p,,)€c.. 


[T+.i(p-«*i_(|p|.»y>  ,  (p.,)«c+. 

T<  .tp-.  .i:(|p|,Wy  c,, 

el(p*it+qy)  +  ei(p*x-qy) ,  (pi,,  s  c_, 
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time  reversal.  The  family  <p_  is  useful  in  the  construction  of  asymptotic 
solutions  for  t  +°°  of  (1.1);  see  [19]. 

The  remainder  of  this  paper  is  organized  as  follows.  §2  contains 
a  construction  of  the  special  solutions  (1.26)  of  Ay  $  ■  A<j>.  In  S3  the 
results  of  §2  are  applied  to  characterize  the  point  spectrum  and  contin¬ 
uous  spectrum  of  Ay.  The  eigenfunctions  and  generalized  eigenfunctions 
of  Ay  are  constructed  in  §4.  In  §5  the  Weyl-Kodaira  theory  is  applied 
to  construct  an  eigenfunction  representation  of  the  spectral  family  of 
A^.  §6  contains  an  analysis  of  the  dispersion  relation  (1.28)  and  the 

^-dependence  of  the  eigenvalues  X^(y)  of  Ay.  In  §7  the  results  of  § §5 
and  6  are  used  to  construct  a  normal  mode  representation  of  the  spectral 
family  of  A.  The  normal  mode  expansions  for  A  are  derived  in  §8.  The 
cases  of  semi- inf inite  and  finite  layers  are  discussed  in  §9.  §10 

contains  concluding  remarks  concerning  applications  and  extensions  of 
the  theory.  A  formulation  of  the  Weyl-Kodaira  theory  appropriate  for 
the  analysis  of  Ay  is  given  in  an  Appendix. 

The  analytical  work  needed  to  derive  and  fully  establish  normal 
mode  expansions  for  a  large  class  of  stratified  fluids  is  necessarily 
intricate  and  lengthy.  This  is  clear  from  examination  of  the  simple 
case  of  the  Pekeris  model  given  in  [17].  Therefore  to  make  the  work 
presented  here  as  accessible  as  possible  the  concepts  and  results  of 
each  section  are  formulated  in  the  first  portions  of  the  sections. 
Detailed  proofs  are  placed  at  the  ends  of  the  sections  and  may  be 
omitted  without  interrupting  the  exposition. 


fOTRmW 


§2.  Solutions  of  the  Equation  A^  0  -  £0. 

The  special  solutions  0j(y,U,X)  (j  -  1,2, 3,4)  described  in  §1 
are  constructed  in  this  section.  Analytic  continuations  of  these  func¬ 
tions  to  complex  values  of  X  are  used  in  §  §  3  and  5  for  the  calculation 
of  the  spectral  family  of  A^.  Hence  the  more  general  case  of  solutions 
of  A^  0  ■  £0  with  ;  £  C  will  be  treated. 

The  equation  A^  $  ■  C0  cannot  have  solutions  in  the  classical 
sense  unless  c(y)  and  p(y)  are  continuous  and  continuously  differentiable, 
respectively.  A  suitable  class  of  solutions  is  described  by  the 
following  definition  in  which  AC(I)  denotes  the  set  of  all  functions 
that  are  absolutely  continuous  with  respect  to  Lebesgue  measure  in  the 
interval  I  -  (a,b)  C  R. 

Definition.  A  function  0  :  I  ■  (a,b)  +  C  is  said  to  be  a 
solution  of 

(2.1)  Ay  <p(y)  =  -c2(y){p(y)(p“l(y)  <p’(y))'  -  u2<PC y)>  -  C0(y) 
in  the  interval  I  (where  0'  ■  d0/dy)  if  and  only  if 

(2.2)  0  6  AC(I) ,  p"V  €  AC (I) 


and  (2.1)  holds  for  almost  all  y  €  I. 

The  following  notation  will  be  used  in  the  definition  and 
construction  of  the  special  solutions  0j(y,U,O-  For  each  <  >  0 


(2.3) 


LOO  -  U  j  Re  ;  <  k1/2  }, 

R(k)  -  t  Re  C  >  <I/2  >. 


n 

I 
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(2.3  cont.)  ^(O  -  R(tc)  n  U  |  ±  Im  ?  >  0}. 

The  deflnlclons  (1.24),  (1.25)  will  be  extended  as  follows. 

' 

q+(y,S)  -  (5  c“2(±»)  -  y2)1/2 

(2.4)  -it/ 4  <  arg  q±(y,0  <  tt/4  ■  z,  €  R(c(±»)y) 
q±(y.O  -  i  q±(y.O 

and 

q+(y.O  *  (y2  -  ;  c~2(±«))1/2 

(2.5)  -tt/4  <  arg  q+(y,0  <  it/4  ■  C  6  L(c(±«)y) 

q+(y,c)  -  -i  q+(u»e). 

*”  * 

The  results  of  this  section  will  now  be  formulated. 

Theorem  2.1.  Under  hypotheses  (1.3),  (1.4)  on  p(y),  c(y)  there 
exist  functions 

(2.6)  :  R  x  R+  x  (L(c(°°)y)  u  R(c(«)y))  C,  j  -  1,2, 

(where  R+  ■  {y  |  y  >  0})  such  that  for  every  fixed  (y,c) 

€  R+  x  (t.(c(»)y)  u  R(c(°°)y)),  $j(y,y,£)  is  a  solution  of  (2.1)  for  y  g  R 
and  j  ■  1,2  and 

<J>i(y.y.c)  ■  “P  {q4.(^»S)y}[i  +  o(i)] 

(2.7)  ■  y  +a>, 
p_1(y)  <>J(y,y,;)  ■  p"1(°°)  q^.(y,?)  exp  {q^.(y,<;)y}[i  +  o(i)] 


and 
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<M y.U.O  a  exp  {-q|(y,OyHl  +  0(1)] 

(2.8)  y  -*•  +oo. 

P_1(y)  ^(y.y.C)  *  -p-1(°°)  q|(y,c)  exp  {-q^(y,Oy}[i  +  o(i)] 
Similarly,  there  exist  functions 

(2.9)  4>j  :  R  x  R+  x  (L(c(-«)p)  U  R(c(-<»)y))  -►  C,  j  -  3, A, 

such  that  for  every  fixed  (y,?)  €  R+  x  (L(c(-®)y)  U  R(c(-»)y),  4>j(y,y,0 
is  a  solution  of  (2.1)  for  y£  R  and  j  =*  3,  A  and 

<t>3(y,U,C)  *  exp  {q^(u,£)y}[l  +  o(l)] 

(2.10) 

p_1(y)  4>J(y,y»C)  -  p-1  (-00)  q.I(y,c)  exp  {q2(y,s)y}[l  +  o(i)] 

and 

<t>4(y»P.^)  *  exp  {-q^(y,?)y>[l  +  o(l)] 

(2.11) 

p  1  (y)4>4(y*y.O  *  -p  1  (-<,0)ql(y,c)  exp  {-q^(u,^)y)[i  +  o(i)] 

The  following  three  corollaries  describe  the  dependence  of  the 
solutions  <|>j(y,u,5)  on  the  parameters  y  and  £. 

Corollary  2.2.  The  functions  4>j(y,y,C)  satisfy 

(2.12)  dy  p-1^’  £  c[r  x  U  {(y,?)  |  £  e  L(c(»)y)}J 


for  j  -  1,  2  and 


20 


(2.13)  p‘V  £  c[r  x  [I  {(p,C)  ]  C  €  L(c(—)u)}] 

J  J  ^  y>0  i 


for  j  -  3,  4.  Moreover 


(2.14) 


♦i.  P-1^’  e  c[r  X  [J  {(y,?)  |  c  6  R+(c(«o)y)}), 
y>0  > 

<j>2»  P  (Jij  S  C  R  x  U  t(U  ,C)  |  5«  R"(c(“)y)}|, 
y>0  -1 


/  j  >  p  y  ^  ^  Vi|R  ^  U  t(u.C)  |  C  6  R-(c(-«)y)}) , 

*■  y>0  > 


p‘l*i  €  Cl 


P  e  c[r  X  (J  ( (y,c)  I  C  S  R+(c(-°°)y) }]  . 

1  y>0  > 

Corollary  2.3.  For  each  fixed  (y ,y)  €  R  x  r+  the  mappings 

(2*!5)  5  -*•  ^(y»y>0>  S  P-1(y)  <J>j(y»M»S) 

are  analytic  for 


(2.16) 


j  -  1.  C  6  L(c(°°)y)  u  R+(c(»)y)int, 

j  -  2,  5  6  L(c(»)y)  u  R“(c(»)y)int, 

j  -  3,  c  €  L(c(-»)y)  u  R'(c(-»)y)int, 

3  ■  4,  ?  £  L(c(— )y)  u  R+(c(-»)y)int, 


where  R±(K)int  -  R(k)  n  (c  |  ±  im  ?  >  o} . 


Corollary  2.4.  The  asymptotic  estimates  for  ^  and  p-I<J>j  of 
Theorem  2.1  hold  uniformly  for  (y,£)  in  any  compact  set  such  that  for 


Mb  " 
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j  -  i,  rt  c  |J  {(y.O  |  c  e  L(c(°°) y)  u  R+(c(°°)y)}, 


y>0 


(2.17) 


j  -  2,  r2  c  [J  {(y,?)  |  5  e  L(c(®)y)  u  R“(c(»)y)}, 


y>0 


j  =  3,  r3  c  U  t(v,5>  I  C  <=  L(c(-°°)u)  u  R“(c(-®)y)}, 


y>0 


j  =  4,  r„  c  (_J  {(y,C)  |  ?  e  L(c(-“)y)  u  R+(c(-®)y)}. 

y>0 


The  special  solutions  <j>j(y,y,?)  are  not,  in  general,  uniquely 
determined  by  the  asymptotic  conditions  (2.7),  (2.8),  (2.10),  (2.11). 
Indeed,  if  Re  q|(y,A)  >  0  (resp..  Re  q^(y,A)  <  0)  it  is  clear  that  any 
multiple  of  <p2  (resp.,  (j)!)  can  be  added  to  (^(resp.,  <|>2) .  A  similar 
remark  holds  for  <J>3  and  <f>4 .  However,  for  each  (  £  C  a  sub-dominant 
solution  (one  with  minimal  growth  at  y  =  00  or  y  =  -“>)  is  unique.  In 
particular,  since 


Re  q+(y,C)  >  0  for  (  £  L(c(±®)y) 

(2.18) 

Re  q+(y,£)  >  0  for  £  e  R”(c(±°°)y) 

Re  q;(y,?)  <  0  for  £  e  R+(c(±°°)y) 

one  can  prove 

Corollary  2.5.  The  solution  tp2  is  uniquely  determined  by  (2.8) 
for  all  £  e  L(c(°°)y)  u  R_(c(°°)y).  Similarly,  03  is  uniquely  determined 
by  (2.10)  for  £  €  L(c(-°°)y)  u  R  (c(-°°)y) ,  <f>i  is  uniquely  determined  by 
(2.7)  in  R+(c(°°)y)  and  $4  is  uniquely  determined  by  (2.11)  in  R+(c(-®)y). 


f 


A. 


- 


zac 


When  Re  £  ■  c2(±°°)y2  Theorem  2.1  provides  no  information  about 


the  asymptotic  behavior  for  y  -*■  ±°°  of  solutions  of  <j>  »  However, 

positive  results  can  be  obtained  by  strengthening  hypothesis  (1.4).  The 
following  extension  of  a  known  result  [11,  p.  209]  will  be  used  in  §3. 

Theorem  2.6.  Assume  that  p(y)  and  c(y)  satisfy  hypothesis  (1.3) 

and 


(2.19) 


| P (y)  -  p(°°)|  dy  <  “,  y2  | c(y)  -  c(«)  |  dy  <  ». 
o  a 


Then  there  exist  functions 


(2.20) 


fj  :  R  x  -  R.  i  =  1.2, 


such  that  for  every  y  €  R+  the  pair  4>j(y,y),  02(y»M)  is  a  solution  basis 


for  Ay  0  -  c2(°°)y2<t>. 


<f>1  (y>y)  =  i  +  o(i) 


(2.21) 


p_1(y)  4>J(y.y)  -  o(i) 


y  -*■  +”, 


and 


<My>u)  =  p(°°)y[i  +  °(D] 


(2.22) 


p-1(y)  4»,'(y.y)  -  l  +  o(i) 


y  -*■  +®. 


Lagrange's  formula  for  A  may  be  written 


(2.23) 


y2 


Ay  <J>  -  <p  Ay  ip>  c~2  (y)  p_1(y)  dy  -  [#]  (y2)  -  [#](yj) 


■Jr  >  li*‘- 
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where 


(2.24)  t#](y)  =  4>(y) (p-1  (y)  ^ ’ <y) )  -  >Ky)(p-1(y)  0’(y))- 


In  particular,  if  <j>  and  ip  are  solutions  of  0  =  £0,  A^  0  =  £0  on  an 
interval  I  then  [ Cy)  *  const,  on  I  and  [00] (y)  =  0  on  I  if  and  only 
if  <p  and  0  are  linearly  dependent  there.  By  combining  these  facts  and 
Theorem  2.1  one  can  show  that 


(2.25) 


[<M*»y>£)  <t>2 C •  »u» O ]  =  -2  p-1(“)  qj(y,0, 
[<l>3(*,y.O  4>4 ( • , y.Ol  -  -2  p-1(-°°)qUy,0 , 


which  imply 

Corollary  2.7.  The  pair  0x(y,y,O>  02(y>P«£)  is  a  solution 
basis  for  A^  0  =*  £0  for  all  (y,£)  €  R+  x  (L( c(°°)y)  u  R(c(°°)y)). 
Similarly,  the  pair  03(y,y,£)>  04(y,y,^)  is  a  solution  basis  for  all 
(y,s)  6  R+  x  (L(c(-°°)y)  u  R(c(-“)y)). 

This  completes  the  formulation  of  the  results  of  §2  and  the 
proofs  will  now  be  given.  The  method  of  proof  involves  replacing 
Ay  0  =  ?<(>  by  an  equivalent  first  order  system.  The  latter  can  be 
regarded  as  a  perturbation  of  the  corresponding  limit  systems  for 
y  ■*  ±°°.  In  this  way  integral  equations  are  established  for  solutions 
with  prescribed  asymptotic  behavior  for  y  -►  <*>  or  y  -*■  -<*>  and  these 
equations  are  solved  by  classical  Banach  space  methods.  This  technique 
for  constructing  solutions  with  prescribed  asymptotic  behavior  is  well 
known  -  see  for  example  [3,  p.  1408]  and  [12,  Ch.  VII]. 

A  first  order  system  equivalent  to  A^  0  =  £0.  If  0(y)  is  any 


solution  of  (2.1)  on  an  interval  I  and  if 


=  p~  (y)  <J>'(y) 


then  , 


(2.27) 


ip2  €  AC (I)  (cf.  (2.2))  and 

4>J(y)  -  P(y)  i|J2(y) 
M (y)  *  P-1 (y)[u2  - 


S  C_2(y)]  (y) 


for  almost  every  y  €  I.  Thus  the  column  vector  ip(y)  with  components 
ipi  (y) .  4<2  (y)  is  a  solution  of  the  first  order  linear  system 

(2.28)  i^'(y)  *  M(y,p,0  tp(y) 

where 


(2.29) 


M(y,U,S)  « 


0  p(y) 

P_l(y)[P2  -  C  c"2(y)]  0 


Conversely,  if  ip  €  AC(I)  is  a  solution  of  (2.28),  (2.29)  and  if  <J>(y) 
=  'Pj  (y)  then  is  a  solution  of  (2.1).  The  solutions  of  Theorem  2.1 
will  be  constructed  by  integrating  (2.28),  (2.29). 

The  limit  system  for  y  -*•  +°°  and  its  solutions.  By  replacing 
P(y)>  c(y)  in  (2.28),  (2.29)  by  p(«),  c(°°)  one  obtains  the  system 

(2.30)  *'(y)  -  M#(u,c)  4»(y) 


e 


(2.31) 


M0(y,C) 


o 


p(°°) 


p~  (°°)[y2  -  5  c~  (°°)] 


M0(P,5)  has  distinct  eigenvalues  q_|_(y,?)>  -q^(y,?)  for 
?  €  L(c(°°)y)  u  R(c(“)y).  The  columns  of 


(2.32) 


B(u.O  = 


p'l(°°)  qj(U,S) 


-p~J(»)  q|(U,C) 


are  corresponding  eigenvectors.  Hence 


(2.33) 


M0(y,O  B(y,C)  -  B(u,i;)  D(y,c) 


where 


(2.34) 


D(y.O  * 


<^(U,C)  0 

0  -q|(y,i;) 


System  (2.30),  (2.31)  may  be  integrated  by  the  substitution 


(2.35)  “  B(y,0  2- 

It  follows  that  z'(y)  -  D(y,S)  z(y),  whence 


(2.36) 


and  therefore 


zL(y)  -  Cj  exp  {q|(y,£)  y) 
z2(y)  “  c2  exp  (-q.J.(y,£)  y> 
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(2.37) 


<My) 

<l>2  (y) 


Cj  exp  {q^.(y,?)y}  +  c2  exp  {-q_J.(y , ^)y } 

p-1(°°)  q+(U,S)(ci  exp  {q_{.(li,Oy}  -  c2  exp  {-q|(y,?)y}) 


where  c: ,  c2  are  constants  of  integration. 

Application  of  perturbation  theory.  System  (2.28)  may  be 
regarded  as  a  perturbation  of  the  limit  system  (2.30).  Thus  if 
N(y,U»C)  is  defined  by 


(2.38) 


M(y,u,0  =  M„(u,C)  +  N(y,p,0 


then 


(2.39) 


N(y,u, C) 


0  a3(y) 

y2  a2(y)  +  C  a3(y)  0 


where 


(2.40) 


a3(y)  “  p(y)  -  pO*) 

a2(y)  -  p-1(y)  -  p  1(°°) 

a3(y)  ■  -[p_l(y)  c-2 (y)  -  p-1(®)  c-2(°°)] 


Note  that  each  of  these  functions  is  in  Lj(y„>“)  for  every  y0  €  R.  For 
a3  (y)  this  is  part  of  hypothesis  (1.4).  For  a2(y)  and  a3(y)  it  follows 
from  (1.3)  and  (1.4).  For  example,  one  can  write 


(2.41)  a2(y)  -  (-p-1 (®)  p-1(y))(p(y)  -  p(w)) 


which  exhibits  as  a  product  of  a  bounded  measurable  function  and  a 
function  in  Lj(y0,®). 
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On  combining  (2.28),  (2.38)  and  making  Che  substitution  (2.35), 
one  finds  that  (2.28)  is  equivalent  to  the  system 

(2.42)  z'(y)  -  D(y,?)  z(y)  +  E(y,y,0  z(y) 
where 

(2.43)  E(y,U,C)  -  B_1  (y.O  N(y,y.O  B(y,5) 

has  components  that  are  in  ^ j (y 0 , oo) .  Solutions  of  (2.42)  will  be  con¬ 
structed  which  are  asymptotically  equal,  for  y  -*•  +°°,  to  the  solutions 
(2.36)  of  z’  -  D(y.O  z. 

Proof  of  Theorem  2.1.  The  proof  will  be  given  for  the  function 
4> t  only.  The  remaining  cases  can  be  proved  by  the  same  method.  Solu¬ 
tions  of  (2.42)  are  related  to  the  corresponding  solutions  of  (2.1)  by 

♦  "  Z1  +  Z2 

(2.44) 

p"V  ■  P-1  (°°)  q|(Zj  -  z2),  q|  *  q|(y,C). 

Thus  <))  will  be  a  solution  of  (2.1)  that  satisfies  (2.7)  if  z  is  a  solu¬ 
tion  of  (2.42)  that  satisfies 

(2.45)  zx  -  exp  {q_J_  y}  Hj,  z2  -  exp  {q|  y)  n2 
and 

(2.46)  n,(y)  -  1  +  o(l),  ri2(y)  *  o(l)  for  y  -*■  ®. 

Equations  (2.42)  and  (2.45)  imply 
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(2.47) 


n{  *  eu  nx  +  e12  n2, 

n2  3  “2  q+  ^2  +  E21  n2  +  E22  n2» 


and  hence  by  integration 


nx(y)  -  Cj  +  |  E.(y')  ruty')  dy' 

Jya 

(2.48) 

n2(y)  -  exp  {-2  y}  c2  +  |  exp  (-2  q|(y  -  y’)}  E  (y')  n^y’)  dy' 

1  yx 

where  Cj,  c2>  yg,  yx  are  constants  and  the  summation  convention  has  been 
used  (j  is  summed  over  j  ■  1,2). 

Construction  of  d>i  for  C  6  L(c(»)u).  By  (2.18),  Re  q^(u,c)  >  0 
for  all  £  €  L(c(°°)y).  Thus  to  construct  a  solution  of  (2.47)  that 
satisfies  (2.46)  it  is  natural  to  choose  cx  *  1,  c2  *  0,  yfl  *  +»  and 
yx  finite  in  (2.48).  This  gives  the  system  of  integral  equations 


nt(y) 

(2.49) 

n2(y) 


1  -  r  Eij(y,)  nj(y,>  dy> 

|  exp  {-2  q_J_(y  -  y')}  E^y')  rtj(y’)  dy' 


y  >  yx- 


It  is  natural  to  study  system  (2.49)  in  the  space 


(2.50) 


X  -  CB([y,,»),c2) 


of  two-component  vector  functions  of  y  whose  components  are  continuous 
and  bounded  on  yt  <  y  <  <*>.  X  is  a  Banach  space  with  norm 
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(2.51) 

The  system  has  the  form 


Jnl  -  sup  ( l nx <y) |  +  |n2(y)(). 
y>yt 


(2.52) 


n(y)  -  n°  +  H  K(y.y')  n(y')  dy\  y  >  yx, 
Jyi 


where  n(y)  and  rj 0  are  column  vectors  with  components  (r^  (y) ,n2(y))  and 
(1,0),  respectively,  and  the  matrix  kernel  K(y,y')  is  defined  by 


(2.53) 


Kjj(y.y')  -  j 


(2.54) 


K2j(y,y') 


o  .  Yi  <  y’  <  y» 

^  -E^Cy'),  yx  5  y  <  y' 
exp  (-2  q|(y  -  y')}  E2j(y'),  yt  <  y’  <  y, 

o  »  yi  l  y  <  y'. 


and  j  ■  1,2.  The  conditions  Ejjt  €  Li  >°°)  and  Re  q|  >  0  imply  that  the 
operator  K  defined  by  (2.52),  (2.53)  and  (2.54)  maps  X  into  itself.  To 
show  that  K  is  a  bounded  operator  in  X  and  estimate  its  norm  note  that 


(2.55) 

|(Kn)j(y) I 

5  Inl  f  (|e  (y)|  +  |E  (y)|)  dy 
J7l  j  J 

for  j  - 

1, 2  and  all  y  >  yt . 

It  follows  that 

r  2 

(2.56) 

IKI 

-  I  dy‘ 

iyl  j »k*l  JK 

In  particular,  since  €  L^yj,00)  for  every  yQ  €  R,  (2.56)  Implies 
that  II  KB  <  1  for  every  sufficiently  large  yt .  For  such  a  value  of  yt 


m.itwwftwir 
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Che  equation 

(2.57)  n  -  n°  +  Kn 

has  a  unique  solution  n  e  X  given  by 


(2.58) 


l  *“0°. 

m“0 


Moreover,  (2.57),  or  equivalently  (2.49),  implies  that  n^y),  n2(y) 
satisfy  (2.47)  for  y  >  yJ4  These  functions  then  have  unique  continuations 
to  solutions  of  (2.47)  for  all  y  e  R,  by  the  classical  existence  and 
uniqueness  theory  for  linear  systems. 

Of  course,  n x  and  n2  are  functions  of  p  and  £  as  well  as  y 
because  q^  and  the  depend  on  these  variables.  The  solution  of 
Theorem  2.1  will  be  defined  by 


(2.59)  ^(y.y.O  -  exp  {^(vuOyKnj  (y,P,0  +  n2(y,u,s)) 

To  complete  the  proof  that  <(>!  is  the  desired  function  on 
R  *  R+  *  L(c(°°)u)  it  is  only  necessary  to  verify  that  (2.46)  holds  for 
each  (U,S)  €  R+  x  L(c(°°)y) .  It  is  clear  from  (2.49)  that 


,00  2 


(2.60) 


|nj(y)  -  i|  <  flnl  l  |e  . (y ’ )  |  dy'  -  o(i),  y  -*■  +». 
Jy  j-l  2 


For  n2,  (2.49)  implies  that 


1^2 (y) |  < 

(2.61) 


n  i 


exp  {-2q^(y-y')> 


2 


l 

j-l 


I e2  j (y ’ ) I dy  * 


+ 


.00  2 

l  |E  . (y') |dy' 

y2  j-i  2 
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for  every  y2  >  yx  and  every  y  >  y2.  Hence  for  any  fixed  y2  one  has 

r“  2  . 

(2.62)  lim  sup  |n,(y)l  <  Uni  l  |E  .(y')|  dy' 

r-  jy2  j»i  22 

because  Re  >  0.  Since  y2  in  (2.62)  is  arbitrary  it  follows  that 
n2(y)  -  o(l). 

Construction  of  4>i  for  C  S  R(c(°°)u).  By  (2.3)  R(c(«)y)  has  the 
decomposition 

(2.63)  K(c(»)u)  -  R+(c(»)y)  u  R“(c(«)p) int . 

Moreover,  for  £  €  R~(c(<»)y)  L  one  has  Re  q^_(p,£)  >  0  and  hence  the 
construction  of  the  preceding  case  is  valid.  In  the  complementary  case 
where  £  €  R+(c(°°)y)  one  has  Re  q^_(y  ,£)  <  0  and  it  is  permissable  to  take 
Cj  ■  1,  c2  *  0,  y0  ■  yx  «  «®  in  (2.48).  The  resulting  system  of  integral 
equations 


(2.64) 


nx(y) 


1  - 


j  Ej j (y’)  nj(y')dy' 


l2(y)  - 


00 

exp  {-2q_|_(y-y ' ) }  E  (y')  n.,(y')dy' 
Jy  J  J 


y  >  yx 


again  defines  an  equation  (2.57)  in  the  Banach  space  X.  Moreover, 

|exp  (-2q^(y-y')|  <  1  for  y  <  y'  <  »  and  (2.56)  is  again  valid.  It 
follows  that  for  yx  large  enough  (2.64)  has  a  unique  solution  given  by 
(2.58).  The  solution  has  a  unique  continuation  to  a  solution  of  (2.47) 
on  the  interval  y  €  R.  The  validity  of  the  asymptotic  condition  (2.46) 
is  obvious  from  (2.64);  cf.  (2.60). 
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Proof  of  Corollary  2.2.  Again  the  proof  will  be  given  for 
only.  Note  that  by  (2.39) 

(2.65)  N(y,U.O  -  N3(y)  +  y2  N2(y)  +  £  N3(y) 

where  the  components  of  (y)  are  in  Lj(y0,«°)  for  j  -  1,  2,  3  and  every 
y0  £  R.  Thus  by  (2.43) 

E(y,y,C) 

(2.66) 

-  B~l (U,C)N1 (y)B(y,c) +  y2B-1 (y,ON2 (y)B(y,s) +  £B-1 (y ,C)N3 (y)B(y,;) 

Proof  of  (2.12)  for  <fci.  Note  that  q|(y,0.  and  hence  also 
B(y, O  and  B-1(y,C)  are  continuous  functions  on  the  set 

(2.67)  U  ((y,C)  |  (6  L(c(“)y) }. 

y>0 

Thus  by  using  the  estimate  (2.56)  for  the  operator  K  *  K(y,i;)  in  X  one 
can  show  that  for  each  compact  subset  T  of  the  set  (2.67)  and  each  6  >  1 
there  is  a  constant  y!  *  y^T.S)  such  that,  taking  yl  -  y^T.S)  in  the 
definition  of  K(U,0,  one  has 

(2.68)  IK(y,01  <  5  for  all  (y.O  e  T. 

Hence  the  series  (2.58)  converges  uniformly  in  X  for  (y,?)  e  T  which 
implies  the  continuity  of  <j>j  and  p  1  <|> *  on  the  set  [ y 2 ,°°)  *  F .  Their 
continuity  on  R  *  T  then  follows  from  the  classical  theorem  on  the 
continuous  dependence  of  solutions  of  initial  value  problems  on 
parameters.  This  implies  the  result  (2.12)  for  <$j  because  T  was  an 
arbitrary  compact  subset  of  the  set  (2.67). 


I 
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Proof  of  (2.14)  for  <pi .  The  method  used  in  the  preceding  case 
is  applicable  to  the  operator  K(y,£)  in  X  defined  by  (2.64). 

Remark  on  Corollary  2.2.  The  argument  given  above  can  also  be 
used  to  show  that 


(2.69) 


v  p'S;  €  c 


R  *  U  Uu.C)  I  C  6  R“(c(»)y)int} 
y>0 


However,  the  continuity  of  pj  and  p  1 4> J  on  the  set 


(2.70)  R  x  U  I  R(c(°°)y)} 

y>0 

cannot  be  asserted  since  the  constructions  for  £  6  R+(c(°°)y)  and 
£  €  R-(c(®)y)int  are  different.  Indeed,  continuity  of  4>j  on  the  set 
(2.70)  is  not  to  be  expected  since  ^  is  not  uniquely  determined  when 
£  e  R_(c(»)y)int. 

Proof  of  Corollary  2.3.  The  components  of  the  matrix-valued 
function  E(y,y,£)  are  analytic  functions  of  £  £  L(c(°°)y)  u  R+(c(°°)y) int 
for  fixed  values  of  y,  y.  Hence  the  uniform  convergence  of  the  Neumann 
series  (2.58)  on  compact  subsets  of  this  set,  which  follows  from  the 
proof  of  Corollary  2.2,  implies  the  validity  of  Corollary  2.3  for  d>  1 . 

The  remaining  cases  can  be  proved  by  the  same  method. 

Proof  of  Corollary  2.4.  The  proof  will  be  given  for  the  case  of 
and  (y,£)  in  a  compact  subset  T  of  the  set  (2.67).  The  remaining 
cases  can  be  proved  similarly.  $j(y,y,£)  was  defined  by  (2.59)  and  the 
functions  rij(y,y,£)  satisfy 


EM*  ^ 


Hi(y,P.O 

(2.71) 

n2(y.u,c) 


1  *  -J  EjjCy'.y.C)  HjCy'.U.^)  dy', 

exp  {-2q|(p,C)(y-y')}  E^y'.y.O  n.j(y\y,S)  dy' 


It  must  be  shown  that  these  integrals  tend  to  zero  when  y  ■*  <*>,  uniformly 
for  (Vi,?)  €  T.  Now  (2.58)  and  the  estimate  (2.68)  from  the  proof  of 
Corollary  2.2  imply  that  ln(*,U.?)l  5  d  “  *$)  1  for  *11  (U»C)  S  r.  It 
follows  that  for  fixed  y'  >  y  >  yj(r,6)  one  has 


(2.72) 


Ekj(y',u,0  Hj (y'.y.O  I  1  iE(y\y,0«  in(y'.w.c)l 


<  (1  -  6)_1  IE(y’,U,?)I 


for  (y.O  €  r.  Now  the  continuity  of  B(y,£)  implies  that  there  is  a 
Y  -  y(0  such  that 

(2.73)  I B(Ut O I  IB_1(y, 6)1(1  +  y2  +  |c|)  <  Y  for  (y.O  G  T 


It  follows  from  (2.66)  that 


3 

(2.74)  |E(y',y,C)l  <  Y  l  I N,  (y ' )l  ,  (y  ,0  €  T. 

j-1  2 

Combining  (2.71),  (2.72)  and  (2.74)  gives 


(2.75) 


|n! (y.y.O  -  i|  <  y(i  -  6)"1 


.00  3 

l  «N.(y’)l  dy' 
7  j-1  J 


i 

1 


for  all  y  >  yt(r,6)  and  (y,5)  G  T.  Since  each  Nj  6  L2(y(|,00),  (2.75) 
implies  that  ni(y.y.S)  -  1  *  o(l)  uniformly  for  (y,c)  e  T. 
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The  case  of  r|2(y,P,C)  is  more  complicated.  Note  that 


|n2(y,u.O  |  <  Yd-<$) 

(2.76) 


.00  3 

exp  (-2  Re  q ' (u, O (y-y ’) >  \  UN . (y ' ) H  dy'. 

y  j-l  2 


Now  the  continuity  of  q_J_(P,C)  and  the  definition  of  L(c(“)y)  imply  that 
there  is  a  <  m  <(0  >  0  such  that 


(2.77)  2  Re  q|(u,?)  >  tc  >  0  for  all  (y.O  e  T . 

Combining  (2.76),  (2.77)  one  has,  if  Yi  *  y(1~<5)  *> 


|n2(y,u.O  I 

(2.78) 

<  Yx 


T2  exp  {-x(y-y') }  1  IN.  (y ' )I dy '  +  [  f  IN  (y')Idy' 
(Jy  j-1  J  Jy2  j-1  3 


00  3 

y2 


_  'i 


.CO  3 


.00  3 


exp  {-<(y-y2)  l  IN.(y')Idy'  +  £  ||N  (y')lldy' 

j-1  J  Jy2  j-i  J  J 


for  all  y2  >  y^T.iS),  y  >  y2  and  (y,C)  €  T . 

Now  let  e  >  0  be  given  and  choose  y2  =  y2(e,r,S)  >  yi(r,5)  such  that 


(2.79) 


and  hence 


Yi 


f°o  3 

l  IN, (y')I  dy'  <  e/2, 

Jy2  j-1  2 


(2.80) 


|n2(y,u.O  |  <  Y i 


exp  {-<(y-y2)> 


Yi 


3 

l  IN, (y ’ ) n  dy '  +  e/2 
j-1  2 


for  all  y  >  y2(e,I\6)  >  yj(r,5)  and  (u,C)  €  T.  Finally,  choose  a 
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y3(e,r,6)  >  y2(e,F,6)  such  that 


(2.81) 


Yi 


exp  {-<(y-y2)} 


'*1 


3 

l  i  N .  (y ' )  II  dy'  <  e/2 
j«!  J 


for  all  y  >  y3(e,r,5).  It  follows  from  (2.80)  and  (2.81)  that 
|n2(y,P,C)|  <  e  for  all  y  >  y3(e,r,<S)  and  (y,;)  e  T;  i.e.,  n2(y,U,C) 

=  o(l)  uniformly  for  (y,£)  €  F. 

Proof  of  Corollary  2.5.  It  will  be  shown  that  (^(y.y.S)  is 
uniquely  determined  by  (2.1)  and  (2.7)  when  C  £  R+(c(°°)y).  The  other 
cases  are  proved  similarly. 

Assume  that  for  some  £  6  I^(c(<=o) y)  there  are  two  solutions  of 
(2.1),  (2.7).  Then  their  difference  (j>(y)  would  satisfy  (2.1)  and 
<p(y)  *  o(l),  p" 1  (y)  <J>'  (y)  =  o(l)  because  Re  q|  <  0  for  t,  e  R+(c(°°)y). 
It  follows  that  the  corresponding  pair  n^y).  n2(y)»  defined  by  (2.44) 
and  (2.45),  would  necessarily  satisfy 


(2.82) 


nx(y) 


n2(y)  = 


Elj(y')  (y')  dy’ 


exp  {-2q_J_(y-y ' ) }  E2j(y')  ^(y')  dy' 


since  |exp  (-2q_J_(y-y ')  |  <  1  for  y  <  y'.  But  (2.82)  is  equivalent  to  the 
equation  n  *  Kn  in  X.  If  y3  is  chosen  so  large  that  ||KI  <1  then  q  =  Kn 
has  the  unique  solution  n(y)  =  0  for  y  >  y3 .  The  unique  continuation  of 
this  solution  of  (2.47)  is  then  zero  for  all  y  G  R.  Thus  <Hy)  =  0  for 
y  S  R,  which  proves  the  uniqueness . 

Proof  of  Theorem  2.6.  The  equation  <p  *  c2(°°)  y2$  is  equi¬ 
valent  under  the  mapping  (2.26)  with  the  system  (see  (2.31),  (2.38)) 
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(2.83) 


<P’  55  P(°°)ipz  +  BjCy)!^ 
<l»2  a  B2  (y)\p1 


where  Bj  (y)  =  ax  (y)  =  p(°°)  -  p(<*>)  and 


(2.84) 


B2(y)  =  V2  p-1(y)(l  -  c2(°°)c  2  (y) ) 

=  y2  P_1(y)  C~2(y)[c(y)  +  c(=°)][c(y)  -  c(°°)] 


It  follows  from  hypotheses  (1.3)  and  (2.19)  that 

(2.85)  Bj (y) ,  B2(y),  y  B2 (y) ,  y2  B2 (y)  €  Lj (yQ , «) 

for  every  yQ  €  R  and  every  y  >  0. 

Construction  of  4>i  •  Application  of  the  variation  of  constants 
formula  to  the  system  (2.83)  gives  the  integrated  form 


4ll  (y)  =  Cj  +  p(°°)  c  2y  +  {p(°°)(y-y')  B2(y')  ^j(y')  +  B^y')  ^2(y')}dy' 


(2.86) 

Cy)  =  c2  + 


B2(y’)  ipj(y’)  dy\ 


y0 


Now  <p1  will  satisfy  4>2  “  c2(®)y2cj>j  and  the  asymptotic  condition  (2.21) 
provided  that  m  4>lt  ^2  =  P  satisfies  (2.86)  and 


(2.87) 


<i>i  (y)  *  1  +  o(l) ,  \p2  (y)  =  o(l) ,  y  +  ». 


To  construct  such  a  solution  take  cx  *  1,  c2  =  0  and  y0  =  00  in  (2.86). 
This  gives  the  system 


38 


(2.88) 


My)  -  1 

<l>  2(y)  = 


(p  C00)  (y-y  ’ )  b2  (y ' )  4^!  (y  * )  +  b 


B2<y')  ^i(y')  dy' 


or 


(2.89) 


ip(y)  -  4>°  + 


K(y.y')  <Ky')  dy',  y 


where  iKy)  and  \p°  have  components  4»i  (y) ,  \p2(y)  and  1,  0 
K(y,y')  =  0  for  y  >  y'  and 


Ki i (y»y*)  =  p(°°)(y'-y)  Bz(y') 


(2.90) 


Kia(y,y’)  -  -B,(y’) 


K21(y,y’)  «  -B2(y’) 


K22(y.y')  =  o 


for  y  <  y'.  As  in  the  proof  of  Theorem  2.1,  one  has 


(2.91) 


and  hence 


I  (B|>)j(y)|  1  141 


(|Kjl(y,y,)|  +  |k.2(5 


I  (K 40 x (y)  |  <  I W 

<  I41I 


p(°°)  [  (y+y ') |b2 (y') jdy'  + 

jy 

2p (°°)  [  y'  |B2(y')|dy'  + 


i(y')  My’>)  dy’ 

>  yi , 

respectively. 


r,y')|)  dy' 

f  |Bl(y,)|dy* 

y 

1  Bj  (y  * )  |  dy  *  • 


(2.92) 


1 
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and 


(2.93) 


W2(y)|  <  1101! 


I B2 (y ’ ) I  dy' 


In  particular, 


*00  *oo  -00 

(2.94)  IK1  <  2p(°°)  I  y’|B2(y')|dy'  +  |Bl(y')  |dy*  +  |B2(y')|dy' 

*  XT  J  V  XT 


Thus  (2.85)  implies  that  K.  is  contractive  in  X  for  y  large  enough. 

Hence  (2.88)  has  a  unique  solution  on  [yj,00)  which  can  be  continued  as  a 
solution  of  (2.83)  to  all  y  e  R.  Moreover,  (2.88),  (2.92)  and  (2.93) 
imply  that  |^(y)  -  l|  -  IW^y)!  =  o(l)  and  |02(y)|  -  |W,(y)| 

*  o(l).  In  fact  (2.85)  implies  that  |02(y)j  *  o(y~2).  Thus  (2.87)  is 
satisfied. 

Construction  of  02.  02  will  satisfy  A^02  *  c2(°°)u24>2  and  the 

asymptotic  condition  (2.22)  provided 


(2.95) 


<My)  *  02  (y)  *  p(°°)y  m(y) 


02(y)  =  P_1(y)02(y)  =  n2(y) 


where 


(2.96) 


nx(y)  -  l  +  o(i),  n2(y)  -  l  +  o(i),  y  -*■  «. 


Substituting  in  (2.86)  with  Cj  =  0,  c2  -  1  and  y#  -  »  gives,  after 
simplification, 
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My)  "  l  - 

(2.97) 

My)  "  1  - 


{p(°°)  (y  '~y-  Cy')  +  p'^^y'^jCy'JnjCy'JJdy' 


p(oo)y'B2(y ')nj  (y')dy' . 


or 


(2.98)  M(y)  -  n°  +  f  K(y,y')  n(y')  dy',  y  >  yx 

Jy  i 

where  H(y)  and  r\°  have  components  Hj(y)»  n2(y)  and  1,1  respectively, 
K(y,y’)  =  0  for  y  >  y'  and 


Ku(y,y')  a  P(°°)(~y'  +  y  1y'2)B2(y') 


(2.99) 


Kl2(y.y')  *  -p'1  My"1  Bj(y') 


Kji(y,y’)  *  -p(°°)  y'  B2(y') 


K22(y,y')  -  o 


for  y  <  y'.  It  follows  from  (2.91)  that 


|  (Kn)  j  (y)  1  <  l  oi 
(2.100) 
and 


P( 


°)  |  (y,+y”ly,2)|B2(y,)|dy’  +  p"l(»)y-l[  |Bj(y')|dy’ 

*  V  *  \r 


(2.101) 


j(Kn)2(y)|  <  inn  p(»)  y'  |B 2(y')|  dy’ 


In  particular, 
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(2.102) 


U  Kl  <  p(«) 


>00  rCO 

y' |B2(y')  |dy'  +  p^y"1  y ,2  j  B2  (y ')  |  dy  ' 


+  p~1(«)y“1 


|Bj (y') | dy  *  +  p(“>) 


B2 (y ') |dy ' 


Hence  (2.85)  implies  that  K  is  contractive  in  X  for  yl  large  enough  and 
a  unique  solution  is  obtained  as  in  the  preceding  case.  Finally, 

(2.100)  and  (2.101)  together  with  (2.85)  imply  that  (2.96)  is  satisfied. 

The  existence  of  the  special  solutions  i^,  <p2  has  thus  been 
proved.  Their  linear  independence  follows  directly  from  (2.21),  (2.22). 
Proof  of  Corollary  2.7.  This  was  verified  by  (2.25). 
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53.  Spectral  Properties  of  A^. 

The  results  of  §2  are  used  in  this  section  to  derive  precise 
results  concerning  the  location  and  nature  of  the  spectrum  of  A^ .  The 
notations  a(A^) ,  afl(A^),  ac(A^)  and  ae(A^)  will  be  used  to  denote  the 
spectrum,  point  spectrum,  continuous  spectrum  and  essential  spectrum  of 
A^,  respectively.  The  definitions  of  [8,  Ch.  X]  will  be  used.  In 
particular,  ac(Ay)  is  a  closed  set  and  °e(Ay)  is  the  set  of  all  non¬ 
isolated  points  of  a(A^).  Note  that  the  properties  of  A^  described  by 
(1.19)  imply  that  cr(A^)  c  [c^y2,“). 

The  Point  Spectrum  of  A  .  Theorem  2.1  and  its  corollaries  imply 
the  following  three  lemmas  concerning  0„ (A^) . 

Lemma  3.1.  For  all  y  >  0, 

(3.1)  °o(AU)  C  [c^2*c2(“^2l  • 

Lemma  3.2.  For  all  y  >  0, 

(3.2)  a(A^)  n  [c2y2,c2(«)y2)  c  a„(AlJ). 

Moreover,  (^(A^)  is  either  a  finite  set  (possibly  empty)  or  a  countable 
set  with  unique  limit  point  cz(°°)y2. 

Lemma  3.3.  The  eigenvalues  of  A^  that  lie  in  the  interval 
[c2y2,c2(°°)y2)  are  all  simple. 

The  possibility  that  c2(°°)y2  €  a0(A^)  is  not  excluded  by  the 
hypotheses  (1.3),  (1.4)  alone.  Criteria  for  c2(°°)y2  £  a„(A^)  are  given 
below. 

It  will  be  convenient  to  use  a  notation  that  permits  a  unified 
discussion  of  the  cases  of  finite  and  infinite  point  spectra  og(A^). 
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The  number  of  eigenvalues  in  [c2y2,c2 (°°)y2)  will  be  denoted  by  N(y)  -  1. 
Thus  N(y)  is  an  extended  integer-valued  function  of  p  >  0  (1  <  N(y) 

<  +°°) .  The  eigenvalues  of  in  [c^y2,c2(°°)y2) ,  arranged  in  ascending 
order  will  be  denoted  by  A^(y)  •  1  <  j  <  N(y)  .  Thus 

(3.3)  c2y2  <  Ax(y)  <  A2(y)  <  *•*  <  c2(»)y2. 

The  corresponding  eigenfunctions  are 

(3.4)  i|>k(y,U)  *  ak(y)  <j>2(y,u,Xjc(y)) ,  k  *  1,2,... 

where  a^(y)  >  0  is  chosen  to  make  B^k(*»U)l  ■  1. 

The  Continuous  and  Essential  Spectra  of  A^.  Lemma  3.2  implies 

that  o  (A  )  C  [c2(°°)y2 ,°°)  .  Moreover,  a  (A  )  and  a  (A  )  are  closed  and 
6  u  c  p  c  y 

Oc(A^)  C  (A^)  (8,  Ch.  X].  The  characterization  of  these  sets  will  be 

completed  in  §5  by  showing  that  (c2(®)y2,°°)  C  a  (A  ) .  These  facts  imply 

c  p 

Theorem  3.4.  For  all  y  >  0, 

(3.5)  Cc(Ay)  *  °e(Ay)  “  [c2(°°)y2 ,“) . 

A  direct  proof  of  Theorem  3.4  can  be  given  by  using  the  special 
solutions  of  §2  and  a  criterion  of  Weyl;  see  [3,  p.  1435]. 

It  is  known  that  the  bottom  point  in  the  essential  spectrum  of  a 
Sturm-Liouville  operator  A  can  be  characterized  by  the  oscillation 
properties  of  the  solutions  of  A<J>  ■  A0  [3,  p.  1469].  For  the  operator 
Ay  the  characterization  is  described  by 

Corollary  3.5.  The  equation  A  $  *  A<£  is  oscillatory  (every 
real  solution  has  infinitely  many  zeros)  for  every  A  >  c2(°°)y2.  The 
equation  is  non-oscillatory  (every  real  solution  has  finitely  many  zeros) 
for  every  A  <  c2(°°)y2. 
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These  results  for  follow  directly  from  Theorem  2.1. 

The  Point  Spectrum  of  (continued) .  The  equation  A^  $  »  X$ 
may  or  may  not  be  oscillatory  for  X  -  c2(°°)y2.  This  property  is  shown 
below  to  provide  a  criterion  for  aQ (A^)  to  be  finite.  The  basic  tool  in 
establishing  such  criteria  is  the  classical  oscillation  theorem  of  Sturm. 
A  version  suitable  for  application  to  A^  may  be  formulated  as  follows. 

Let  I  *  (a,b)  be  an  arbitrary  interval  (-“  <  a  <  b  <  +*)  and 
consider  a  pair  of  equations 

(3.6)  Lj 4>  =  (Pjl(y)4>') '  +  Qj (y)4>  -  0,  j  -  1,2, 

where  Pj  (y)  and  Qj  (y)  are  defined  and  real  valued  for  almost  every  y  6  I, 
Pj (y)  >  0  for  almost  every  y  €  I  and  Pj ,  Qj  are  Lebesgue  integrable  on 
compact  subsets  of  I  (j  *  1,2).  A  solution  of  (3.6)  on  I  is  a  function 
<(>  €  AC(I)  such  that  Pj 1 *  €  AC(I)  and  (3.6)  holds  for  almost  all  y  €  I. 
Such  solutions  are  uniquely  determined  by  the  values  <Ky„)  “  c0, 

Pj 1  (y0)4>'  (y0)  *  ci  at  anY  Point  y0  e  Pairs  of  equations  (3.6)  such 
that 

(3.7)  Pi(y)  <  P2(y),  Qt(y)  <  Q2(y)  for  almost  all  y  e  I 

will  be  considered.  When  (3.7)  holds  the  operator  L2  is  said  to  be  a 
Sturm  majorant  of  operator  Lx,  and  the  operator  Lt  is  said  to  be  a  Sturm 
minorant  of  operator  L2,  on  I.  Sturm's  theorem  may  now  be  formulated 
as  follows. 

Theorem  3.6.  Let  4>j  (y)  $  0  be  solutions  of  L j j  *  0  on  I 
(j  ■  1,2)  and  assume  that  yt  and  y2  are  successive  zeros  of  4>2  (y)  in  I, 
with  yt  <  y2.  Moreover,  let  L2  be  a  Sturm  majorant  of  Lj  on  (yt>y2). 

Then  <$>2(y)  has  at  least  one  zero  in  [y j , y 2 ) .  In  addition,  if  either 


i . 
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Q3(y)  <  Q2(y)  or  px (y)  <  P2(y)  and  Q2(y)  Nona  subset  of  (yj,y2) 
having  positive  measure  then  02(y)  has  a  zero  in  (yx,y2). 

The  special  solution  <fc3(y,p,X)  is  real  valued,  tends  to  zero 
exponentially  when  y  -*■  and  has  finitely  many  zeros  when  X  <  c2(®)p2. 
Theorem  3.6  will  be  shown  to  imply 

Corollary  3.7.  If  X3  <  X2  <  c2(«)p2  then  <t>3(y,y,X2)  has  at 
least  as  many  zeros  as  <J>3  (y,p,Xx ) . 

It  will  be  convenient  following  [3,  p.  1473]  to  introduce  the 

sets 


(3.8)  Ik  =  Ik(P)  -  (X  |  <h(y,p,X)  has  exactly  k  zeros},  k  =  0,1,2, ••*. 

Note  that  by  Corollary  3.5  each  1^  C  (-»,c2(®)p2] .  The  point  c2(°°)p2 
may  or  may  not  be  in  one  of  the  sets  1^.  Corollary  3.7  implies  that 
each  1^  is  an  interval  and  Ik  lies  to  the  left  of  Ik+i  for  k  ■  0,1,2, •••. 
It  is  important  for  the  analysis  of  a„ (A^)  to  know  that  the  intervals 
Ifc  <j>  for  k  ■  1, 2,  •  •  •  ,N(p)  -  1.  This  is  a  corollary  of  the  following 
fundamental  oscillation  theorem. 


Theorem  3.8.  If  0O(A^)  f  <p  then  for  k  *  l,2,**»,N(p)  -  1  the 
eigenfunction  i^k(y,p)  has  precisely  k  -  1  zeros. 

Corollary  3.9.  If  a0(A^)  /  <p  then 


(3.9)  Ik  -  (Xk(p),Xk+1(p)],  k  -  0,1, •••.N(p)  -  2 

(where  X0(p)  =  -») .  Moreover,  if  N(p)  <  ®  then  (X^  (U)  »c2(°°)p2) 

C  ^(P)-!’ 

Corollary  3.10.  The  number  of  eigenvalues  that  satisfy 
Xk (p)  <  X  <  c2(<»)p2  is  equal  to  the  number  of  zeros  of  $3(y,p,X). 
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Criteria  for  the  Finiteness  of  (A^) .  The  principal  criterion 
for  o0(A^)  to  be  finite  is  described  by 

Theorem  3.11.  a0(A^)  is  finite  if  and  only  if  the  equation 

A^  $  =  c2(°°)y2<f>  is  non-oscillatory  on  R.  Hence  a„(A^)  is  infinite  if 
and  only  if  A^  <J>  ■  c2(°°)y2<j)  is  oscillatory  on  R. 

It  is  shown  below  that  Theorem  3.11  is  a  consequence  of  Theorem 
3.8  and  Stum's  comparison  theorem. 

Corollary  3.12.  If  c(»)  <  c(-<*>)  then  o0(A^)  is  finite  if  and 
only  if  <f>3(y,y,c2(°°)y2)  has  only  a  finite  number  of  zeros. 

Specific  criteria  for  the  finiteness  of  a0(A^)  will  now  be 
obtained  by  deriving  criteria  for  A^  $  =  c2(°°)y2<}>  to  be  non-oscillatory 
and  using  Theorem  3.11.  It  will  be  assumed  that  c(°°)  <  c(-°°)  so  that 
»  c2(«)y2<{>  is  non-oscillatory  in  neighborhoods  of  y  *  -<®.  Cases 
for  which  c(°°)  -  c(-°°)  may  be  treated  by  applying  non-oscillation 
criteria  at  both  y  »  °°  and  y  - 

A  criterion  for  a„(A^)  to  be  finite  is  provided  by  Theorem  2.6. 
For  under  the  conditions  of  the  theorem  A^  tp  =»  c2(°°)u2<J>  has  a  solution 
basis  <j>lf  <p2  satisfying  (2.21),  (2.22).  It  follows  that  the  equation  is 
non-oscillatory.  This  implies 

Theorem  3.13.  If  p(y) ,  c(y)  satisfy  (1.3),  (1.4),  c(°°)  <  c(-®) 

and 


(3.10) 


c(y)  -  c (°°)  |  dy  < 


00 


then  o0(Ay)  is  finite  for  every  y  >  0. 

Alternative  criteria  for  the  finiteness  of  a0(Ay)  can  be  derived 
by  constructing  Sturm  majorants  of  A  $  *  c2(°°)y2<J>  that  are  non-oscillatory 


rT7*p»"ll'#OW  - 


u  uiJJW! 
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and  using  Theorems  3.6  and  3.11.  Similarly,  criteria  for  o0(A^)  to  be 
infinite  can  be  derived  by  constructing  Sturm  minorants  that  are 
oscillatory.  Several  criteria  of  this  type  will  be  given. 

The  equation  4>  =*  can  be  written 

(3.11)  (p-1  (y) 4> ' ) '  +  P_1(y)(X  c"2(y)  -  u2)<J>  -  0. 

In  particular,  for  X  *  c2(°°)y2  one  has 

(3.12)  (p~1(y)<P') '  +  p-1(y)U2(c2(»)c"2(y)  -  l)<j>  -  0. 

The  first  factor  p  1 (y)  in  (3.12)  is  unimportant  for  the  oscillation 
properties  of  the  equation.  Replacing  it  by  gives  the  majorant 

(3.13)  <j>"  +  PM  P-1  (y)U2  (c2(°°)c~2  (y)  -  1)*  «  0. 

Each  non-oscillatory  Sturm  majorant  of  (3.13)  gives  a  criterion  for  the 
finiteness  of  afl(A  ).  Since  solutions  of  (3.13)  are  non-oscillatory  on 
any  interval  (-°°,yQ)  when  cO”)  <  c(-°°),  it  is  enough  to  construct 
majorants  of  (3.13)  on  intervals  (y0,°°).  An  obvious  non-oscillatory 
majorant  for  (3.13)  is  4>"  =  0.  Thus  ct0(A^)  is  finite  for  every  u  >  0  if 
there  is  a  y0  such  that  c2(°°)c  2(y)  -  1  <  0  for  all  y  >  y„;  that  is, 

(3.14)  c(y)  >  c(°°)  for  y  >  y„. 

This  means  the  graph  of  c  =  c(y)  lies  above  or  on  the  limit  line  c  =  c(®) 
in  a  neighborhood  of  y  ■  00 .  Weaker  hypotheses  that  include  this  case 
can  be  derived  by  comparing  (3.13)  with 

(3.15)  4>"  +  a  y~2  <J>  -  0 

which  is  oscillatory  on  (yo»°°)  if  a  >  1/4  and  non-oscillatory  if  a  <  1/4. 


49 


Oscillation  theorems  based  on  (3.15)  were  first  given  by  A.  Kneser  (9); 
see  [3,  p.  1463].  Comparison  of  (3.12)  with  (3.15)  gives 
Theorem  3.14.  If  c(°°)  <  c (-°°)  and 

(3.16)  lim  sup  y2(c  2(y)  -  c  2(®))  <  0 

y-*oo 

then  o^A^)  is  finite  for  all  p  >  0.  Conversely,  if 

(3.17)  lim  inf  y2(c  2(y)  -  c  2(°°))  >  0 

y-*co 

then  there  exists  a  p0  >  0  such  that  O0(A^)  is  infinite  for  every  p  >  p0. 

Mote  that  the  criterion  (3.16)  includes  (3.14)  as  a  special  case. 
Note  also  that  sufficient  conditions  for  (3.16)  or  (3.17)  to  hold  are  the 
existence  of  constants  y0,  K  and  e  >  0  such  that 

(3.18)  c(y)  >  c(«)  -  K  y"2~e  for  y  >  y0 
or 

(3.19)  c(y)  <  c(°°)  -  K  y”2  for  y  <  y0. 


respectively.  In  particular,  a0(A^)  is  finite  for  all  u  >  0  if  c(y) 
approaches  c(°°)  from  below  sufficiently  rapidly. 


Criteria  that  Guarantee  a  (A  )  j4  <f>.  Such  criteria  may  be  derived 


by  constructing  Sturm  minorants  for  A  $  ■  c2(°°)u2<|>  whose  solutions  have 
zeros.  If  the  minorant  has  solutions  with  infinitely  many  zeros  then 


CT0(A^)  is  infinite.  If  the  minorant  has  a  solution  with  finitely  many 
zeros  then  it  can  be  shown  that  <J>3(y,y,c2(°°)u2)  has  at  least  as  many 


zeros  and  one  may  use  the  following  refinement  of  Theorem  3.11. 


Theorem  3.15.  If  A  $  ■  c2(«,)y2<J>  has  a  solution  having  a  finite 
number  k  of  zeros  on  R  then  the  part  of  a )  below  c2(«>)p2  is  finite  and 
has  at  least  k  -  1  and  at  most  k  +  2  points. 


m 


To  apply  the  method  in  cases  where  o0(A^)  is  finite  consider 
first  the  case  p(y)  *  const,  so  that  (3.12)  becomes 

(3.20)  0>"  +  u2(c2(°°)c“2(y)  -  l)<fr  -  0 

Note  that  c2(°°)c  2(y)  -  1  >  c2(“)c02(y)  -  1  for  all  y  6  R  if  and  only  if 

(3.21)  c(y)  <  c0(y)  for  all  y  €  R. 

If  c<j(y)  can  be  chosen  in  such  a  way  that 

(3.22)  0>"  +  u2(c2(co)c~2(y)  -  1)0>  *  0 

has  a  solution  on  R  with  k  zeros  then  ao(A^)  will  have  at  least  k  -  1 
points  by  Theorem  3.15.  In  this  way  one  can  prove 

Theorem  3.16.  Let  p(y)  -  const,  for  all  y  e  R  and  assume  that 
there  is  a  constant  cQ  >  cffl  and  an  interval  I  *  [a,b]  with  b  >  a  such 
that 


(3.23)  c(y)  <  c0  <  c(“0  <  c(-°°)  for  all  y  €  I. 

Then  O0(A,  )  ^  0>  for  all  sufficiently  large  y.  In  fact,  N(y)  -*■  when 

p  -►  0°. 

Theorem  3.16  can  be  proved  by  comparing  c(y)  with  a  suitable 
piece-wise  constant  function  c0(y)  that  satisfies  (3.21).  An  analogue 
of  Theorem  3.16  can  be  proved  in  the  general  case  where  p(y)  +  coft^t.  by 
making  the  change  of  variable  y  +  n  in  (3.12),  where 


(3.24) 


n 


y 

p(y')  dy'. 


The  details,  which  are  elementary  but  lengthy,  are  omitted. 


This  completes  the  formulation  of  the  results  of  §3  and  the 


proofs  will  now  be  given.  Note  that  Lemma  3.1  is  an  immediate  conse¬ 
quence  of  Theorem  2.1  which  implies  that  for  A  >  c2(“)p2  the  equation 
A^  <J>  =  A0  has  no  solutions  in  3C(R) . 

Proof  of  Lemma  3.2.  The  resolvent  of  is  an  integral  operator 
in  3f(R)  [3,  XIII. 3] 


(3.25) 


(Ay  -  O'1  f(y) 


I  G  (y.y’.O  f(y')  c_2(y')  p-1(y')  dy'. 
JR  M 


Gy(yty'«?)>  the  Green's  function  of  A  ,  is  known  to  have  the  form 
[3,  p.  1329] 


(3.26) 


Gy(y,y',S) 


tO-J 


I  ^(y)  ^(y') »  y  <  y'» 
^(y)  G-o^y').  y  >  y’» 


where  0^  and  are  non-trivial  solutions  of  A^  0  =  £0  that  are  in 
L2 (O,00)  and  L2(-®,0),  respectively.  Thus  for  £  e  L(c(°°)y)  C  L(c(-®)y), 
<P.„m  $2  *  0_o „  *  $3  and  one  has 


(3.27) 


Gy(y, y'.O  -  [0203J 


tj>3 (y ,u, ?)  <t>2 (y *  ,u,c) ,  y  <  y\ 
02(y»y»c)  4,3(y,»y»?) »  y  >  y’. 


It  follows  from  Corollary  2.3  that  G^(y,y',c)  is  meromorphic  in 
L(c(°°)]j)  with  poles  at  the  zeros  of 

(3-28)  F(p,s)  -  [02(*,p,c)  03(*  ,p,C)] . 


As  remarked  in  §1,  these  are  precisely  the  eigenvalues  of  A^  that  are 
less  than  c2(°°)ij2.  Their  only  possible  limit  point  is  c2(°°)vi2  since 
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F(U,0  is  analytic  in  L(c(°°)u)  by  Corollary  2.3.  These  results  imply  the 
two  statements  of  Lemma  3.2. 

Proof  of  Lemma  3.3.  This  follows  from  Theorem  2.1  which  implies 
that  Ap  <{>  =  X<{>  always  has  at  least  one  solution  that  is  not  in  TC(R) . 

Proof  of  Theorem  3.4.  It  was  remarked  above  that  ac(Ay)  C  Oe(A^) 

C  [ c2  (°°)U2 ,°°) .  Hence  to  prove  (3.5)  it  is  enough  to  show  that 

(c2  (°°)p2 ,°°)  C  ac(A^).  This  may  be  done  by  constructing  a  characteristic 

sequence  for  each  X  €  (c2  (°°)u2 ,°°) ;  i.e.,  a  bounded  sequence  {^(y)}  in 

3f(R)  such  that  each  <J>n  €  D(A^)  and  (A^  -  A ) 4>n  0  in  Jf(R)  but  { 4>n}  has 

no  convergent  subsequences.  Indeed,  a  suitable  sequence  has  the  form 

<i>n(y)  =  Cn(y)  <My>L,X)  where  5n  6  D(A0),  £n(y)  =  1  for  |y|  <  n, 

supp  C  [-n-l,n+l]  and  £^(y)  and  (p  1(y)C^(y))'  are  bounded  for  all  y 

and  n.  Such  a  sequence  {t^}  can  be  constructed  but  the  details  are 

lengthy.  They  will  not  be  given  here  since  the  inclusion  (c2 (°°)y2 ,°°) 

C  0  (A  )  is  proved  in  §5. 
c  y 

Proof  of  Corollary  3.5.  For  X  t  c2(°°)u2  every  solution  of 
A^  <j>  =  X<J>  is  a  linear  combination  of  <j)j(y,)j,X)  and  <J)2(y,u,X)  (Corollary 
2.7).  It  follows  from  Theorem  2.1  that  every  real  solution  with 
X  >  c2(“>)u2  has  infinitely  many  zeros  in  any  interval  (y0,°°).  On  the 
other  hand  for  X  <  c2(=°')y2  Theorem  2.1  implies  that  every  real  solution 
of  A  <J>  ■  X0  is  either  exponentially  large  or  exponentially  small  for 
y  -*•  ±°°.  In  every  case  $(y)  has  constant  sign  outside  of  some  interval 
t-y0,y  ]  and  hence  can  have  only  finitely  many  zeros. 

Proof  of  Theorem  3.6.  Results  equivalent  to  Theorem  3.6  are 
proved  in  [6,  Ch.  XI]  under  the  additional  hypothesis  that  the  P^  and 
are  continuous.  The  same  method  will  be  shown  to  be  applicable  under 
the  hypotheses  of  Theorem  3.6.  The  method  is  to  study  the  phase 
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plane  curves 

(3.29)  (C.n)  =  (PT1(y)<J)J!(y),<J>j(y)),  y  e  i,  j  =  1,2, 

defined  by  the  solutions  4>j(y)  and  to  transform  to  polar  coordinates 
(Prufer  transformation).  Thus  (3.29)  can  be  written 

(3.30)  (C.n)  =  (rj(y)  cos  9j(y)»rj(y)  sin  6j(y)),  y  €  I,  j  =  1,2. 

Moreover,  the  curves  (3.29)  cannot  pass  through  the  origin  because 
<(> .  ( y)  ?  0.  Thus  Tj  (y)  >  0  and  0^  (y)  is  uniquely  defined  by  continuity 
and  its  value  at  the  point  yx  E  I.  Finally,  0,  E  AC(I)  and  (3.6) 
implies  that  0^  is  a  solution  of  the  first  order  equation 

(3.31)  0j(y)  =  Pj(y)  cos20j (y)  +  Qj(y)  sin20j (y) ,  y  £  I. 

To  prove  the  first  statement  of  Theorem  3.6  note  that  one  can 
assume  without  loss  of  generality  that  <*>i  (y)  >  0  for  yt  <  y  <  y2  and 
4>2(y1)  >  0.  Thus  0j  (y)  (j  =  1,2)  may  be  defined  as  the  unique  solutions 
of  (3.31)  such  that  01(y1)  =  0  and  0  <  02(yt)  <  IT-  It  follows  that 
0  <  0i  (y)  <  7T  for  yt  <  y  <  y2  and  81(y2)  =  ir.  It  must  be  shown  that 
4>2(y)  has  a  zero  in  [y15y2).  If  <j>2(yi)  =  0  there  is  nothing  to  prove. 

If  <f)2(yl)  >  0  then  0  <  d2(yl)  <  it  and  it  follows  from  (3.31)  and  (3.7) 
that  02 (y)  >  0j(y)  for  all  y  >  yj  (see  [6,  p.  335]).  In  particular, 
02(y2)  >  0!  (y2)  =  7T  whence  by  continuity  02(yQ)  =  and  therefore 
$2(y0)  a  0  for  some  yg  e  (yj,y2). 

To  prove  the  second  statement  of  Theorem  3.6  it  is  only  necessary 
to  remark  that  if  Qt (y)  <  Q2 (y)  or  Pj (y)  <  P2 (y)  and  Q2 (y)  #  0  on  a 
subset  of  (yj,y2)  having  positive  measure  then  02(y2)  >  tt  even  if 
02(yt)  =*  0;  see  [6,  p.  335].  This  completes  the  proof. 


Proof  of  Corollary  3.7.  The  function  <{>3(y,M,A)  is  a  solution  of 
equation  (3.11).  Thus  if  4>j  (y)  3  <J>3  (y ,  U ,  )  ,  j  =  1,2,  then 

(3.32)  (P_1  (y)4>j ) '  +  p'1(y)(Xjc"2(y)  -  u2)^  =  0,  j  =  1,2. 

These  equations  have  the  form  (3.6)  with  Pj (y)  =  p(y)  and  Qj (y) 

3  P_1 (y) (^jc_Z (y)  ~  M2).  Hence  P3(y)  =  P2 (y)  and  Q3  (y)  <  Q2 (y)  for  all 
y  G  R,  since  p(y)  and  c(y)  are  always  positive,  and  the  second  part  of 
Theorem  3.6  is  applicable.  It  follows  that  if  yt  <y2  <  •••  <  yk  are 
the  zeros  of  ^(y)  =  4>3(y,U»^i)  then  02(y)  =  <|>3(y,u»X2)  has  k  -  1  zeros 
in  the  interval  (yx,yk).  Hence  it  will  be  enough  to  show  that  4>3(y,u»X2) 
also  has  a  zero  in  (-“, y3].  To  verify  this  apply  Lagrange's  formula 
(2.23)  to  4>j(y)  and  <J>2(y)  in  (— °®, y 3 ] .  This  is  possible  because  4>3(y) 
and  $2(y)  are  exponentially  small  at  y  =  The  result  can  be  written 


(X2 


fyi 

-  V 

—00 


(y)  $2(y)  c-2(y)  p-1(y)  dy 


(3.33) 


{<!>!  <J>2  -  <t>2  (Jjj }  C-2  p-1  dy 


=  <J>2(yj)  (p'^yj)  4>i (y i ) } 


since  (l)1(y1)  3  0  and  <f>j  and  $2  vanish  at  y  *  Now  suppose  that  <P2(y) 

has  no  zero  in  (— °°, y t ] .  Then  02  (yt)  >  0  because  <}>2(y)  >  0  near  y  =  -» 
by  Theorem  2.1.  Moreover,  $3(y)  >  0  for  -<*>  <  y  <  y3  and 
P  X(y3)  ^'(yj)  <  0  because  yj  is  the  first  zero  of  <J>j .  Thus  the  right 


hand  side  of  (3.33)  is  negative.  But  the  left  hand  side  is  clearly 
positive.  This  contradiction  completes  the  proof. 
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Proof  of  Theorem  3.8.  For  regular  Sturm-Liouville  operators 
A<t>  =  -4>"  +  q(y)4>  on  finite  intervals  the  oscillation  theorem  goes  back 
to  Sturm.  For  singular  operators  in  L2(0,°°)  such  that  q(y)  -*■  +°°  when 
y  -*•  +°°  it  was  first  proved  by  H.  Weyl.  More  recently  the  result  was 
proved  by  B.  M.  Levitan  and  I.  S.  Sargsjan  [11,  p.  201]  under  more 
general  conditions  on  Q(y)  that  guarantee  that  the  solutions  of  A0  *  \<p 
are  non-oscillatory  on  0  <  y  <  00  for  all  X  6  R  (and  hence  a(A)  is 
discrete).  It  will  be  shown  here  that  the  method  of  Levitan  and  Sargsjan 
is  applicable  to  the  case  of  Theorem  3.8. 

The  method  of  [11]  is  to  regard  the  Sturm-Liouville  problem  for 
A  on  0  <  y  <  ®  as  a  limit  of  regular  problems  for  Aon0<y<b<°°  and 
to  study  the  behavior  of  the  eigenvalues  and  eigenfunctions  as  b  +  ®. 


Here  the  operator  in  JC(R)  will  be  regarded  as  a  limit  of  the  regular 

—  2  —  1 

operator  in  3C(a,b)  =  L2(a,b;c  (y)  p  (y)  dy) ,  -°°  <  a  <  b  <  °°,  defined 

by  A^  and  the  boundary  conditions  p(a)  =  <{> (b)  *  0.  The  corresponding 

operator  in  3C(a,b)  will  be  denoted  by  A  .  The  limit  a  -«•  will  be 

p ,  a ,  d 

studied  first. 

The  operator  A^  is  more  general  than  the  operator  studied  in  [11] 
However,  examination  of  the  proofs  in  [11]  reveals  that  nothing  is  used 
but  the  Sturm  comparison  theorem,  the  convergence  of  the  eigenvalues 
when  b  -*■  00 ,  the  continuity  and  asymptotic  properties  of  the  solution 
<Ky»X)  of  A<p  =  A<j>  that  satisfies  the  boundary  condition  at  y  *  0  and  the 
non-oscillatory  character  of  A$  =  X<j>  in  a  X-interval  containing  the 
point  spectrum.  All  of  these  properties  have  been  established  for  A^. 

The  solution  of  A^  <j»  ■  X<j>  that  satisfies  <}>(b)  *  0,  p  '  (b)  p'(b) 

*  1  will  be  denoted  by  0^(y,X).  For  X  +  c2(-°°)y2,  ^(y.X)  is  a  linear 
combination  of  $3(y,p,X)  and  (^(y.p.X)  and  hence  has  the  regularity 
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properties  of  Corollaries  2.2  and  2.3.  The  eigenvalues  of  A  are 

U  *  &*  D 

the  roots  of  the  equation  ( a , A )  =  0.  They  will  be  denoted  by 

X  ,  k  =  1,2, ••*,  with  the  convention  that  X,  ,  <  X,  ,  .  The 

D  1C)  cl)  D  K.t  1  D 

corresponding  eigenfunctions  ik  K(y)  “  4>,  (y , X.  .  )  have  precisely 

k  -  1  zeros  by  the  classical  oscillation  theorem.  For  tho.  din  of 

operators  considered  here  this  result  can  be  proved  by  the  method  of 

(k) 

[11,  p.  17].  The  zeros  of  i ik  ,  (y)  will  be  denoted  by  y.  , 

*£)  cl  )  D  J  j  3  y  D 

1  <  j  <  k  -  L 

The  operator  in  3C(-«>,b)  defined  by  and  the  boundary  condition 
<(>(b)  =  0  will  be  denoted  by  A^  The  methods  used  to  study  a(A^)  above 
can  be  used  to  show  that  a(A^  n  [c^y2 , c2 (-°°)y2)  (C  o0  (A^  ^))  is 

finite  or  countably  infinite  with  unique  limit  point  c2(-°°)y2.  The 
number  of  eigenvalues  in  [c^y2 ,  c2  (-°°)y2)  will  be  denoted  by  N(y,b)  -  1 
(<  +oo)  in  analogy  with  the  notation  for  A^,  and  the  eigenvalues  will  be 
denoted  by  Xfc  (X^  <  X^+1  .  The  eigenfunctions  for  A^  are 

\,b(y}  3  *b(y-xk,b)* 

The  proof  of  the  oscillation  theorem  for  A  ,  by  the  method  of 

y ,  d 

[11]  will  now  be  outlined.  First, 


(3.34)  lim  X,  ,  ■  X,  ,  for  1  <  k  <  N(y,b). 

a-*— oo  *’a’D  1S’D 

This  follows,  for  example,  from  the  convergence  of  the  Green's  functions. 
It  follows  that 


(3.35)  lim  ik  (y)  *  i k  ,  (y)  for  -oo  <  y  <  b, 

a->— 00  K*D 

uniformly  on  bounded  subsets  of  (-°°,b].  The  proof  of  the  oscillation 
theorem  given  in  [11]  is  based  on  the  following  three  lemmas. 


57 


Lemma  3.17.  For  each  k  =  1, 2, • * • ,N(y,b)  -  1  and  each  fixed 
b  S  R  one  has 


(3.36) 


(k) 

SUP  yk-i  a  b  • 
a<b  *  1,a,D 


Lemma  3.18.  For  each  k  *  1, 2, • • • ,N(y ,b)  -  1,  each 
j  *  l,2,***,k  -  1  and  each  fixed  b  €  R  one  has 


(3.37) 


inf 

a<b 


(  (k) 
j+i ,a,b 


,<k>  J  >  0- 

J ,a,bj 


Lemma  3.19.  For  each  k  *  1, 2, • • * ,N(y ,b)  -  1  and  each  fixed 
b  €  R  one  has 


(3.38) 


(Xk+l,a,b  "  Xk,a,b)  >  °* 


The  proofs  of  these  lemmas  and  the  oscillation  theorem  for  A 


U,b 


are  the  same  as  those  given  in  [11,  pp.  202-4]  and  will  not  be  repeated 
here. 


The  proof  of  Theorem  3.8  may  now  be  completed  by  regarding  A  as 


a  limit  of  A  .  for  b  -*•  00  and  repeating  the  argument  given  above.  The 


solution  of  A^  <{>  ■  A<J>  that  satisfies  the  condition  of  square  integrability 


at  y  ■  is  d>3  (y » u , A )  and  is  non-oscillatory  for  all  y  €  R  when 
\  <  c2(°°)y2  (Corollary  3.5).  The  remainder  of  the  proof  follows  as 
before. 

Proof  of  Corollary  3.9.  Theorem  3.8  implies  that  ^+1 (h)  6  1^ 
for  k  ■  0,l,2,***,N(y)-2.  Moreover,  a  continuity  argument  based  on 
Corollary  2.2  shows  that  the  intervals  1^  have  the  form  1^  ■  (ajc»ajc+1l 
where  a0  <  <  •••  <  (see  [3,  p.  1475]).  Thus  to  prove  (3.9) 

it  will  suffice  to  prove  that 
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(3.39) 


F(y,ak)  -  ($2(*,y,ak)  4>3  (•  ,y,ak)  ] 


for  k  *  l,2,***,N(y)  -  1.  This  proof  will  be  based  on  the  following  two 
lemmas . 

Lemma  3.20.  Let  X0  <  c2(-°°)y2  and  let  y0  be  a  zero  of 
$3 (y.U»X0) •  Then  to  each  sufficiently  small  e  >  0  there  corresponds  a 
6  >  0  such  that  for  | X  —  X 0  J  <  5  the  function  $3(y,y,X)  has  exactly  one 
zero  in  the  interval  |y  -  y  |  <  e. 

This  result  follows  from  Corollary  2.2  and  the  fact  that 
P_1(y)  4> '  (y)  cannot  vanish  at  a  zero  of  a  non-trivial  solution  of 
<p  ■  X<{>.  For  a  proof  see  [11,  p.  16]. 

For  X  €  i  let  y3(X)  <  y2(X)  <  •••  <  yk(X)  denote  the  zeros  of 
4>3(y»y»X).  Then  each  y^(X)  is  uniquely  defined  for  X  e  Uk>j  \  and 
one  has 

Lemma  3.21.  Each  of  the  functions  y^(X)  is  continuous  and 
strictly  monotone  increasing. 

The  continuity  follows  immediately  from  Lemma  3.20.  The  strict 
monotonicity  follows  from  the  proof  of  Corollary  3.7. 

Proof  of  Corollary  3.9  (concluded).  (3.39)  will  be  proved  by 
contradiction.  Assume  that  F(y,ak)  f  0  and  note  that  for  X  <  c2(°°)y2 
one  has 


(3.40)  4>3 (y »y»X)  -  c(y,X)  ^(y.y.X)  +  c’(y,X)  4>2(y,y,X), 
by  Corollary  2.7.  Moreover,  Theorem  2.1  implies  that 

(3.41)  c(y,X)  -  P(“)  [<i>24>3 1/2iq+(y  ,X)  -  p(«)  F(y,X)/2iq+(y,X)  . 

Thus  c(y,ak)  i  0  and  by  continuity  (Corollary  2.2)  there  is  an  interval 
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|X  -  a^J  <  6  in  which  c(y,X)  ^  0.  It  follows  from  (3.40)  and  the 
uniformity  of  the  asymptotic  estimates  of  Theorem  2.1  (Corollary  2.4) 
that  there  is  an  M  >  0  such  that 


(3.42)  |  <J)3  (y,w, X)  |  >  1  for  all  y  >  M  and  |X  -  a^|  <  6. 


Note  that  Lemma  3.21  implies 

(3.43)  lim  y.(X)  *  y.(a.),  j  -  l,2,**»,k. 

***il 

Now  consider  y^+i (X)  which  is  defined  for  X  a^.  For  a^  <  X  <  a^  +  6 
one  has  y  (X)  <  M  by  (3.42).  It  follows  that  the  limit 


(3.44) 


lim  yfc+i(X) 
Wk 


exists  and  y  >  But  <j>3(y,y,ak)  *  0  by  continuity  and  hence 

y  »  y^C^) .  But  this  implies  that  for  6  >  0  small  enough  and 
ak  <  X  <  ak  +  ^  every  neighborhood  |y  -  I  <  e  contains  two  zeros 

of  4>3(y,U,X)  in  contradiction  to  Lemma  3.20.  This  completes  the  proof 
of  (3.9). 


The  last  statement  of  Corollary  3.9  follows  from  the  proof  of 

Theorem  3.8.  Indeed,  if  N(y)  <  00  and  N(y)  <  N(y,b)  then  N(y)  <  N(y,b) 

and  X  ,  .  is  defined.  In  this  case 
N(y) ,b 


(3.45) 


since  otherwise  X^^  -*■  X0  <  c2(°°)y2  for  some  subsequence  (b(k)}, 

which  would  imply  that  X0  was  an  additional  eigenvalue  of  A^.  If 
N (y )  -  N(y,b)  <  00  the  same  argument  can  be  applied  to  the  operator  A  , 
Proof  of  Corollary  3.10.  This  follows  immediately  from 


Corollary  3.9. 
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Proof  of  Theorem  3.11.  It  will  suffice  to  prove  the  second 
statement  of  the  theorem.  To  this  end  let  ^(y.p.A)  be  the  solution  of 
A  <j>  ■  X<{>  that  satisfies  ^(b.y,!)  *  0,  p-1(b)  <Ji£(b,u,X)  »  1.  Then 
0^(y,U,X)  and  p  1 (y)  <^(y,y,A)  are  continuous  functions  of  (y,X)  e  R2. 

Now  assume  that  A  $  ■  c2(°°)y2<t>  is  oscillatory.  Then  $jJ(y,u,c2(00)u2) 
has  infinitely  many  zeros.  It  follows  by  the  method  used  to  prove 
Lemma  3.20  that  the  number  of  zeros  of  4>^(y,y,X)  tends  to  infinity  as 
\  -*■  c2(°°)y2.  But  then  the  same  is  true  of  <J>3(y,y,X),  by  Theorem  3.6, 
and  it  follows  from  Corollary  3.10  that  aQ(A^)  is  infinite.  To  prove 
the  converse  note  that  if  a0(A^)  is  infinite  then  Theorem  3.6,  applied 
to  the  kth  eigenfunction  and  any  solution  of  <j>  =  c2(°°)y2<{>  implies 
that  $  has  k  -  2  zeros.  Since  k  is  arbitrary  it  follows  that 
A  $  ■  c2(°°)u2<t>  is  oscillatory. 

Proof  of  Corollary  3.12.  This  follows  immediately  from  Theorem 
3.11.  The  hypothesis  c(°°)  <  c(-°°)  is  needed  only  to  ensure  that 
<S>3(y.y,c2(«>)y2)  is  defined. 

Proof  of  Theorem  3.13.  This  follows  immediately  from  Corollary 
3.12  and  Theorem  2.6. 

Proof  of  Theorem  3.14.  To  prove  the  first  half  of  the  theorem 
it  will  be  shown  that  condition  (3.16)  implies  the  existence  of  a  non- 
oscillatory  majorant  for  equation  (3.13)  for  every  y  >  0.  This  implies 
that  (3.12),  i.e.,  A^  <j>  *  c2(°°)y2$,  is  non-oscillatory  for  every  y  >  0 
and  the  finiteness  of  a0(A^)  follows  from  Corollary  3.12. 

To  construct  a  majorant  for  (3.13)  note  that  (3.16)  implies  that 
for  every  e  >  0  there  is  a  y,  ■  y0(e)  such  that 

(3.46)  y2(c2(°°)c  2 (y)  -  1)+  <  e  for  all  y  >  y0(e), 
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where  ct+  ■  Max  (a,0) .  It  follows  that  for  every  y  >  0  there  is  a 
y0  “  y0(M)  such  that 

pM  P_1(y)  y2  (c2(°°)c-2 (y)  -  1)  <  Pjj  P_I(y)  y2  (c2(“)c-2(y)  -  1)+ 

(3.47)  <  PM  p”1  y2(c2(°°)c-2(y)  -  1)+ 

<  l/4y2  for  all  y  >  y0 (y) . 


Hence  for  any  y  >  0  one  has 

(3.48)  PM  p-1(y)  y2(c2(°°)c-2  (y)  -  1)  <  l/4y2  for  all  y  >  y0(y). 

It  follows  on  comparing  (3.13)  with  (3.15)  with  a  =  1/4  that  (3.13)  is 
non-oscillatory  on  y0(y)  <  y  <  °°.  It  is  non-os cilia to ry  on  -»  <  y  <  y„(y) 
for  any  y  >  0  because  c(°°)  <  c(-°°).  This  proves  the  first  half  of 
Theorem  3.14. 

To  prove  the  second  half  it  will  be  shown  that  (3.17)  implies 
the  existence  of  a  y0  >  0  such  that  <(>  ■  c2(°°)y2<j>  is  oscillatory  for 

every  y  >  y0.  The  result  then  follows  from  Theorem  3.11.  To  this  end 
note  that  if  £  satisfies 

(3.49)  0  <  £  <  lira  inf  y2(c2(°°)c  2  (y)  -  1) 

y~*oo 

then  there  is  a  yt  >  y0(s)  such  that 

(3.50)  y2(c2(“)c_2(y)  -  1)  >  £  for  all  y  >  y0(£). 

In  particular,  given  any  a  >  1/4  there  is  a  y„  >  0  such  that 

(3.51)  0  <  pM  p~!a/y2  <  lira  inf  y2 (c2 (“) c-2 (y)  -  1). 

y^oo 


H 

ii 
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It  follows  that  there  is  a  yQ  ■  yQ (a)  such  that 

(3.52)  y2(c2(“)c-2  (y)  -  1)  >  pM  p~la/yj  for  all  y  >  y0(a). 

This  implies  that 

(3.53)  Pm  p^1  y2(c2(°°)c  2(y)  -  1)  >  a/y2  for  all  y  >  y0(a). 

Hence,  comparison  of 

(3.54)  <J>"  +  pm  p”1  y2(c2(°°) c~2(y)  -  1)4)  =  0 

and  (3.15)  with  a  >  1/4  implies  that  (3.54)  is  oscillatory.  But  (3.54) 
is  a  Sturm  minorant  of  (3.12);  i.e.,  <j>  *  c2(°°)p2<J>,  provided  y  >  y0 . 

Hence  the  latter  is  oscillatory  for  all  y  >  yQ . 

Proof  of  Theorem  3.15.  This  result  is  proved  in  [3,  p.  1481]  for 
Sturm-Liouville  operators  with  smooth  coefficients.  The  proof  is  based 
on  the  oscillation  theorem  (Theorem  3.8),  Sturm's  comparison  theorem  and 
the  continuous  dependence  of  the  zeros  of  solutions  of  if  =  A<f>  on  X 
(Lemma  3.20).  Hence  it  extends  immediately  to  the  operator  A^. 


§4.  Generalized  Eigenfunctions  of  A^. 

The  eigenfunctions  i|>k(y,y)  corresponding  to  the  point  spectrum 
of  A^  were  constructed  in  the  preceding  section.  In  this  section  the 
special  solutions  <fK  (j  =  1,2, 3,4)  of  §2  are  used  to  construct  gener¬ 
alized  eigenfunctions  of  A^  corresponding  to  the  points  of  the  continuous 
spectrum.  These  functions  will  be  used  in  §5  to  construct  the  spectral 
family  (11^(1)}  of  A^  and  to  prove  that  ac(A^)  *  ae(A^)  ”  [ c2 (“)y2 ,») . 

To  construct  the  generalized  eigenfunctions  ip. (y ,y, A)  ,  ^+(y>y»A) 
described  in  §1  recall  that  the  special  solutions  <jK(y,y,A)  are  defined 
for  all  real  A  £  c2(±°°)y2  and  the  pairs  (J>j ,  <j>2  and  $3,  <j>4  are  solution 
bases  for  A^  <J>  =  A<{>  (Corollary  2.7).  It  follows  that 


(4.1) 


<*>.  -  c.3<j>3  +  Cji>04,  j  =  1,2, 


=  Cjl*l  +  Cj2^z’  j  = 


The  coefficients  c.,  =  c.,  (y,A)  can  be  calculated  by  means  of  the 

J  k.  j  K. 

bracket  operation 


(4-2)  [^$k](y,A)  =  (•  ,y ,A)  4>k(*,y,A)] 

of  Lagrange's  formula.  Indeed,  by  forming  the  brackets  of  equations 
(4.1)  with  $4,  $3,  02  and  <f>j  in  succession  and  using  the  asymptotic 
forms  of  Theorem  2.1  one  finds 


(4.3) 


(~2iq_)  Cj  3  =  p(-~)  [<*>j4>4  ] 


j  *  1.2, 


(2iq_)Cj4  =  p(-«)  [ 4> j 4> 3  1 
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(4.3  cont.) 


(“21^) Cj  l  =  p(<=°)[0j02] 


(21^)  Cj  2  *  p(°°)  ] 


j  =  3,4. 


In  particular.  Corollary  2.2  implies  chat  each  c.,  (y,X)  is  a  continuous 

J 

function  for  X  t  c2(±°°)y2.  These  relations  will  be  used  to  determine 
the  generalized  eigenfunctions  of  A  .  The  notation 


(4.4) 


A  =  A(y)  =■  {X  |  c2(-«)y2  <  X}, 

AQ  *  Afl(y)  -  {X  |  c2(“)y2  <  X  <  c2(-®)y2} 


will  be  used.  Note  that  AQ  #  (})  only  if  c(°°)  <  c(-°°). 

The  Spectral  Interval  A.  The  generalized  eigenfunctions  of 
are  the  bounded  solutions  of  the  differential  equation  A^  <f>  =  Xc}) .  For 
X  €  A,  Theorem  2.1  and  the  relations  (4.1)  imply  that  all  the  solutions 
are  bounded.  It  will  be  shown  that  the  functions 


(4.5) 


^+(y»y,X)  =  a+(y,X)  <}>„  (y,y  ,X) 
'I'.Cy.MtX)  =  a_(y,X)  0j(y,y,X) 


have  the  asymptotic  forms  described  in  §1.  The  completeness  in 
TI^ (A)  3C(R)  of  these  functions  will  be  proved  in  15.  The  pair  (<f>2,i)>3), 
which  provides  an  alternative  basis,  will  not  be  treated  explicitly  here. 
It  may  be  shown  to  correspond  to  the  second  family  {<j>_(y,p,q)  |  q  >  0} 
described  at  the  end  of  §1. 

It  follows  from  (4.5),  (4.1)  and  Theorem  2.1  that  the  asymptotic 


behavior  of  is  given  by 


1 


I 


c41  eiq+y  +  c  ,  e-lq+y  +  o(l),  y  -*•  +«>, 


(4.6)  ^(y.P.X)  =  a+ 


riq-y  +  o( 


The  equivalence  of  (4.6)  and  Che  asympcotic  form  (1.34)  of  §1  follows 


Lemma  4.1.  For  all  p  >  0  and  X  €  A  the  coefficients  c4  (y,X), 
c42(p,X)  satisfy 

(4*7)  p  1  (°°)  <4  |c42i2  *  p  1(°°)  q+  lc41l2  +  p-1(-*)  q_. 

In  particular,  c42(y,X)  f  0. 

The  proofs  of  Lemma  4.1  and  subsequent  lemmas  are  given  at  the 
end  of  the  section. 

The  asymptotic  forms  (1.34)  and  (4.6)  coincide  if  the  coefficients 

satisfy 


c+  =  a+  C|»2»  c+  a+’  c+  R+  "  a+  C4  i ' 


In  particular,  the  first  relation  and  (4.3)  imply  that 

(4.9)  «+<M.A)  *  a+(M,X)  • 

The  normalizing  factor  a+(y,X)  will  be  calculated  in  §5.  The  factors 
R+(P>^)»  T+(U,X)  of  (1.34)  are  independent  of  the  normalization.  Indeed, 
on  combining  (4.8),  (4.9)  and  (4.3)  one  finds 


T  (y  X)  =  - ^iqj. (p,?X) _ 

+^U>  ;  pWt^KtiA) 


(4.10) 
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The  Spectral  Interval  Ap.  For  X  £  A„ ,  Theorem  2.1  and  the 
relations  (4.1)  imply  that  the  only  bounded  solutions  of  <j>  =  X<f>  are 
multiples  of  4>,(y>M«X).  It  will  be  shown  that 


(4.18) 


iMy.U.X)  =  a.  (y  ,A)  <J>3  (y,y, A) 


has  the  asymptotic  form  (1.33).  Indeed,  (4.18),  (4.1)  and  Theorem  2.1 


imply  that 


(4.19)  n(y,y,A)  =  a. 


+  c  ,  e  iq+y  +  o(l),  y  ■+  +oo. 


eqIy[l  +  o(l)] 


The  equivalence  of  (4.19)  and  (1.33)  follows  from 


Lemma  4.3.  For  all  y  >  0  and  1  6  A(  one  has 


(4.20) 


(4.21) 


c32(y,A)  -  c31(u,X)  1  o. 


Comparison  of  (1.33)  and  (4.19)  gives 


C0  *  a0  C32’  C0T0  "  a0>  C 0  R0  =  3 0  C3 1 


Solving  for  cQ,  TQ  and  RQ  and  using  (4.3)  gives 


(4.22) 


(4.23) 


(4.24) 


C  (U  X)  =  a  fu  X)  J  VU 

c0UJ,a;  a0iv,\)  2iq+(vi,A) 


T0(U,X) 


R0 (d  *  A) 


2iq4.(y,A) 

P(°°) [4>34>! 3 (u,X)  ’ 


[4>34>j  1  (y,A) 


The  denominator  [$3$j]  #  0  by  Lemma  4.3  and  relations  (4.3). 
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Finally,  note  that  the  conservation  laws  (1.46)  hold;  i.e., 

(4.25)  q±  1 R± | 2  +  q-  | T± | 2  =  q+  for  all  X  £  A, 

(4.26)  | R0 |  =1  for  all  X  €  A0 . 

In  fact,  relations  (4.25)  are  equivalent  to  relations  (4.7)  and  (4.13) 
of  Lemmas  4.1  and  4.2,  as  may  be  seen  by  combining  (4.7),  (4.8)  and 
(4.13),  (4.14).  Similarly,  relation  (4.26)  follows  froip  Lemma  4.3 
because  (4.21)  implies  that  RQ  =  c31/c32. 

Proof  of  Lemma  4.1.  Relation  (4.7)  can  be  verified  by  calculating 
I  ^1+^4 3  in  two  ways,  using  the  asymptotic  forms  of  <J>  as  y  ■+•  <»  and  y  •*  -<*>. 
Note  that  for  X  £  A  one  has  04  =*  $3  by  the  uniqueness  theorem  (Corollary 
2.5)  and  Theor.m  2.1.  Hence,  calculating  [  4>4 1  at  y  =  -<*>  gives 

(4.27)  [4>4<t>J  *  [<f>„<j>3]  “  2i  p  1  (-“) q_ , 

by  (2.25).  Next,  relations  (4.1)  and  Theorem  2.1  give  (with  the 
notation  c.c.  for  complex  conjugate) 

[<M<J  =  ~  c-c- 

(4.28) 

•  (c41  0!  +  C42  $2){"c41  p_14)J  +  c42  D-1^}  -  C.c. 

=  (c4I  eiq+y  +  c42  e-iq+y) 

x  tc41p  1 (”) (-iq+e~iq+y)  +  c42p-1 (®) (iq+eiq+y) }  -  c.c.  +  o(l) 
-  P  1(°°)(-iq+){  |c41  |2  -  c41c42e2iq+y  +  74lc42e“2iq+y-  |  c4  2  J  2  }  -  c  .c .  +  o(l) 

-  -2iq+  p"1^)  {|c4J2  -  |c4J2}. 
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Combining  (4.27)  and  (4.28)  gives  (4.7). 

Proof  of  Lemma  4.2.  (4.13)  can  be  verified  by  calculating 
[<J>  (J)^]  in  two  ways,  in  analogy  with  the  proof  of  Lemma  4.1.  It  can  also 
be  derived  directly  from  (4.7)  and  the  relations  (4.3). 

Proof  of  Lemma  4.3.  Note  that  for  X  e  A0  the  uniqueness  theorem 
(Corollary  2.5)  and  Theorem  2.1  imply  that  <p3  is  real  valued  and  ^  =  <p2. 
Hence  relations  (4 . 3)  imply 

(4.29)  c31  =  p(»)  [4>34>2]/2iq+  =  p(<*>)  [4>3<J)1]/2iq+  =  c32. 

Moreover,  if  c32  =  0  then  c31  =  0  by  (4.29)  and  hence  by  (4.3)  one  has 

(4.30)  *  [<4>34>2]  -  0. 

But  this  would  imply  that  4> 3  and  <J>2  are  linearly  dependent  which 
contradicts  (2.25).  Hence  cJ2  ^  0. 
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§5.  The  Spectral  Family  of  A  . 

The  eigenfunctions  and  generalized  eigenfunctions  of  §4  are  used 
in  this  section  to  construct  the  spectral  family  (II^(X)}  of  A^.  The 
construction  if  based  on  the  Weyl-Kodaira  theory  as  presented  in  the 
Appendix.  Note  that  the  operator  A^  has  the  form  (A.l)  with  I  ■  R, 
p(y)  *  p(y),  q(y)  =  y2  P-1(y)  and  w(y)  =  c-z(y)  p-1(y).  It  is  clear  that 
p,  q  and  w  satisfy  (A. 2),  (A. 3),  (A. 4)  when  p(y)  and  c(y)  are  Lebesgue 
measurable  and  satisfy  (1.3). 

It  will  be  convenient  to  decompose  R  into  the  disjoint  union 


R  =  Ad  U  {c2(«°)y2}  U  A0  u  (c2  (-°°)vi2 }  u  A 


if  c(°°)  <  c(-°°)  and 


R  =  Ad  u  (c2 (»)u2}  u  A 


if  c(°°)  »  c(-°°)  where 


(5.3) 


Ad  “  Ad(y)  ”  C-00^2  (°°)y2), 


and  A0  and  A  are  defined  by  (4.4).  The  spectral  measures  of  the 
components  of  (5.1)  and  (5.2)  will  be  studied  separately. 

The  Spectral  Family  in  A.  The  spectral  measure  11^  (A)  of  inter¬ 
vals  A  =*  (a,b)  C  A  will  be  calculated  by  applying  the  Weyl-Kodaira 
theorem  to  A  in  A.  The  solution  pair 


(y»A)  -  <j>4(y>y>*) 


<l>,  (y.A)  -  d>,  (y,y,X) 


rmm 
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will  be  used  to  obtain  a  spectral  representation  in  terms  of  the  gener¬ 
alized  eigenfunctions  ip+(y,y,A)  defined  by  (4.5).  The  normalizing 
factors  a+(u,A)  will  be  chosen  after  the  matrix  measure  for  (ipltip2)  has 
been  determined.  Note  that  the  pair  (\p1,\p2)  satisfies  the  hypotheses  of 
the  Weyl-Kodaira  theorem.  Indeed,  (A. 11)  follows  from  Corollary  2.2  and 
(A. 12)  from  Lemma  4.1. 

The  Weyl-Kodaira  theorem  implies  that 

(5.5)  anu(A)ffl 2  =  |  f.(\)  fk(A)  mjk(dA),  A  C  A, 

for  all  f  G  3f(R)  where  (m^k(A))  is  the  spectral  measure  on  A  associated 
with  the  basis  (5.4)  and 


(5.6) 


VX) 


fM  _ 

lim  iMy,A)  f(y)  w(y)  dy, 
M-«°  -M  J 


the  integrals  converging  in  Lz(A,m).  Thus  to  complete  the  determination 
of  II  (A)  for  A  C  A  it  is  only  necessary  to  calculate  (m  (A)}.  Now 

M  J  «C 

IIy(A)  can  be  calculated  from  the  resolvent 
(5.7)  Ry(c)  -  (Ay  -  o'1 

by  means  of  Stone's  theorem  (see,  e.g.,  [16,  p.  79]).  For  A  C  A  the 
theorem  takes  the  form 


iny(A)fU 2 


lim  -^-r  \  (f»[R.  (A+ie)  -  R  (A-ie)]f)dA 

e-*-0+  Z7T1  p  y 


because  iJ0(Ay)  n  A  ■  <j>  by  Lemma  3.1.  Moreover,  Ry(£)  is  an  integral 
operator  in  3C(R)  whose  kernel,  the  Green's  function  of  A,  ,  can  be 


represented  by  the  analytic  continuation  of  the  basis  (5.4)  into  the 
C-plane.  This  procedure,  whose  details  are  presented  in  the  proofs  at 
the  end  of  the  section,  leads  to 


Theorem  5.1.  For  all  f  £  JC(R)  and  all  A  C  A  the  spectral 
measure  II  (A)  satisfies 


(5.9) 


SiyA)f82 


f  {A2(m,A)  |fj (A) I2  +  Af(y,A)  |f2 (A) | 2  >  dA 

*  A 


where 


,  q+(u,A) 

(5.10)  A2 (u» A) - = -  . 

rr  p(i«>)l[<pl<plt]j2 

Corollary  5.2.  The  matrix  measure  (m  (A))  for  the  basis  (5.4) 

3 

is  given  by 


(5.11) 


!  (A) 


■'A 


A2(y,A)dA,  m22 (A) 


f  A2(y,A)dA 
;A 


and  mJ2(A)  *  m21(A)  *  0  for  all  A  C  A. 

These  results  suggest  an  appropriate  choice  of  the  normalization 
factors  a+(li,A)  of  (4.5).  Note  that  if  instead  of  the  basis  (5.4)  one 
takes  (4.5)  then  (5.9)  becomes 


(5.12) 


nny(A)f»2 


where 
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(5.13) 


f+(y,A)  *  11m 


•M  _ 

<My,y»A)  f(y)  »(y)  <iy 

-m  _ 


converge  in  the  space  L2(A,m')  corresponding  to  (<|>+,^_).  This  suggests 
the  choice  |a+|2  =  A2  or 

(5.14)  a±(y,A)  -  ei0± (y »A)  A+(y>X) 

where  9+(y,A)  is  an  arbitrary  real  valued  continuous  function.  The 
matrix  measure  {mj^}  for  (’I'+j’JO  is  independent  of  the  choice  of  the 
phase  factors  exp  {i  9+(y,A)}  and  one  could  take  0+(y,A)  =  0.  However, 
it  will  be  more  convenient  to  choose  9+(y,A)  in  §8  in  a  way  that  simpli¬ 
fies  the  asymptotic  form  of  the  normal  mode  functions  <t>+(y,p,q). 

Theorem  5.1  and  the  above  remarks  imply 

Corollary  5.3.  If  the  basis  (i|>+(y,y, A) ,  ^_(y,y,A))  is 
normalized  by  (5.14)  then  for  all  K  A  one  has 


(5.15) 


#IIM(A)fll2 


•L 


|f+(y,A)| 


I  f  _  (y » A)  |  2 )  dA 


and  the  matrix  measure  (ml^A))  for  0J>+,^_)  is  given  by 

(5.16)  mj t (A)  -  mJ2(A)  -  | A [ 

and  mJ2(A)  =  m21(A)  =  0  where  | A |  is  the  Lebesgue  measure  of  A. 

The  Spectral  Family  in  Ap.  The  spectral  measure  of  intervals 
AC  A(  will  be  calculated  by  applying  the  Weyl-Kodaira  theorem  to  A^ 
Aq  and  the  solution  pair 
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(5.17) 


^i(y.X)  -  4> 3 (y , m, X) , 


^  02 (y » x>  -  (J>x  (y , U, X) . 

The  function  <J>3  is  chosen  to  obtain  a  representation  in  terms  of  the 
generalized  eigenfunction  \p0  defined  by  (4.18).  The  second  function 
could  be  replaced  by  any  independent  solution  of  Ay  $  •  \<p.  The  pair 
(5.17)  satisfies  (A. 11)  by  Corollary  2.2  and  (A. 13)  by  Lemma  4.3. 

Calculation  of  the  spectral  measure  in  Afl  by  the  method  described 
above  leads  to 

Theorem  5.4.  For  all  f  €  3C(R)  and  all  A  C  A0  one  has 


(5.18) 


any(A)fB2 


f  A2 (u,A)  |f x (X) I  2  dA 
* 


where 


q. (u,A) 

(5.19)  A2(y,X)  - - . 

TT  p(»)  |  ]  |  2 

Hence  the  matrix  measure  (m^(A))  associated  with  the  basis  (5.17)  is 
given  by 

(5.20)  mu(A)  -  f  Aj (y ,A)dA 

}  A 

and  m12(A)  =  m21(A)  =  m22(A)  *  0  for  all  A  C  A„ . 

On  replacing  (5.17)  by  the  basis  Oj>0  (y,y,A) ,  <frl(y,y,X))  and 


defining  the  normalizing  factor  by 
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(5.21)  a0(u,A)  -  e1  9°(U,X)  A0(y,X) 

where  0o(y,X)  is  an  arbitrary  real  valued  continuous  function  one 
obtains 

Corollary  5.5.  If  a0(y,X)  is  defined  by  (5.21)  then  for  all 
A  C  one  has 

(5.22)  All  (A)fl2  -  !  |f0  (U,X)  | 2  dX 

^  '  A 


where 


(5.23) 


f0(U,X) 


_ 

lim  ip0  (y,p,X)  f  (y)  w(y)  dy. 
M-*°°  J-M 


In  particular,  the  matrix  measure  (m^(A))  for  the  pair  O|>0,$,,)  is 
given  by 

(5.24)  mJ^A)  *  lAl 

and  mj 2(A)  «  mj t (A)  *  m22(A)  *  0  for  all  A  C  A„ ,  and  the  integral  in 

(5.23)  converges  in  L2(A0). 

The  Spectral  Family  in  A^.  The  portion  of  a(A^)  in  A^  was  shown 
in  §3  to  be  the  set  of  eigenvalues  {X^(y)  |  1  <  k  <  N(y)}.  Moreover, 
each  Xk(y)  is  a  simple  eigenvalue  with  normalized  eigenfunction  \^k(y , yi) 
defined  by  (3.4),  and  corresponding  orthogonal  projection  P  ,  defined  by 

UK 

(5.25)  Pyk  f (y)  =»  (*k(*,y),f)  \(y,u). 

Hence,  recalling  that  by  convention  II^(X)  «  II  (X  .+  0),  one  has 


77 


(5.26) 


nu(X)  f (y)  =  l  OI>k(*,u),f)  ^k(y,u),  X  e  Ad 


\(y)<X 


The  no  Cation  in  (5.26)  denotes  a  summation  over  all  indices  k  such  that 
X,  (y)  <  X.  The  sum  in  (5.26)  is  finite  for  all  X  6  A,. 

K  "  Q 

The  Spectral  Measure  of  the  Points  c2(°°)y2  and  c2(-°°)y2.  When 
c(°°)  <  c(-°°)  one  always  has  IT^({c2  (-°°)y2})  =  0  because  in  this  case  the 
special  solutions  <j>i  (y.y,X) ,  <J>2  (y,y,X)  are  defined  for  X  =  c2(-<»)y2  and 
Theorem  2.1  implies  that  A^  <p  ~  c2  (— °°) yi2c|)  has  no  solutions  in  3f(R).  The 
point  X  =  c2 («) y2  may  be  an  eigenvalue  of  A  .  In  each  case  this 
question  must  be  decided  by  determining  the  behavior  of  solutions  of 
A  <f>  *  c2(<»)y2<f>  for  y  -►  ±°°.  Theorem  2.6  gives  simple  sufficient  condi¬ 
tions  for  II^({c2  (°°)y2})  3  0.  For  simplicity,  it  will  be  assumed  in  the 
remainder  of  the  report  that  c2(®)y2  £  afl  (A^)  .  In  cases  where  c2(°°)y2 
is  an  eigenvalue  a  corresponding  term  must  be  added  to  the  eigenfunction 
expansion. 

The  Eigenfunction  Expansion  for  A  .  Combining  the  representa¬ 


tions  of  II  obtained  above,  one  finds  the  representation 


N(y)-1 

iny(X)fll2  =  l  H(X  -  Xk(y))  jfk(y) 


k=l 


(5.27) 


H(X  -  X’)  | f  0  (y , X ’ ) | 2  dX' 


H(X  -  X ' ) ( | f+(y ,X ' ) | 2  +  | f_(y ,X  ') | 2)  dX ' 


where  H(X)  ■  1  for  X  >  0,  H(X)  =  0  for  X  <  0, 
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(5.28)  fk(p) 


JR 


^(y.y)  f(y)  C~2(y)  p-1  (y)  dy,  1  <  k  <  N(y), 


(5.29)  f0(y,A)  =  L, (A 


2'“0 
M-*» 


fM  _ 

)-lim  ^0(y,y,A)  f (y)  c~2(y)  p-1(y)  dy, 
-k»  J-m 


and 


(5.30)  f+(y,A)  =  L, (A)-lim 


M-wo  " 


M 


f(y)  c“z(y)  p-1  (y)  dy. 


In  particular,  on  making  A  -+■  <*•  one  obtains  the  Parseval  relation  for  A^: 

(|f+(y,A)|2  +  |  f_(y  -A)  | 2 )  dA  . 


N(y)-1  A 

(5.31)  Ilf#  2  =■  f  |  f  k  (y)  |  2  + 

k=l 


I  f0  (y»A)  I  2dA  + 


Ja 


Thus  the  correspondence 

(5-32>  f  ^  4'yf  =  (f+(y.*),f_(y,*),f0(y,‘),f1(y),f2(y),***) 

defines  an  isometric  mapping  4^  of  3C(R)  into  the  direct  sum  space 

(5.33)  L2(A)  +  L2 (A)  +  L2(A0)  +  CN(m)_1 
and  one  has 

Theorem  5.6.  S'  is  a  unitary  operator  from  Jf(R)  to  the  space 

(5.33) . 

This  result  will  be  shown  to  follow  from  the  Weyl-Kodaira  theorem 
and  the  corresponding  properties  of  the  partially  isometric  operators 
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Corollary  5.8.  The  inverse  isometries  ¥  ,  0  <  k  <  N(y), 

are  given  by 


('*',. +f+)  (y)  =  X(R) 
U_  •  <5-0 


(5.39) 


fM 

)-lim 

,M-°°  ‘  c2  (- 


(^foHy)  -*«o 
6-0 


°°)y2+6 

c2(-°°)y2-<$ 


iMy.y.A)  f+(X)dX, 


/•C 

-lin  ^0(y,y,A)  f0(X)dX, 

0  •'c2(00)y2+6 


(^yl.fk)(y)  =  fk  it»k(y,y),  l  5  k  <  N(y). 


The  spectral  property  of  the  unitary  operator  V  is  described  by 
Corollary  5.9.  For  every  f  €  D(A^)  one  has 

(5.40)  ¥u  Auf  =  (Xf+(p,X),Af_(y,X),Xf0(p,X),X1(y)f1(y),“*). 

This  completes  the  formulation  of  the  results  of  §5  and  the 
proofs  will  now  be  given. 

Proof  of  Theorem  5.1.  The  integral  representation  of  the 
resolvent  will  be  used.  Thus,  as  in  the  proof  of  Theorem  3.2, 


(5.41) 


Ry(0  f(y)  « 


Gy(y,y',0  f(y’)  w(y')  dy\  ;  £  a(Ay), 


_  2  _  2 

where  w(y)  =  c  (y)  p  (y) ,  and  the  Green's  function  has  the  form 
(5.42)  Gu(y,y',0  =  [<)>„(  •  ,Z)<P_„(  *  ,0  ]"14>_00(y<^)  <l>00(y>,0 

where  y<  =  y<(y,y')  =  Hin  (y,y'),  y>  =  y>(y,y')  =  Max  (y,y’)  and 
4>+00(y,O  are  non-trivial  solutions  of  0  =  £<j>  such  that  <{>_„(  *,0 
£  L2(-°°,0)  and  (J00(*,?)  e  L2(0,°°).  To  identify  these  solutions  note 
that  by  (2.4) 


T~  T 


Im  q+(y,£)  >  0  for  £  6  R+(c(-°°)y)  int 


(5.43) 


and 


(5.44) 


Im  q+(jj,5)  <  0  for  £  £  R_(c(-°°)y) 


int 


eJyqi  .  eii  y  q±  =  e;y  im  q+±i y  Re  q+  _ 


It  follows  from  Theorem  2.1  that 


(5.45) 


and 


(5.46) 


^(y,  O  =  ^(y.u,  0 
=  <t\(y»u»o 


)•  for  £  £  R+(c(-°°)y)inC 


4>00(y>C)  =  4>2(y»u,C) 


^Jy.C)  =  <t>3(y»u,Q) 


for  Q  6  R  (c(-°°)y) 


int 


Now  the  functions  appearing  in  (5.45)  and  (5.46)  have  continuous  exten¬ 
sions  to  R+(c(-°°)u)  and  R  (c(-°°)y),  respectively,  by  Corollary  2.2. 
Indeed,  (5.42),  (5.45),  (5.46)  and  Corollary  2.2  imply 
Lemma  5.10.  For  all  X  £  A  one  has 

(5.47)  G)J(y,y'  ,A+iO)  =  (y^U.  A)  <j>j  (y>,u»X) , 


(5.48) 


G  (y,y',A-iO)  =  [<J>2<f>3]QN  4>3(y<,U,A)  4>2(y>,u,A), 


u 


(A) 


and  the  limits  are  uniform  on  compact  subsets  of  R  x  R  x  A. 

To  prove  Theorem  5.1  it  is  clearly  sufficient  to  verify  (5.9) 


for  the  functions  f  of  a  dense  subset  of  )C(R) .  It  will  be  convenient 


to  use  the  subset 


(5.49)  5<com(R)  =  3C(R)  G,  {f  )  supp  f  is  compact}. 

Note  that  for  f  S  3CCom(R)  the  integrand  in  Stone's  formula  (5.8)  can 
be  written 

(f,[Ry(X+ie)  -  Ry(A-ie)]f) 

(5.50) 

{Gy(y,y',A+ie)  -  Gy(y,y ’ ,A-ie) }  f(y)  f(y')  w(y)  w(y')  dydy' 

Lemma  5.10  implies  that  for  all  f  £  3fcoin(R)  this  expression  tends  to 
a  limit 


1 

2iri 


(5.51) 


(f,Hy(A)f) 


Hy(y,y',A)  f (y)  f(y')  w(y)  w(y')  dydy', 


uniformly  on  compact  subsets  A  C  A,  where 


(5.52) 


Hy(y,y',A) 


2^  {Gy(y,y',A+iO)  -  Gy  (y,y' ,A-iO)}. 


It  follows  that 


(5.53) 


n  ny  (A)  fii 2 


(f,H  (A)f)  dA 


for  all  f  S  •Kcom(R).  The  proof  of  Theorem  5.1  will  be  completed  by 

'it  — 

calculating  Hy(y,y',A).  Note  that  the  well-known  property  Ry(£)  *  Ry(;;) 

implies  that  Gy(y',y,£)  -  Gy(y,y',£)  and  hence 


Similarly 


(5.61)  <}>2  (y '  ,y ,  A)  =  ^(y’.y.A) 
and 

(5.62)  <l>3 (y»U, A)  =  c31(y,X)  ^(y.y.A)  +  c34(y,X)  4>4(y»U»X) 
where 

(5.63)  c31  -  [<D30>4]/[<j>1<DJ,  c34  =  [<D3<D1]/[(J)4<J.1]. 

Combining  these  relations  and  Lemma  5.10,  one  finds  for  X  6  A,  y  <  y' 

Gu(y,y' ,x+i0)  =  [<J)14>4]-1  <t>4(y,M,A)  (y  ’ , y , X) 

(5.64) 

*  [ <t>! ]  1{c2J  4>4(y,y,A)  4> t (y ' ,u.  A)  +  c^  <J>4(y,y,X)  <f4(y’,u,X)} 

and 

Gy (y»y ' » X-iO)  =  [<()2<t)3]_1  <f>3 (y»M. X)  <t»2(y',h,X) 

(5.65) 

=  [^2^3 1_1  3 x  <t>i(y»y»A)  4>x(y',y,X)  +  c34  <(>4(y,y,A)  ^(y'.y.A)}. 

Combining  (5.52),  (5.64)  and  (5.65)  gives 


Hy(y,y',X) 

(5.66) 


2-rri 


My.u.A)  <j)4(y',y»X)  -  4»1(y,y,x)  <j>1(y,,y,X) 


£iO_ 


2T3  • 


+  (ifil  '  HAll  'My,',,X) 


To  calculate  the  coefficients  in  (5.66)  recall  that  <|>  (y,y,A)  *  <j>2  (y,y,A) , 
4>j(y,U,X)  -  <J>i»(y»Ut^)  for  X  €  A.  It  follows  from  (5.60),  (5.63)  and 


(2.25)  Chat 


(5.67) 


(5.68) 


— E<t>  1  <J>2  ' 

f  ..  2ip_1  (»: 

umji 

2 

-  [03<J)4 

1  2ip_1  (-< 

1 C4^!  $4  ] 

2  ‘  |[<!>A 

=  2iTi  A2 , 


(5.69) 


C2  1  C34 


..  IMJ  [4>24>3]  (Mi][<M3] 

Thus  (5.66)  can  be  written 

Hu(y,y',X)  =  a_£  (J)4(y,p,X)  <j>4(y '  ,y, X)  +  a2  (y, m, X)  4>j (y ’ »p. X) 


(5.70) 


A+  'l'i(y.X)  ip  1  (y '  ,X)  +  Af  ip2(y,X)  ip2(y',\). 


This  was  proved  for  y  <  y'.  However,  both  sides  of  (5.70)  are  Hermitian 
symmetric  and  it  therefore  holds  for  all  (y,y')  e  R2.  Multiplying 
(5.70)  by  f (y)  f(y')  w(y)  w(y')  and  integrating  over  R2  gives 


(5.71) 


(f,Hy(X)f)  =  A_£  j f !  ( X )  I  2  +  A 2  |f2(X) 


and  combining  (5.53)  and  (5.71)  gives  (5.9).  This  completes  the  proof 
of  Theorem  5.1. 

Proof  of  Corollary  5.2.  It  follows  from  (5.5)  and  (5.9)  that 


(f,nu(A)g)  =  j  f j (X)  gj (X)  mjk(dX) 


(5.72) 


{A*  ft(X)  g, (X)  +  A2  f_(X)  g2(X)}  dX 

A 


for  all  f,g£  3f(R) .  Moreover,  the  Weyl-Kodaira  theorem  implies  that 


w 
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?  ■  (f i » f 2 )  and  8  =  (g1>i2)  can  be  arbitrary  vectors  in  L2(A,m).  Now 

A 

all  vectors  f  with  bounded  Borel  functions  as  components  are  in  L2 (A,m) . 
Thus  the  first  of  equations  (5.11)  can  be  obtained  from  (5.72)  by  taking 

A  A  A  ^ 

fj  =  gj  =  X^»  the  characteristic  function  of  A,  and  f2  =  g2  =  0.  The 
remaining  equations  of  Corollary  5.2  are  obtained  similarly. 

Proof  of  Theorem  5.4.  The  proof  follows  that  of  Theorem  5.1. 
Equations  (5.41),  (5.42)  for  the  resolvent  are  still  valid.  However, 
instead  of  (5.43)  one  has 


(5.73) 


Im  q+(p.O  >  0  for  ?  6  R+(c(»)y) int, 
Im  q+(p,;)  <  0  for  ?  £  R~(c (»)y)int. 


Im  q_(p,5)  <  0  for  Z  e  R(c(°°)y)  n  L(c(-°°)y). 


Thus  by  Theorem  2.1 

^(y.S)  =  <!>1(y»p»5) 

(5.74)  .  for  Z  €  R+(c(°°)y)int  n  L(c(-«)y), 

“  4>3 Cy.y 

and 

^(y.c)  =  $2(y,y,c) 

(5.75)  •  for  Z  £  R“(c(»)y)int  n  L(c(-°°)y). 

^(y.s)  “  M y,u.O 

Hence  Corollary  2.2  implies 

Lemma  5.11.  For  all  X  S  A0  and  y  <  y'  one  has 

(5.76)  G^(y,y'  ,X+iO)  -  C  ]  1  <J>3(y,y,X)  (^(y'.U.A), 

(5.77)  Gy(y,y'  ,X-iO)  -  <J>3(y,y,A)  <J>2  (y*  ,p,X)  , 
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and  the  limits  are  uniform  on  compact  subsets  of  R  x  R  x  . 

Proceeding  to  the  calculation  of  H^(y,y',X)  one  has 

(5.78)  ■  for  z  e  R(c(“)y)  n  L(c(-“>)u) 

qUy.T)  =  ql(y,?) 

whence 

(5.79)  4>i(y.y.O  **  <l>2(y»u,D  for  R+(c(“)u)  n  L(c(-»)p) 
and 

(5.80)  4>3 (y.P.O  =  4>3(y,y.?)  for  ;  e  R(c(“)y)  n  L(c(-”)u). 
It  follows  that 


(5.82)  c21  =■  “d  c23  =  [4>2<l>1  ]/ [4>34>1  ]  • 

These  relations  and  Lemma  5.11  imply  that  for  X  6  A0  and  y  <  y', 

G/y.y'.x+io) 

(5.83) 

*  [‘t>i<J>3r1{c2i  4>3(y,y,X)  ^(y'.U.X)  +  c23<{>3(y,u,X)  <j>3 (y ' ,y,X) } 

end 

Gy(y,y' ,X-iO)  -  td>2d>3 ]-1  $3(y.u,X)  <j>j (y ’ ,y,A) . 


(5.84) 


Thus  (5.85)  can  be  written 


(5.88) 

for  y  <  y 
for  all  f 

(5.89) 


Hu(y,y\A)  =  a£  4>3 (y,y, A)  <J>3(y T ,u,X)  =*  A2  iMy>*)  (y ' ,X) 

'  and  hence  for  all  (y,y’)  €  R2.  It  follows  by  integration  that 
€  Hcom(R)  one  has 


(f,Hy(A)f)  -  A|  ! f0 (X) 1 2 ,  ASA,. 


Combining  (5.89)  and  (5.53)  gives  (5.18).  Finally,  (5.18)  implies 
(5.20)  by  the  argument  used  to  prove  Corollary  5.2.  This  completes  the 
proof  of  Theorem  5.4. 


Proof  of  Theorem  5.6.  4*  is  isometric  by  (5.31),  (5.32).  Hence 
to  prove  that  4*  is  unitary  it  is  only  necessary  to  prove  that  it  is 
surjective.  But  this  is  an  immediate  consequence  of  Theorem  5.7, 
equations  (5.38).  The  latter  are  implied  by  the  Weyl-Kodaira  theorem. 

Proof  of  Theorem  5.7  and  Corollaries  5.8  and  5.9.  These  results 
are  direct  consequences  of  the  Weyl-Kodaira  theorem. 
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§6 .  The  Dispersion  Relations  for  the  Guided  Modes. 

The  eigenvalues  of  A^  determine  the  relation  between  the  wave 
number  |p|  and  frequency  u)  =  ,  or  dispersion  relation,  for  the 

guided  mode  functions  ^(x.y.p).  The  functions  Ak(p)  also  appear  in  the 
definition  (3.4)  of  ^k(y,p) .  The  purpose  of  this  section  is  to  provide 
the  information  concerning  the  p-dependence  of  Ak(p)  and  i|>k(y,p)  that  is 
needed  for  the  spectral  analysis  of  A  in  §§7-8. 

The  domain  of  definition  of  the  function  A^(p)  is  the  set 

(6.1)  0k  =  {p  |  N(p)  >  k  +  l},  k  =  1, 2, • • • . 

Note  that  0^  is  not  empty  if  and  only  if  1  <  k  <  N0  where 

(6.2)  N0  *  sup  N(p)  <  +°°. 

p>0 

Clearly,  if  N0  <  +°°  then  N  -  1  is  the  maximum  number  of  eigenvalues  of 
A^  for  p  >  0.  If  N0  *  +o°  then  either  a0(Ay)  is  infinite  for  some  p  >  0 
or  a0(A^)  is  finite  for  all  p  >  0  and  N(p)  -*■  00  when  p  ■+  °».  Theorem  3.14 
implies  that  both  cases  occur.  The  principal  result  of  this  section  is 
Theorem  6.1.  For  1  <  k  <  NQ  the  set  0 k  is  open  and  Ak  :  -*■  R 

is  an  analytic  function. 

The  proof  of  this  result  given  below  is  based  on  analytic  per¬ 
turbation  theory  as  developed  in  [8]. 

The  curves  A  *  A^(p)  »  U  £  0^,  can  never  meet  or  cross  because 
each  eigenvalue  is  simple  and  the  corresponding  eigenfunction  if>k(y,p) 
has  exactly  k  zeros  (Theorem  3.8).  Thus  for  1  <  k  <  N0  -  1  one  has 
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(6.3)  c2p2  <  Ak(y)  <  Ak+i(y)  <  c2(°°)p2,  p  6  Ofc+i . 
Moreover,  if  (?k  is  unbounded  then  (6.3)  implies 

(6.4)  c2  <  lim  inf  p  2  A,  (p)  <  lim  sup  p  2  A.(p)  <  c2(°°). 

y-*»  |j-*x> 

In  particular,  if  0k  is  unbounded  then 

(6.5)  lim  A,  (p)  =»  +“. 

p-«o 

A  related  property  is  given  by 


Theorem  6.2.  For  1  <  k  <  N0  the  function  Ak(p)  is  strictly 
monotone  increasing;  i.e.,  for  all  P1#  p2  €  0  one  has 

(6.6)  Ak(Pi)  <  Ak(p2)  when  Pi  <  p2. 

The  proof  of  (6.6)  given  below  is  based  on  a  variational 
characterization  of  Ak(p) . 

By  Theorem  6.1,  ^k  is  open  and  is  therefore  a  union  of  disjoint 
open  intervals.  Hence  the  curve  A  =  Ak(p)  consists  of  one  or  more 
disjoint  analytic  arcs.  It  is  interesting  that  these  arcs  can  terminate 
only  on  the  curve  A  =  c2(°°)p2.  More  precisely,  one  has 

Corollary  6.3.  Let  p0  be  a  boundary  point  of  $k.  Then 

(6.7)  lim  Ak(p)  =  c2(»)p2. 

U-’-Po 

It  is  clear  from  Theorem  6.2  and  (6.3)  that  the  limit  in  (6.7) 
exists  and  does  not  exceed  c2(°°)p2.  The  equality  (6.7)  is  proved  below. 

The  result  (6.4)  can  be  improved  by  strengthening  the  hypotheses 
concerning  c(y) .  A  result  of  this  type  is 


Hiti  |  i.MP— 


Theorem  6.4.  Let  cffi  <  c(“)  and  assume  that  for  each  e  >  0  there 
is  an  interval  1(e)  C  R  such  that 

(6.8)  c(y)  <  cm  +  e  for  all  y  6  1(e). 

Then  N0  =  +®,  0^  is  unbounded  for  each  k  >  1  and  XR(y)  ~  c^y2  when 
M  -*■  00  in  0^;  i.e. , 

(6.9)  lim  y-2  X,  (y)  =  c2  . 

V^oo  k  m 

The  analyticity  of  XR(y)  and  Corollary  2.2  imply  the  continuity 
of  the  eigenfunctions  ^(y.y).  More  precisely,  one  has 

Corollary  6.5.  For  1  <  k  <  N0  the  function  ^(y.y)  satisfies 

(6.10)  V  P_1  6  C(R  x  0k). 

This  completes  the  formulation  of  the  results  of  §6. 

Proof  of  Theorem  6.1.  The  analytic  perturbation  theory  of 
[8,  Ch.  VII]  will  be  used.  Note  that  the  operator  may  be  defined  for 
all  Vi  £  C  by  (1.16),  (1.18)  and  is  a  closed  operator  in  5C(R) .  Moreover, 

the  domain  D(A^)  is  independent  of  y  and  is  a  Hilbert  space  with  respect 

to  the  norm  defined  by 

(6.11)  'Kcnj  =  ii^(R)  +  «<n2<R)  +  kp“V>’i£(r). 

It  follows  that  y  -*■  is  holomorphic  in  the  generalized  sense.  Indeed, 
in  the  definition  of  [8,  p.  366]  one  may  take  Z  =  D(A)J)  (independent  of 
y)  and  define  U(y)  :  Z  -*■  3C(R)  to  be  the  identification  map.  Then  U(y) 
is  bounded  holomorphic  (in  fact,  constant)  and 

V (y ) 4>  -  Ay  U(y)<j>  -  -c2p  [(p-1<J>')'  -  p-1y2<]>] 


(6.12) 


is  holomorphic  for  all  tp  e  Z.  Thus  is  holomorphic.  It  follows  by 
[8,  p.  370]  that  each  Ak(p)  has  a  Puiseux  expansion  at  each  point 
P0  6  0^.  But  each  X^(u0)  is  simple  (Lemma  3.3)  and  hence  the  Puiseux 
series  can  contain  no  fractional  powers  of  p  -  pQ  [8,  p.  71].  Thus 

is  in  fact  analytic  at  each  pfl  €  0 This  proves  both  statements  of 
Theorem  6.1. 

Proof  of  Theorem  6.2  and  Corollary  6.3.  The  eigenvalue  Ak(p) 
can  be  characterized  by  the  variational  principle  [3,  pp.  1543-4] 

(6.13)  A,  (u)  =  inf 

where  denotes  the  set  of  all  k-dimensional  subspaces  of  3f(R). 
Moreover,  D(A^)  is  independent  of  p  and 

(6.14)  (A $,$)  *  <  i  4>*  <y)  1 2  +  u2  |$(y)|2)  p_1(y)  dy 

M  R 

for  all  $  S  D(A^)  [8,  p.  322].  Hence  if  <  p2  then 

(6.15)  (A  <*> , 4>)  <  (A  $,$) 

Pj  p2 

for  all  4>  £  D(A  )  =  D(A  ).  In  particular,  (6.13)  and  (6.15)  imply 
Ul  P2 

that  if  px ,  p2  S  0k  then 

(6.16)  Ak(p,)  <  Ar(p2) 

which  proves  the  weak  monotonicity  of  Ak>  It  will  be  convenient  to  use 

(6.16)  to  prove  Corollary  6.3  before  proving  the  strong  raonotonicity . 

To  prove  Corollary  6.3  note  that  (6.16)  and  (6.3)  imply  that  the 
limit  in  (6.7)  exists  and  does  not  exceed  c2(°°)p0.  But  if 


sup 

$eAH)(A  ) 
II  $11  =1  M 


(Ay<p,q>7 


|> 
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lim  Ak(y)  =  A?  <  c2(°°)y2  then  Corollary  2.2  implies  that  y0 ,  X?  satisfy 

(6.17)  F(y0,X“)  =  0; 

see  (3.28).  Moreover,  F(y,A)  is  analytic  at  Uo>  A£  hy  Corollary  2.3. 

It  follows  that  A®  is  an  eigenvalue  of  and  hence  y0  €  0^  by  Theorem 
6.1.  This  contradicts  the  assumption  that  U0  is  a  boundary  point  of  0^. 

To  prove  that  each  A^(p)  is  strictly  mono tonic  in  0^  two  cases 
will  be  considered.  First,  if  At,  X2  are  in  the  same  component  of  0^, 
say  (a,b)  C  0^,  then  X^(Ui)  =  X^(y2)  would  imply  that  A^(p)  =  const,  in 
[M1,U2]  and  hence  in  (a,b),  since  Ak(y)  by  Theorem  6.1.  But  this 
contradicts  Corollary  6.3  since 

(6.18)  c2(°°)a2  =  lim  A,  (y)  <  lim  A,  (y)  =  c2(°°)b2 

V*a  y+b 

In  the  second  case  yx  and  y2  lie  in  different  components  of  0^,  say 
Mi  6  (aj.bj)  C  0^  and  y2  £  (a2»b2)  C  0^  with  bt  <  a2 .  In  this  case,  by 
the  preceding  argument  one  has 

(6.19)  Ak(p1)  <  c2(»)b2  <  c2  (°°) a|  <  Xk(y2) 
which  completes  the  proof. 

Proof  of  Theorem  6.4  (sketch).  The  proof  is  based  on  the  method 
proposed  for  the  proof  of  Theorem  3.16  and  the  variational  principle 
(6.13).  Note  that  the  hypothesis  (6.8)  and  Theorem  3.16,  generalized  to 
non-constant  p(y),  imply  that  N(y)  -*•  °°  for  y  -*■  <*>.  Hence  N0  =  +°°  and  each 
0k  is  unbounded. 

To  prove  (6.9)  choose  piece-wise  constant  functions  c0(y)  and 
c°(y)  such  that 
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(6.20) 


c0(y)  <  c (y)  <  c°(y)  for  all  y  e  R, 


(6.21) 

(6.22) 


c0(y)  =  cm  <  c0(°°)  on  an  interval  Ifl,  and 


+  £  =  c(y)  <  c  (“)  on  an  interval  1(e) 


The  notation 


(6.23) 


3fn(R)  =  L2(R,c~2(y)  p  1  (y)  dy) 


jf°(R)  =  L2(R,c°  2  (y)  p  1  (y)  dy) 


will  be  used.  The  three  spaces  7f(R) ,  JCQ  (R)  and  JC°  (R)  have  equivalent 
norms.  In  particular,  if  II  $11 0  and  II  <J>1  0  denote  the  norms  in  ~K0  (R)  and 
5C° (R) ,  respectively,  then  by  (6.20) 


(6.24) 


U  cf>M  u  <  H  <  I  <J)II  Q . 


Now  note  that  the  variational  principle  can  be  formulated  in  the 
homogeneous  form 


(6.25) 


where 


X,  (p)  =  inf  sup 


Mes,  <t>eVinD(A  )  upu2 

*  &o  w 


(6.26) 


(Ay4>* 4>)  /R(|<t>'  j  2+U2  !$|  Z)P~  (y) dy 

U  4>ll 2  II  pll 2 


Both  the  numerator  in  (6.26)  and  D(A^)  are  independent  of  c(y).  Hence, 
if  A0^andA®  denote  the  operators  corresponding  to  p(y) ,  co(y)  and  z(y) 
c°(y)  respectively,  then  (6.24)  and  (6.26)  imply 
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(6.27) 


(Aqu<M)  (Ay<M) 

HI  l  l$l2  1 4>l 


for  all  <j>  S  D(AQy)  *  D(A^)  ■  D(A^)  such  that  <j)  ^  0. 
principle  (6.25)  and  (6.27)  imply  that 


(6.28)  Xok(y)  <  Xk(u>  <  X®(y) 

for  all  sufficiently  large  y,  where  f^ok(^)}  and  (X 
eigenvalues  of  Afly  and  A^,  respectively. 

The  proof  of  (6.9)  can  now  be  completed  by 
calculation  that  (see  [17]) 


(6.29) 


lio  y~2  X  .  (y)  ■  c* 

0lC  m 

lim  y"2  X£(y)  »  (cm  +  e) 


It  follows  from  (6.28),  (6.29)  that 


(6.30)  c*  <  lim  inf  y”2  X.  (y)  <  lim  sup  y  2  X. 

m  "  y-K»  k  ‘  y-«o 

Equation  (6.9)  follows  because,  by  hypothesis,  e  > 
Proof  of  Corollary  6.5.  The  result  (6.10) 
Corollary  2.2,  the  relation 


(6.31)  4>k(y,y)  -  My»y.*k<w>>/,,M*»' 


»♦) 

02 

The  variational 


k(y) }  are  the 
showing  by  direct 


2 

(y)  <  (cm  +  e)2. 

0  is  arbitrary, 
is  immediate  from 

,Xk(y))l 


and  Theorem  6.1. 


'mmmim  fa os  bunk-jot  films® 


§7.  The  Spectral  Family  of  A. 

The  acoustic  propagator  A  was  defined  in  §1  and  shown  to  be  a 
selfadjoint  non-negative  operator  in  the  Hilbert  space  3C.  In  this 
section  the  spectral  family  {II(\)  |  A  >  0}  of  A  is  constructed  by  means 
of  the  normal  mode  functions  {'P+*'P_»'Po»'P1»  *** K  The  method  of  construc¬ 
tion  is  to  use  Fourier  analysis  in  the  variables  x  €  R2  to  reduce  A  to 
the  operator  A|p|  and  then  to  use  the  spectral  representation  of  lI|p|(X) 
developed  in  §5.  The  construction  is  given  in  Theorems  7. 1-7. 4.  In  the 
remainder  of  the  section  the  proofs  of  Theorems  7. 1-7. 4  are  developed  in 
a  series  of  lemmas  and  auxiliary  theorems. 

The  formal  definitions  of  the  normal  mode  functions  i|»+(x,y,p,A) , 
ipQ(x,y, p,A)  and  ip^(x,y,p)  were  given  in  SI,  equations  (1.29)- (1.31)  and 
(1.36)-(1.41) .  The  definitions  were  completed  by  the  construction  of 
the  special  solutions  4>j(y,W,X)  in  §2  and  of  the  normalizing  factors 
a+(y,A),  a0(y,A),  a^y)  in  §5.  The  construction  of  11(A)  will  be  based 
on  these  normal  mode  functions  and  Che  corresponding  generalized  Fourier 
transforms.  Formally  the  latter  are  the  scalar  products  of  functions 
f  €  X  with  the  normal  mode  functions.  The  following  notation  will 
be  used. 

(7.1) 

(7.2) 

(7.3) 
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f+(p,A)  -  [  ip  (x,y,p,A)  f (x,y)  c~2(y)  p_1(y)  dxdy, 
Jr3 


f0(p,A)  -  [  'Pj (x,y,p,A)  f (x,y)  c“2(y)  p_1(y)  dxdy, 
JR3 


fk(p)  *  ^(x.y.p)  f (x,y)  c  z(y)  p  (y)  dxdy,  k  >  1. 


LOO 


Of  course  these  Integrals  need  not  converge  since  the  normal  mode  func¬ 
tions  are  not  in  X.  Instead,  they  will  be  interpreted  as  Hilbert  space 
limits  as  in  the  Plancherel  theory  of  the  Fourier  transform.  This 
interpretation  will  be  based  on  the  following  three  theorems. 

Theorem  7.1.  If  f  6  L^R3)  then  the  Integrals  in  (7.1),  (7.2), 
(7.3)  are  absolutely  convergent  for  (p,A)  €  ft,  (p,A)  €  ft#  and  p  €  (y 


respectively,  and 


(7.4) 


f±  6  c(ft) ,  f0  e  c(ft0),  fk6  coy,  k  >  l. 


For  each  f  €  X  and  M  >  0  define 


(7.5) 


fM(x,y) 


f(x,y)  if  |x|  <  M  and  |y|  <  M, 


0  if  |x|  >  M  or  |y|  >  M. 


It  is  clear  that  f M  -►  f  in  X  when  M  Moreover,  fjj  €  X  n  L,(R3)  and 


one  has 


Theorem  7.2.  For  every  f  €  Jf  and  M  >  0, 


(7.6) 


€  L2  (0),  €  L2  (n„),  fjjk  e  Lj(y),  k  >  1, 


and  the  Parseval  relation  holds: 


(7.7) 


^m'x"  1  fM+l  i2  («)  +  '^M-'LjW)  +  k^0  'WLj^) 


The  relation  (7.7)  suggests  the  introduction  of  the  direct  sum 


space 


(7.8) 


Nb-1 

X  -  Lj (ft)  +  Li (ft)  +  2  L 2 ' 

k-0 
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X  is  a  Hilbert  space  with  norm  defined  by 


(7.9) 


|h'£ 


lh+'L 


(0) 


lh  iT2 


(£)) 


No-1 

I 

k-0 


lh, 


'l. 


«k>’ 


see  [3,  p.  1783].  Theorem  7.2  Implies  that  for  each  f  €  X  and  M  >  0, 
the  sequence  ?M  -  (^‘^M-’ W*Mi’ 6  *  and 

(7.10)  »Vx  '  '¥)(• 


For  arbitrary  f  €  X  the  generalized  Fourier  transforms  associated  with  A 
are  defined  by 

Theorem  7.3.  For  all  f  €  X,  {fM>  is  a  Cauchy  sequence  in  X,  for 
M  -*■  and  hence 

(7.11)  lim  t,  -  f  -  (f.,f  ,?,?,•••) 

M^-M  +  -  o  i 

exists  in  X.  In  particular,  each  of  the  limits 


(7.12) 


f+  -  L2(R)-lim  f„+ 
f0  «  L2(nB)-lim 

M-*» 

fk  -  L2(flk)-lim  k  >  1, 


exists  and  the  Parseval  relation 


(7.13) 


lfl£  -  Ifl^  -  |f+I 


L2(«) 


+  If  I 


-  L.CD 


No-1 

L  l^v^r  (ft,  } 

k*0  K 


holds  for  every  f  e  X. 
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Theorem  7.3  associates  with  each  f  6  X  a  family  of  generalized 
Fourier  transforms  f  ■  (f+,f_,f # ,f t , • • •)  6  X  such  that 

(7.14)  f+(p,A)  -  L2 (£2) -11m  [  [  'J>+(x.y.p.A)  f (x,y)c-2 (y)p-I(y)dxdy, 

M4«  * -M  *  1 x | <K  " 


(7.15)  f0(p,A)  -  L2(fl0)-lim  [  f  i{>0 (x.y.p.A)  f (x,y)c~2 (y)p  ‘(y)dxdy, 

M-*»  ‘  -M  ;  Ixl^l 

(•M  t  _ 

(7.16)  f.(p)  *  L2(%)-llm  ip.  (x.y.p)  f  (x,y)c-2  (y)p“l  (y)dxdy, 

M-*»  J -M  *  |x| <M 

k  >  1. 

It  is  easy  to  verify  that  if  f  €  X  n  Lj (R3 )  then  the  functions 
(f+,f_,f0,f1,***)  defined  by  Theorems  7.1  and  7.3  are  equivalent  and  hence 
the  notation  is  unambiguous.  A  construction  of  the  spectral  family 
Ol(y)}  based  on  these  functions  is  described  by 

Theorem  7.4.  For  all  f,g  €  3C  and  all  real  y  >  0,  H(y)  satisfies 
the  relation 


(f.IIOOg)  -  [  H(y-A)  (f ,  (p,A)  g.(p»A)  +  l_(p,A)  g_(p,A)  dpdA 


(7.17) 


+  f  H(y-A)  f0(p,A)  g.  (p,A)  dpdA. 

Na-l  r  - - 

+  I  H(y-A.(|pj))  f.(p)  g.(p)  dp 

k-1  \ 


where  H(y)  »  1  for  y  >  0  and  H(y)  »  0  for  y  <  0. 

The  remainder  of  §7  presents  the  proofs  of  these  theorems.  The 
proof  of  Theorem  7.1  will  be  based  on 
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(7.18) 


Lemma  7.5.  The  normal  mode  functions  satisfy 
<J>±(x,y,p,A)  6  C(R3  x  Q), 


ipQ  (x»y»p»A)  €  C(R3  x  Qq), 

^k(x,y,p)  e  C(R3  x  l^),  k  >  1. 

Moreover,  for  each  compact  set  K  C  Q  there  exists  a  constant  such 
that 


(7.19)  |<M*,y.P,*)|  <  \  for  all  (x.y)  6  R3  and  (p,X)  €  K. 

Similarly,  for  each  compact  K  C  there  exists  a  constant  such  that 

(7.20)  |i|/0 (x,y,p,A)  |  <  MR  for  all  (x.y)  S  R3  and  (p,A)  6  K, 

and  for  each  k  >  1  and  compact  K  C  £1  there  exists  a  constant  such 
that 


(7.21)  |i|>k(x,y,p)  |  <  for  all  (x.y)  e  R3  and  p  6  K. 

Proof  of  Lemma  7.5.  To  prove  (7.18)  note  that,  by  (1.31),  (1.36) 

(7.22)  Vx'?’P’X)  “  <2ir>’1  eip’x  a+<lpU>  ^<y*lpl*x> 

where  a+(p,A)  is  defined  by  (5.14).  The  continuity  of  0H(y, |p|,X)  on 
R*fl  follows  from  Corollary  2.2.  The  continuity  of  a+(|p|,A)  on  ft 
follows  from  (5.10),  (5.14)  and  the  assumed  continuity  of  the  phase 
function  0+(lJ,A).  Thus  the  continuity  of  on  R3  x  ft  follows  from 

(7.22) .  The  proofs  for  and  are  similar  and  will  not  be  given. 

The  continuity  of  <J^,  k  >  1,  follows  from  Corollary  6.5. 
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To  prove  (7.19)  for  note  chat  (7.22)  and  Che  continuity  of 
a+(|p|»^)  Imply  chat  it  is  enough  to  prove  the  existence  of  a  constant 
Mg  such  that 

(7.23)  |$4(y, 1 P 1 ,A) |  5  Mg  for  all  y  €  R  and  (p,X)  €  K. 

Now  the  uniformity  of  the  asymptotic  estimates  (2.11)  on  the  compact 
sets  r4  of  Corollary  2.4  implies  that 

(7.24)  <My>  IpM)  •  ®*p  t-iy  <i_(|p|.*)Hi  +  o(i)J,  y  -► 
uniformly  for  (p,A)  €  K.  Hence,  there  exists  a  constant  yR  such  that 

(7.25)  |<My.  |pf\*>  I  i  2  for  all  y  <  -yR  and  (p,A)  S  K. 
Similarly,  using  the  relation 

(7.26)  My.y.A)  *  c41(y,A)  <f>j(y,li,A)  +  c42(u,A)  4>2(y,y,A) 

from  (4.1),  (4.2),  (4.3),  the  continuity  of  c41(u,A)  and  c42(y,A)  and 
the  uniformity  of  the  asymptotic  estimates  for  <t>2  when  y  -►  +°°,  one 
finds  that  there  exist  constants  yg,  such  that 

(7.27)  |<My»  IpI  »A)  |  <  Mg  for  all  y  >  y^  and  (p,A)  €  K. 

Finally,  the  continuity  of  $4(y, |p|»A)  on  R  «  SI,  which  follows  from 
Corollary  2.2,  implies  the  existence  of  a  constant  M£  such  that 

(7.28)  I <>4  (y.  |p|  .A)  I  <  M£  for  -yR  <  y  <  y£  and  (p,A)  e  K. 

Combining  (7.25),  (7.27)  and  (7.28)  gives  the  estimate  (7.23)  with 
Mg  -  Max  (2, Mg, Mg)  . 
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The  proofs  of  (7.19)  for  and  of  (7.20)  and  (7.21)  can  be 
given  by  the  same  method.  This  completes  the  discussion  of  Lemma  7.5. 

Proof  of  Theorem  7.1.  Consider  the  function  f+(p,A).  The 
absolute  convergence  of  the  integral  in  (7.1)  for  each  (p,A)  €  ft  follows 
from  (7.19).  To  prove  that  f+  €  c(ft)  let  (p0,A0)  e  ft  and  let  K  C  ft  be 
compact  and  contain  (p0>AQ)  in  its  Interior.  Then  by  Lemma  7.5 

(7.29)  |<J>+(x,y,p,A)  f(x,y)|  <  ^  |f(x,y)|  for  (x,y)  6  R3,  (p,A)  e  K. 

Hence,  the  continuity  of  f+  at  (p„,A0)  follows  from  (7.18)  and  (7.29)  by 
Lebesgue's  dominated  convergence  theorem.  The  continuity  of  f_,  fQ  and 
follows  by  the  same  argument.  This  completes  the  proof  of  Theorem  7.1. 
Relationship  of  A  to  A^.  As  a  preparation  for  the  proofs  of 
Theorems  7.2,  7.3  and  7.4  the  operator  A  will  be  related  to  A|p|  by 
Fourier  analysis  in  the  variables  x  6  R2 .  To  this  end  note  that  if 
u  €  X  then  Fubini's  theorem  implies  that  u(*,y)  €  L2(R2)  for  almost 
every  y  €  R.  Thus  if  F  :  L2(R2)  *►  L2(R2)  denotes  the  Fourier  transform 
in  L2(R2)  then  the  Plancherel  theory  implies  that 

(7.30)  u(p,y)  -  (Fu)(p,y)  -  L2(R2)-lim  (2tt)-1  [  e_iP*x  u(x,y)dx 

M-*00  ■’  |x|<M 

exists  for  almost  every  y  €  R  and 

(7.31)  [  |u(p,y)|2  dp  ■  [  |u(x,y)|2  dx  for  a.e.  y  e  R. 

'  R2  '  R2 


Another  application  of  Fubini's  theorem  gives 
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Lemma  7.6.  u  €  3C  if  and  only  if  u  ■  Fu  €  3C  and  the  mapping 
F  :  X  ■*  X  is  unitary.  In  particular 

(7.32)  lul^  -  lul^  for  all  u  e  3f. 

The  Fourier  transform  of  the  acoustic  propagator  A  will  be 
denoted  by  A.  Thus 

(7.33)  A  -  FaF"\  D(A)  -  F  D(A) . 

A  A 

A  more  detailed  characterization  of  D(A)  is  needed  to  relate  A  to  A.  < . 

I  P I 

It  will  be  based  on 

Lemma  7.7.  Let  u  €  X.  Then  DjU  €  3C  (j  ■  1,2)  if  and  only  if 
Pjfi(p,y)  e  3f  and 

(7.34)  F  DjU  »  Pj  F  u,  j  -  1,2. 

Similarly,  D^u  €  X  if  and  only  if  Dyu  6  JC  and 

(7.35)  F  Dy  u  »  Dy  Fu. 

Proof  of  Lemma  7.7.  The  distributional  derivatives  D^u,  D^u  may 
be  characterized  as  temperate  distributions  on  the  Schwartz  space  S(R3) 
of  rapidly  decreasing  testing  functions  [7].  S(R3)  is  mapped  onto 

itself  by  F.  The  proof  of  (7.34)  is  essentially  the  same  as  in  the 
standard  Plancherel  theory.  To  verify  (7.35)  note  that  the  distribution- 
theoretic  definition  of  D^u  €  X  is 

(7.36)  [  D  u(x,y)  $(x,y)dxdy  ■  -[  u(x,y)  D  d>(x,y)dxdy  for  all  4>  e  S(R3). 

;R3  y  jr3  y 
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Application  of  Parseval's  relation  gives 


(7.37) 


dpdy 


|  (Fu)(F 
JR3 


Dy<t>)dpdy  for  all  <p  <=  S(R3) . 


Now  for  <J>  €  S(R3)  it  is  easy  to  verify  that 
(7.38)  Dy$  -  DyF  <p  -  F  Dy<|>. 

Thus  (7.37)  is  equivalent  to 


(7.39)  (F  D  u)4>  dpdy  *  -  (Fu)  D  $  dpdy  for  all  $  6  S(R3) 

jR3  y  ;r3  y 


which  in  turn  is  equivalent  to  (7.35). 

A 

Application  of  Lemma  7.7  to  A  gives 

Lemma  7.8.  The  operator  A  is  characterized  by  the  relations 


(7.40) 

F  L‘(R3)  -  (u  | 

p,u,  p2u  and 

Dyu  are  in  3f}, 

(7.41) 

D(A)  »  F  L3(R3)  n  {Q  | 

Dy  (p’lDyu) 

-  | p |  2  u  €  Jf),  and 

(7.42) 

&  u  »  -c2  (pDy(p" 

'DyU)  -  |P|2 

u} ,  u  6  D(A)  . 

Proof  of  Lemma  7.8.  These  results  follow  from  application  of 
Lemma  7.7  to  the  definition  of  l£(R3),  D(A)  and  A  -  equations  (1.2), 
(1.8)  and  (1.9). 

Corollary  7.9.  For  all  u  e  D(A)  one  has 

u(p , • )  €  D(A|p|)  and 


(7.43) 
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(7.44)  (A  u)(p,‘)  -  A|p|  u(p,*) 
for  almost  every  p  €  R2. 

Corollary  7.9  is  an  immediate  consequence  of  Lemma  7.6,  Lemma 
7.8  and  Fub ini's  theorem. 

The  Sets  X^°m,  X\  X"  and  X"' .  The  following  subsets  of  X  will 
be  used  in  the  proofs  of  Theorems  7.2,  7.3  and  7.4. 

(7.45)  JC001”  ■  3f  n  (f  |  supp  f  is  compact), 

(7.46)  X'  -  F”10(R3)  -  {f  |  f  -  Ff  G  P(R3)), 

(7.47)  3T'- (f(x,y)  =  fj(x)f2(y)  |  f 2 €  T ^(R2)  ,f 2e jc(R)  ,supp  f2  compact), 

m 

(7.48)  X"'  -  span  X"  -  (f  -  £  aa  fa  I  aa  €  C,fa€  3f"}’ 

a*l 

In  (7.46)  and  (7.47),  P(Rn)  denotes  the  Schwartz  space  of  testing  func¬ 
tions  with  compact  support  [7].  The  sets  X'  and  X”’  are  linear 

submanifolds  of  X  which  are  dense  in  Jf.  Indeed,  it  is  well  known  that 
P(R3)  is  dense  in  L2(R3).  This  fact  implies  that  JfCOm  is  dense  in  X. 

The  denseness  of  X’  in  X  follows  from  that  of  P(R3 )  and  the  unitarity  of 
F.  The  denseness  of  Xn'  follows  from  the  fact  that  X  is  the  tensor 
product  of  L2(R2)  and  JC(R) . 

It  is  clear  that  each  of  the  sets  JCcom,  Jf'  and  Jf"'  is  a  subset 
of  X  n  Lj(R3).  Hence,  for  f  in  one  of  these  sets,  the  transforms  f+,  f„ 
and  f^  defined  by  (7.1)- (7. 3)  are  continuous  functions  by  Theorem  7.1. 

An  alternative  characterization  is  given  by 


L 


rjvni' 
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Lemma  7.10.  If  f  €  UC00111  u  X'  U  Jf"’  then 


(7.49) 


(7.50) 


(7.51) 


f+(p,X) 

f0(p>» 

fk(P) 


My* |p | ,A)  f (p.y)  c-2(y)  p-1(y)  dy, 

R  " 


» 

■My*  I p I  f(p»y)  c~2  (y)  p_1(y)  dy, 

JR 


|  MyTTpTT  f(p»y)  C_2(y)  p-1  (y)  dy,  k  >  1. 


Proof  of  Lemma  7.10.  Equations  (7 .49)-(7 .51)  follow  from 
(7.1)-(7.3)  on  substituting  the  definitions  (1.36)-(1.38)  of  the  normal 
mode  functions  and  carrying  out  the  x- integration.  These  operations  are 
justified  for  f  €  LX(R3)  by  Lemma  7.5  and  Fubini’s  theorem. 

Corollary  7.11.  If  f  €  JC'  and  f(x,y)  *  fx(x)  f2(y)  then 

(7.52)  f±(p,A)  *  f1(p)f2±(|p|,A)  =  [F  <*»  |p  |±f  ]  (p.A)  =  [*  )p  [±Ff  ]  (p.  A) , 

(7.53)  f  „  (p.  A)  =  f1(p)f2Q(  |pj  ,A)  =  [F  *|  i.fKp.A)  =  [<f|p|0  pfHp.A), 
(7 - 54)  ffc(p)  =*  fj(p)f2k(|p|)  =  [F  <J,|pjkf](p)  *  [<I>|p|kFf](p)  ‘ 

These  results  follow  immediately  from  Lemma  7.10  and  the  results 

of  §5. 

The  notation 

(7.55)  R(T,0  -  (T  -  O’1 

will  be  used  for  the  resolvent  of  T.  The  proofs  of  Theorems  7.2,  7.3 
and  7.4  will  be  based  on  Stone's  theorem  relating  R(A,C)  and  the 
spectral  family  of  A,  together  with  the  following  three  lemmas  relating 
A  and  A  |  p  j  . 
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Lemma  7.12.  Let  f  €  3C"'  and  let  *  R(A, £)f  or,  equivalently, 

=  R(A,?)f.  Then 

(7.56)  u^(p,y)  -  |  G|p| (y,y* ,C)  f(p.y')  c~2(y')  p-1(y')  dy' 

for  almost  every  (p,y)  €  R3  where  G^(y,y',i;)  is  the  Green's  function 
for  Ay. 

Proof  of  Lemma  7.12.  It  is  enough  to  verify  (7.56)  for  functions 
f  €  -7 Let  f(x,y)  =  fj(x)f2(y)  be  such  a  function  so  that  f(p,y) 

=  f1(p)f2(y).  Now  €  D(A)  and  hence  by  Corollary  7.9 

(7.57)  ((A  -  C)u?)(p,y)  *  ((A|p|  ~  C)u^)(p,y)  *  fj(p)f2(y) 
for  almost  every  (p,y)  €  R3 .  It  follows  that 

(7.58)  G?(p,y)  -  [R(A|p|,0?1(p)f2l(y)  *  f, (p) [R(A|p| ,5)f2 ] (y) 
which  is  equivalent  to  (7.56)  because 


(7.59) 


[R(A|p| 


G|p|(y,y\<;) 


f2(y’)  c-2(y')  p  1 (y')  dy'. 


Lemma  7.13.  Let  f(x,y)  -  fx(x)  f2(y)  and  g(x,y)  »  gt(x)g2(y)  be 
elements  of  jc"  and  let  £  *  X  +ie  with  X  €  R  and  e  >  0.  Then  for  all 
y  €  R  one  has 


f  (f,[R(A,C)  -  R(A.O]g)dX 

J-i 

(7.60) 

-  |  fj(P>  Si  (P)  |  (f2,[R(A|p|,5)  -  R(A|p|,C)]g2)dXdp. 
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Proof  of  L^nima  7.13.  The  Plancherel  theory  implies  that 
(7.61)  (f,[R(A,5)  -  R(A,?)]g)  -  (f,[R(A,?)  -  R(A,?)]g). 

Combining  this  and  (7. 58)  gives,  by  Fubini's  theorem 


(7.62)  (f,[R(A,0-R(A,0)g) 


g1(p)(f2,[R(A|pj,0-R(A|pj,C]g2)d 


Now  a  standard  estimate  for  the  resolvent  of  a  selfadjoint  operator 
[8,  p.  272]  implies  that 

(7.63)  |fj(p)  |l(p)(f2,[R(A|pj,5)-R(A|p|,C)]g2)|  <  f  I  f  i  (p)g,  (p)  I'  f2» ■  82l 

Thus  integrating  (7.62)  over  “1  <  X  <  y  gives  (7.60)  by  Fubini's  theorem. 

Lemma  7.14.  The  spectral  family  {JI(p> }  satisfies  relation 
(7.17)  of  Theorem  7.4  for  all  f,g  €  3C"' . 

Proof  of  Lemma  7.14.  Stone's  theorem  in  its  general  form  is 
[16,  p.  79] 


lim  -Jr  [  (f,[R(A,A+ie)  -  R(A,  A-ie) ]g)dA 

e-HM-  1  Ja 

(7.64) 

-  (f,  [11(b)  +  Il(b-)  -  n(a)  -  n(a-)]g). 


In  the  present  case  o(A)  C  [0,®)  and  (7.64)  implies 


(7.65)  ^(f,[II(y)  +  II(y-)]f)  -  -sjr  lim  f  (f ,  [R(A,A+ie)  -  R(A,A-ie)  ]f)dA. 

£-*0+  J-1 


Moreover,  since  Il(y-)  -  lim  II(y  -  6)  and  lim  II((y  -  5)-)  ■  II(y— ) , 

5-KH-  6-*0+ 

(7.65)  implies 


•s*  .».*»» 
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1  fU 

(7.66)  (f,n(p-)f)  «  lim  11m  (f ,  [R(A,X+ie)  -  R(A,X-ie)  ]f)dX. 

1  6-*0+  e-KH-  J-i 

If  f(x,y)  *  fJ(x)f2(y)  €  3 C  then  combining  (7.60)  with  g  *  f  and  (7.66) 
gives 


(f,n(u-)f) 

(7.67) 

-^jlim^lim^  J  Jfx(p)|2  (f2,{R(A|p|,X+ie) -R(A|p|,X-ie)]f2)dXdp. 


Now  application  of  the  spectral  theorem  to  A^j  gives 

(7.68)  (f2,  [R(A|p|,X+ie)  -  R(A|p|,X-ie)]f2)  ■  J  (f  2>II|p  |  (dX’)f  2) 


It  follows  that 


(7.69) 


r 

j  (f 2, [R(A|pj , X+ie)  -  R(A|p| ,X-ie) ]f2)dX 

'  y-6 


<j  r 

JR 


2e 

(HW 


dX 


(f2,n|p|(dX')f2). 


Moreover, 


(7.70) 


(X 


-A’)2+e2  ^  -  In  (X-X^)2+ei  dX 


2tt 


for  all  X'  s  R  and  all  e  >  0.  Combining  (7.69)  and  (7.70)  gives 
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/•y-6 

(7.71)  J  (f  2, [R(A|p | ,X+ie)  -  R(A|p|,X-ie)]f2)dX  <  2tt  lf2l 


for  all  p  6  R2  and  all  y  >  0,  6  >  0  and  e  >  0.  In  addition  Stone's 
theorem  applied  to  A<  i  gives 


TTliii  11a  (f2,[R(A,  ,,X+ie)  -  R(A.  >,X-ie)]f  )dX 

<WH-  e-K)+  J-i  |P|  IPI 


(7.72) 


(f2,n|p|(y-)f2). 


Equations  (7.67),  (7.72)  and  the  estimate  (7.71)  imply,  by  Lebesgue's 
dominated  convergence  theorem, 


(7.73)  (f,n(y-)f)  •  J  2  |fj(p)|2  (f2,njp J (y-)f 2)  dp. 


It  follows  by  polarization  that 


(7.74)  (f,n(y-)g)  *  |  2  ?,(p)  gj(P)  (f2,II|p|(y-)g2)  dp 


for  all  f,g  e  X4' .  The  same  argument  applied  to  (7.65)  gives  the  relation 


(7.75)  (f,[II(y)  +  n(y-)]g)  -  j  *  ?2(p)  g2 (p) (f2»  [n j p |  (y)  +  n|p( (y-)  ]g2)dp 


Subtracting  (7.74)  from  (7.75)  gives 


(7.76)  (f.il(y)g)  «  |  2  ft(p)  g1(p)(f2,IIjpj  (y)g2)dp 
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for  all  f,g  S  X".  To  prove  the  relation  (7.17)  for  f,g  e  JC"  note  that 
the  construction  of  11^  given  by  (5.27)  implies,  by  polarization 


2  | p | (U)g2) 

(7.77) 


-  f  H(W  -  X)[£  (|p|,A)  g2.(|p|,X)  +  f2.(|p|,X)  g._(|p|,X)]dX 

JA(|p|)  2+2+  22 


f  ,  H(y  -  X)  f20(|p|,X)  g2#(|p|,X)  dX 
JA0(  P  ) 


N(|p|)-1  - - 

+  l  H(p  -  Xk(|p|))  f2k(|p|)  glk(|p|) 


for  all  p  €  R2  such  that  |p|  >0,  where 


(7.78) 


A ( | p | )  -  (X  |  c2(-®)|p|2  <  X},  and 
A0 (  1  p | )  -  {X  |  c2 (“»)  | p |  2  <  X  <  c2(-«)|p|2}, 


Substituting  this  into  (7.76)  and  recalling  the  definitions  of  ft,  ft0  and 
(k  >  1)  gives  (7.19)  for  f,g  €  K".  The  relation  extends  immediately 
to  all  f,g  €  3f"'  by  linearity.  This  completes  the  proof  of  Lemma  7.14. 

Proof  of  Theorem  7.2.  Let  f  S  3f  and  M  >  0  be  given.  Then  since 
3C"'  is  dense  in  3C  there  exists  a  sequence  {gn)  in  Jf"*  such  that  gn  -*■  fjj 
in  3C.  Note  that  since 

rM  f 

,fM  "  gn'K  "  J  Jr2  I *<*»*>  "  8n(x.y)|2  c"2(y)  p~*(y)dxdy 

+  f,  [  ,  |gn(*»y)|2  c“2(y)  p"‘(y)  dxdy 
*  | y | >M  'R2  n 


(7.79) 
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It  aay  be  assumed  Chat  g^(x,y)  ■  0  for  |yj  >  M.  Now  Leona  7.14  Implies 
chac  Parseval's  relation 


(7.80) 


Igljf  “  Igl^-  *8+l  2  +  ■SJ2  +  ^  lgkl2 


holds  for  all  g  6  JC'" ,  where  g  ■  (g+.g_*g0»g1 . ***) •  On  applying  (7.80) 
to  the  differences  gQ  -  gm  it  is  found  that  gn  -  (gn+>8n_»gn0 *8nl *  ”  *) 
is  a  Cauchy  sequence  in  JC.  Hence  there  exists  a  limit 

(7.81)  lim  g  -  h  -  (h+,h_,h#,hj,*»*)  e  JC, 

n-*°° 

since  JC  is  a  Hilbert  space.  To  complete  the  proof  of  Theorem  7.2  it  will 
be  enough  to  show  that 


(7.82) 


^(p.A)  -  h+(p,A)  for  a.e.  (p,A)  6  (2, 
fM0(p,X)  -  h9(p,A)  for  a.e.  (p,A)  e  £2fl, 
f^p)  *  hk(p)  for  a<e’  P  6  k  -  1’ 


This  clearly  implies  (7.6)  since  h  €  JC.  Moreover,  since  Hilbert  space 
convergence  implies  convergence  of  the  norms,  the  relation  (7.80)  for 


gn  €  JC"'  implies 


Nn-1 

(7.83)  lfMl£  -  lim  |gni£  -  lim  lgnl£  -  lhl£  -  lh+l2  +  lh_|2  +  £  Ih^l 

nr*°°  n-**  k“0 

which  is  equivalent  to  (7.7)  when  (7.82)  holds. 

Relation  (7.82)  will  be  proved  for  ?M| .  The  proofs  for  the 

remaining  cases  are  entirely  similar.  To  prove  that  f^+(p,X)  -  h+(p,A) 

for  a.e.  (p,X)  in  $2  note  that  if  K  is  any  compact  subset  of  (2  then 

f^  S  C(K)  C  L2(K)  by  Theorem  7.1  and 
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(7.84) 


Now  by  Lemma  7.10 


1  *M+  '  VLj(K)  “  ^  1  ^M+  "  ®n+*L2(K) 


(7.85)  -  5Trt.(p.A>  -  J*1  i|»+(y,  I p | ,X)  [fM(p,y)  -  gn(p,y) ]c-2 (y)p_1 (y)dy. 


Hence  by 


7.5  and  Schwarz's  Inequality 


(7.86) 


~  «iH.(p.A)  |  <  Mj^  ^  |^,(p.y)  ~  8n(p,y)  |  c*2(y)p_1(y)  dy 

-M 

<  c*2(y)p"l(y)dy|  |j  jf^p.y)  -  gn(p,y)  1 2c"2(y)p“l(y)dy 


for  all  (p,A)  €  K.  It  follows  that  there  is  a  constant  C  •  C(K,M)  such 
that 


1/2 


(7.87) 


l?ttf  ’  ^tH-'l^K)  -  C  ,gM  -  *n'3f  *  C  ,fM  “  gn'3f 


Since  g  -►  fu  in  X,  (7.87)  implies  that  the  limit  in  (7.84)  is  zero  and 
n  h 

hence  f^  (p,A)  «  h+(p,A)  for  a.e.  (p,A)  e  K.  This  completes  the  proof 
since  KC  SI  was  an  arbitrary  compact  set. 

Proof  of  Theorem  7.3.  To  prove  that  { f^.}  is  a  Cauchy  sequence 
in  X  for  M  -*■  ®,  let  M  >  0  and  N  >  0  be  arbitrary  numbers  and  let  {g^j}, 
{g^}  be  sequences  in  xm  such  that  g^  -►  fM,  g^  -*•  fjj  in  X.  Then,  as 

-wM  ““ 

proved  above,  g  ■*  fu  and  g  -*•  f„  in  X  and  Parseval's  relation  (7.80) 
n  ti  n  n 

M  N 

holds  for  g  -  g  .  Passage  to  the  limit  n  00  gives 


(7.88) 


1  fM  “  V*  "  '*M  “ 
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which  Implies  that  {fM)  is  a  Cauchy  sequence  because  -►  f  in  Jf  when 
M  -►  ®.  Finally,  one  gets  (7.13)  for  arbitrary  f  €  Jf  by  passage  to  the 
limit  in  (7.9). 

Proof  of  Theorem  7.4.  It  will  be  enough  to  prove  (7.17)  for 
f  *  g  €  Jf  since  the  general  case  then  follows  by  polarization.  Now  by 
Lemma  7 . 14 


(g,II(y)g)  -  R(y-X)([g  (p,X) I2  +  |g_(p,X)l2)dpdX 

A  ' 


+  j  H(y-X) |g0(p,X) | 2  dpdX 
No-1  r 

+  l  H(y  -  X.  ( |  p  | ) )  |gk(p)|2dp 
k-1  \ 


for  all  g  €  X”' .  Let  f  €  Jf,  M  >  0  and  let  (gn)  be  a  sequence  in  Jf"'  such 
that  gn  fM  in  X.  Then  it  follows  from  the  proof  of  Theorem  7.2  that 
^  -*•  fM  in  3f.  Replacing  g  by  gQ  in  (7.89)  and  making  n  -►  »  gives  (7.89) 
with  g  *  fM.  If  N0  ■  +»  then  passage  to  the  limit  is  justified  because 
the  right-hand  side  of  (7.89)  is  majorized  by 


No-1 

(7.90)  I  g .  I  2  +  I  g  *  2  +  I  lgj2<«. 

k-0 


Thus  (7.89)  is  valid  with  g  ■  f ^  where  f  €  Jf  and  M  >  0  are  arbitrary. 
Making  M  -*•  ®  and  repeating  the  above  argument  gives  (7.89)  with 
g  »  f  6  Jf,  by  Theorem  7.3. 


Another  proof  may  be  obtained  by  noting  that  the  left-hand  side 
of  (7.89)  is  a  bounded  quadratic  form  on  Jf,  while  the  right-hand  side  is 


a  bounded  quadratic  function  of  g  ■  'i'g  because  of  the  majorization  by 
(7.90).  Thus  (7.17)  follows  from  the  boundedness  of  4*  and  the  fact  that 
(7.89)  holds  for  g  in  the  dense  set  . 
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§8.  Normal  Mode  Expansions  for  A. 

The  normal  mode  expansions  for  the  acoustic  propagator  A  that 
are  the  main  results  of  this  report  are  formulated  and  proved  In  this 
section.  The  starting  point  Is  the  representation  of  the  spectral 
family  of  A  given  by  Theorem  7.4.  The  main  result,  Theorem  8.8,  shows 
that  the  family  is  a  complete  orthogonal  family  of 

normal  modes  for  A.  theorem  8.9  shows  that  it  provides  a  spectral 
representation  of  A.  These  results  are  shown  to  imply  that  the  families 
and  defined  in  §1,  are  also  complete 

orthogonal  families  of  normal  modes  for  A  and  provide  alternative 
spectral  representations. 

The  basic  representation  space  for  A  associated  with  the  family 
»***}  is  the  direct  sum  space 

No-1 

(8.1)  X  -  La<8)  +  L,(R>  +  I  L2(0.) 

k-0  K 

introduced  in  §7.  Theorem  7.3  associates  with  each  f  €  JC  an  element 
f  €  JC.  The  Parseval  relation  (7.13)  implies  that  the  linear  operator 

(8.2)  ¥  :  JC  -*■  JC 
defined  by 

(8.3)  ¥  f  -  f  for  all  f  e  JC 


is  an  isometry;  i.e.. 


(8.4) 


IV  fl^  -  for  all  f  €  JC. 


The  principal  result  of  this  section  is 


(8.5) 


(8.6) 


Theorem  8.1.  The  operator  V  is  unitary;  i.e.. 


V*  V  =  1  in  JC,  and 


V  V*  -  1  in  3f. 


Relations  (8.5)  and  (8.6)  generalize  the  completeness  and 
orthogonality  properties,  respectively,  of  the  eigenfunction  expansions 
for  operators  with  discrete  spectra.  Relation  (8.5)  is  equivalent  to 
(8.4)  and  thus  follows  from  Theorem  7.3.  Relation  (8.6)  is  shown  below 
to  follow  from  the  unitarity  of  the  operator  associated  with 
(Theorem  5.6). 

The  completeness  relation  (8.5)  implies  that  every  f  €  JC  has  a 
normal  mode  expansion  based  on  the  family  »’i'i » ^  •  The 

orthogonality  relation  (8.6)  implies  that  the  space  JC  is  isomorphic  to 
JC  and  thus  provides  a  parameterization  of  the  set  of  all  states  f  €  JC 
of  the  acoustic  field.  These  Implications  of  Theorem  8.1  will  be  devel¬ 
oped  in  a  series  of  corollaries. 

The  normal  mode  expansion  will  be  based  on  the  linear  operators 


(8.7) 


(8.8) 


V+  :  JC  L2  (0) 


\  :  K  L*  (fV  *  0  -  k  <  No  * 


defined  for  all  f  €  JC  by 


(8.9) 


’"■***  '■*'  •*  — im.  btwhimw 
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(8.10)  \  f  -  fk,  0  <  k  <  N0, 

where  l+,  f^  are  defined  as  in  Theorem  7.3.  It  is  clear  that 

(8.11)  V  f  -  ('^f.'f.f.V’V'***) 

for  all  f  e  3f  and,  by  (8.4), 

■  Nor1 

(8.12)  i  4*  f  I  2  +  IY  fl2  +  l  0  f  ■  2  -  «fl2. 

k-0 

In  particular,  each  of  the  operators  V+,  4^  is  bounded  with  norm  not 
exceeding  1.  The  normal  mode  expansion  for  A,  in  abstract  form,  is 
given  by 

Corollary  8.2.  The  family  {4'+,4'_,4'|),4'l,«»«}  satisfies 

No-1 

(8.13)  1  -  4'*  V ,  +  V*  4'  +  l  V 

k=0  k 

where  1  is  the  identity  operator  in  3C  and  the  series  in  (8.13)  converges 
strongly. 

It  will  be  shown  that  (8.13)  is  equivalent  to  the  completeness 
relation  (8.5).  The  orthogonality  relation  (8.6)  will  be  shown  to 
be  equivalent  to  the  relations  described  by 

Corollary  8.3.  The  family  {4'+,4,_,’i'0,'?1  ,•••}  satisfies  the 


relations 

(8.14) 

»±  -  1  in  L2(fl), 

(8.15) 

\  -  1  in  L20y,  0  <  k  <  N0. 
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In  addition,  H»+  Y*  -  0,  I'*  «  0,  ^  4<*  «  0  and  ?k  4fJ  -  0  for 

all  k  and  2,  #  k  such  that  0  <  k,  2,  <  Nfl. 

Relations  (8.14),  (8.15)  imply  that  each  of  the  operators  ¥  , 

4*_,  ^  (0  <  k  <  N0)  is  partially  isometric  [8,  p.  258].  It  follows  that 
the  operators  in  JC  defined  by 


P  »  v* 

+  1  +  1  -t 


(8.16) 


Pk  -  0  <  k  <  No. 


are  orthogonal  projections  in  JC  onto  subspaces 


K  -  P+  X 


(8.17) 


\  -  Pk  Jf,  0  <  k  <  N0 . 


Combining  this  with  Corollaries  8,2  and  8.3  gives 

Corollary  8.4.  {P+,P_,P0 ,Pj , • • • }  is  a  complete  family  of 

orthogonal  projections  in  3C;  i.e.,  the  range  spaces  defined  by  (8.17) 
are  mutually  orthogonal  and 


(8.18) 


space 


(8.19) 


N0-l 

1  -  p+  +  p.  +  i  pt. 

k-0 


The  spaces  (8.17)  are  subspaces  of  JC  and  hence  the  direct  sum 


No-1 

X.  +  JC  +  I  JC. 
k-0 


may  be  identified  with  the  set  of  all 


8.2i 


f 


E 


‘  J 


Ln 


0) 


k=l 


in  3C  such  that 

(8.21) 


f+  e  3f+,  fk  G  \  for  0  <  k  <  N0 


and 


No-1 


(8.22) 


i f+a 2  +  8 f j 2  +  l  « f.i 2  < 


k=0 


With  this  convention.  Corollary  8.4  implies 

Corollary  8.5.  K  has  the  decomposition 


(8.23) 


N0-l 

jc  =  jc,  +  x  +  l  jc 

k=0 


The  definitions  of  the  operators  4'+,  4^  and  equations  (7.14)- 
(7.16)  imply 

Corollary  8.6.  The  operators  ,  4^  have  the  representations 


(8 


.24)  (V+£)(p,\)  -  L2 (ft)-lim  j  ^+(x,y,p, A)  f (x,y)c-2 (y)p-1 (y)dxdy, 

M-*»  1 R^  " 


(8.25)  (4’0f)(p,A)  *  L2(n0)-lim  [  i|;0  (x,y,p,A)  f(x,y)c~2(y)p-1  (y)dxdy, 

M-«o 


(8.26)  (U)(p)  -  L2(fik)-lim  f  \(x,y,p)  f  (x,y)c”Z  (y)p-1  (y)dxdy, 


for  1  <  k  <  N0  where  R^  ■  R3  n  {(x,y)  |  |x|  <  M,  |y|  <  M> , 
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Of  course,  the  family  of  sets  (R^  |  M  >  0}  can  be  replaced  by 
any  family  {K^  |  M  >  0}  of  compact  sets  whose  characteristic  functions 
tend  to  1  almost  everywhere  in  R3  when  M  -*■  ».  Corollary  8.6  implies  a 
similar  representation  for  the  adjoint  operators  ¥*,  f*.  It  is 
formulated  as 

Corollary  8.7.  The  operators  4'*,  have  the  representations 

(8.27)  (<t'*g+)(x,y)  »  3f-lim  [  \p+(xfy,p,\)  g+(p,X)dpdA, 

~  “  M-*°°  J£r  ~ 


(8.28) 


(I'ogoXx.y) 


K^lim 


*Po  (x,y,p,A)  g0  (p,A)dpdA, 


(8.29) 


a*gk)^.y) 


JC-lim 


\(x,y,p)  gk(p)dp. 


for  1  <  k  <  No  where  8**  and  ^  (0  <  k  <  N0)  are  families  of  compact 
subsets  of  (2  and  (2^  whose  characteristic  functions  tend  to  1  almost 
everywhere  in  (2  and  (2^,  respectively,  when  M  *►  00 . 

By  combining  Theorem  8.1  and  Corollaries  8.2,  8.6  and  8.7  the 
following  explicit  formulation  of  the  normal  mode  expansion  if  obtained. 

Theorem  8.8.  Every  f  €  X  has  a  representation 


No-1 

(8.30)  f(x,y)  -  f+(x,y)  +  f_(x,y)  +  J  ffc(x,y), 

k»0 


convergent  in  Jf,  where  f+  6  X+  and  f^,  0  <  k  <  N0 ,  are  given  by 


(8.31) 


f±(x,y) 


JC-lim 


<P±(x,y,p,A) 


f+(p»X)dpdA, 
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(8.32) 


(8.33) 


f0(x,y)  =  3f-lim  ^0(x,y,p,A)  f0(p,A)dpdA, 
M-*»  1 


fk(x,y)  =  3(111,11  m  \(*.y.p)  fk(p)dP. 


M-*» 


for  1  <  k  <  N0  and  f  +  ,  are  defined  by  (7 .14)-(7 .16) .  Conversely,  if 
f  =  (f+,f_, f 0 ,f j , • • •)  is  any  vector  in  3f  then  (8.30)-(8.33)  define  a 
vector  f  €  3C  such  that  f+,  ?k  are  related  to  f  by  (7.14)-(7.16) . 

Theorem  8.9.  The  unitary  operator  4*  defines  a  spectral  repre¬ 
sentation  for  A  in  3f  in  the  sense  that  for  all  f  €  D(A)  one  has 


(8.34)  (4*±A  f)(p,A)  -  X(V±f)(P,X)  *  X  f±(p,A), 


(8.35)  (V  f) (p,X)  -  X(Y0f)(p,X)  -  X  f0(p,X), 


(8.36)  (\A  f)(p)  =  Xk(|p|)(4'kf)(p)  -  Ak(|p|)  fk(p), 
for  1  <  k  <  Nfl . 

Corollary  8.10.  The  complete  family  of  orthogonal  projections 
{P+,P_,Pa,Pj ,•••}  reduces  A;  i.e., 

(8.37)  P+A  C  A  P±,  PkA  C  A  Pk 
for  0  <  k  <  N0  and  if 


(8.38) 


A+  *  A  P+  »  P+  A  P+ 

Ak  -  A  Pk  -  Pk  A  pk.  0  <  k  <  Na, 


i 


126 


denote  the  parts  of  A  In  JC+  and  then 


No-1 

(8.39)  A  «  A ,  +  A  +  I  A. 

k=0 


The  Family  {<t>4-(x,y,p,q)  |  (p,q)  €  R3  -  N}.  It  will  now  be  shown 
how  the  normal  mode  expansion  of  Theorem  8.8  can  be  reformulated  in 
terms  of  the  family  ,iJj2,  •  •  • } .  To  begin  equation  (1.59)  for  the 

normalizing  factor  c(p,q)  will  be  verified. 

First,  recall  that  the  normalizing  factors  a+(y,A),  a0(y,A)  for 
ip+f  i|>0  are  related  to  the  factors  c+(y,A),  c0(y,A)  in  their  asymptotic 
forms,  equations  (1.33)-(1.35) ,  by 

(8.40)  c+  =■  T"1  a+,  c0  *  T"1  a0; 

see  (4.8),  (4.14),  (4.21).  Combining  (8.40)  with  equations  (4.10), 
(4.16),  (4.23)  for  T±,  T0  and  equations  (5.10),  (5.14),  (5.19),  (5.21) 
defining  a+,  a„  gives 


(8.41)  c±(y,A)  -  (4¥p«a,)  >  c0(y,A)  *  (47r  JJ^vU).  * 


i0  1.0 

provided  that  the  phase  factors  e  -,  e  0  are  defined  by 


ji8+(y,A)  ,  i[(p  d>2]/|  t4,3<|,2]  |,  A  €  A(y) , 


(8.42) 


i9o(y,A)  .  i[#3*2]/|[*3#2]|t  x  €  A0(y). 


On  combining  (8.41)  and  the  definition  of  4>+(x,y,p,q) ,  equations  (1.53) 
and  (1.54),  and  the  asymptotic  forms  (1 .33) — (1 . 35)  for  i|>±,  \p0  one 
obtains  the  asymptotic  forms  (1.57),  (1.58)  for  <$>+  with  the  normalizing 
factor  c(p,q)  defined  by  (1.59). 


« 
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To  derive  the  normal  mode  expansion  for  the  family  •'i'a * *  *  * ) 

from  that  for  ic  is  clearly  sufficient  to  restrict 

attention  to  the  subspace 


(8.43) 


JCf  =  JC+  +  JC_  +  JC0 


and  the  corresponding  orthogonal  projection 


(8.44) 


Pf  -  P+  +  P_  +  P„. 


Thus  if  h  €  jf  and  h^  =  P^h  =  h+  +  h_  +  h„  then  Theorem  8.8  implies  that 


(8.45) 


hf(x,y)  -  ^.(x.y.p.X)  h  .(p,A)dpdA  +  <f»0(x,y,p,X)  h  (p.X)dpdX 

r  jq  Jn0 

+  ^_(x,y,p,A)  h  (p,A)dpdA 

JQ 


where  the  integrals  converge  in  X.  Changing  the  variables  of  integration 
in  the  three  integrals  by  means  of  the  mappings  X+,  X„  and  X_, 
respectively,  and  using  the  definition  (1.53),  (1.54)  of  <j>+  gives 


(8.46) 


hf(x,y)  =  |  ^(x.y.p.q)  h+(p,X)c(®)  (2q) 1/2  dpdq 

+  [  <t>+(x,y,p,q)  h0(p,X)c(®)(2q)  1/2  dpdq 

JC„ 

+  [  $+(x,y,p,q)  K_(P, X)c(-“»)  (2 1 q  |  )  1/2  dpdq 

c- 

M  ,  <f+(x,y,p,q)  h.(p,q)dpdq 

J r3-n  *  + 


where  X  ■  X(p,q)  is  defined  by  (1.56)  and 


128 


(8.47) 


A 

h+(p,q) 


(2q)l/2c(»)  h+(p,A(p,q)) ,  (p,q)  €  C+, 
(2q)l/2c(“)  h0(p,A(p,q)),  (p,q)  €  C„, 

(2 1 q | )  1/2  c(-<*>)  h_(p,X(p,q)) ,  (p,q)  €  C_. 


It  is  easy  to  verify  by  considering  the  three  cones  C+,  Co  and  C 
separately  that 


(8.48)  h+(p,q)  -  4>+(x,y,p,q)  h(x,y)  c“2(y)  p_1(y)dxdy 

'R3 

where  the  integral  converges  in  L2(R3).  Moreover,  it  can  be  verified  by 
direct  calculation,  using  the  Parseval  relation  of  Theorem  7.3,  that 

(8.49)  *  *N-*L2(R3) 

for  all  h  €  X.  These  considerations  suggest 

Theorem  8.11.  For  all  h  £  X  the  limit 


(8.50)  h  (p,q)  -  L2(R3)-lim  <J>.  (x,y,p,q)  h(x,y)c-2 (y)p-1  (y)dxdy 

K-*° 

exists.  Moreover,  the  mapping  :  X  L2(R3)  defined  by  $+h  *  fi+  is  a 
partial  isometry  such  that 

(8.51)  *+  $*  -  1  and  \  -  Pf, 
and  the  adjoint  mapping  h^  «  $+  h+  is  given  by 


(8.52)  hf(x,y)  -  JC-lim  4>+(x,y,p,q)  h.(p,q)dpdq. 

1  M-«c  J  (R3-N)m 

Finally,  <t+  is  a  spectral  mapping  for  A  in  the  sense  that  for  all 
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h  €  D(A)  one  has  j 

I 


(8.53) 


(4>+A  h)  (p,q)  -  X(p,q)  4>+  h(p,q) 


where  X(p,q)  is  defined  by  (1.56). 

Note  that  Theorem  8.11  is  simply  a  reformulation  of  Theorem  8.8 
and  not  a  new  theorem. 


The  Family  {<fr_(x,y,p,q)  I  (p,q)  g  R3  -  N).  The  analogue  of 
Theorem  8.11  for  the  family  {<{>_}  will  be  formulated  as 
Corollary  8.12.  For  all  h  €  X  the  limit 


(8.54)  h  (p,q)  -  L2(R3)-lim 

M-+® 

exists  and  the  mapping  :  X  L2(R3)  defined  by  $_h  *  h_  is  a  partial 
isometry  such  that 

(8.55)  4>*  -  1  and  <t>*  -  Pf. 

Moreover, 


4>_(x,y ,p ,q>  h(x,y)c-2(y)p-1  (y)dxdy 


(8.56)  hf(x,y)  *  JC-lim  $  (x,y,p,q)  h  (p,q)dpdq. 

1  M-*»  /rm 

Finally,  for  all  h  €  D(A) ,  one  has 

(8.57)  ($_A  h) (p, q)  -  X(p,q)  4>_h(p,q). 

These  results  are  direct  corollaries  of  Theorem  8.11.  This 
follows  from  the  observations  that  f(p,q)  f(-p,q)  defines  a  unitary 
transformation  in  L2(R3)  while  f  -*■  f  defines  a  unitary  transformation  in 
both  30  and  L2  (R3)  . 


ii 


I 
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This  completes  the  formulation  of  the  results  of  §8.  The  proofs 
will  now  be  presented. 

The  proofs  of  Theorem  8.1  and  Corollaries  8.2  and  8.3  will  be 
based  on  a  lemma  concerning  bounded  linear  operators  from  a  Hilbert 
space  into  a  direct  sum  space.  To  formulate  it  let  X  and  (k  €  N) 
denote  Hilbert  spaces,  where  N  is  a  finite  or  denumerable  set,  and 
define 


(8.58)  If-  l  JC. 

k£N 

Elements  of  X  will  be  written  g  ■  (gk)  where  g^e  ^  and 

(8.59)  Igl2  -  l  IgJu  <  “ 

kSN 

if  II  *  I  k  is  the  norm  in  3^.  If 

(8.60)  B  :  X  +  X 

is  a  bounded  linear  operator  then 

(8.61)  Bf  -  { (Bf )fc)  -  (Bkf} 
where 

(8.62)  Bk  :  X  * 

is  a  bounded  linear  operator.  With  this  notation  one  has 

Lemma  8.13.  The  adjoint  operator  B*:  X  X  is  given  by 

(8.63)  B*g  -  l  B*  gk,  g  -  {g.  }  €  X, 

kSN 


« 


* 
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where  Che  series  converges  in  Jf.  Hence 

(8.64)  B*B  -  l  B*  B,  , 

k€N 

where  the  series  converges  strongly,  and 

(8.65)  B  B*g  -  |  Bk  B*  gj|,  g  *  {gk>  €  X, 

where  the  series  converge  in  3K^. 

Proof .  If  N  is  finite  then  (8.63)  follows  directly  from  the 
definition  of  B*.  If  N  is  denumerably  infinite  then  (8.63)  holds  for 
all  g  €  X  with  finitely  many  non-zero  components  gk.  But  any  g  6  Jf  can 
be  approximated  in  X  by  such  vectors  and  B*  is  bounded.  The  convergence 
in  X  of  the  series  in  (8.63)  follows.  Equation  (8.64)  follows  on 
applying  (8.63)  to  the  vector  (8.61).  (8.65)  follows  on  applying  B  to 

the  vector  (8.63). 

Corollary  8.14.  B*B  *  1  in  X  if  and  only  if 

(8.66)  l  B*  B.  -  1  in  X, 

kfEN  K 

the  series  converging  strongly. 

Corollary  8.15.  B  B*  ■  1  in  Jf  if  and  only  if 

(8.67)  Bk  B*  -  6jk  for  all  k,j  S  N 

where  6jk  *  1  if  j  *  k  and  6jk  ■  0  if  j  +  k. 

Proof  of  Corollaries  8.14  and  8.15.  (8.66)  is  an  Immediate 

consequence  of  (8.64)  and  B*B  ■  1.  (8.67)  follows  from  B  B*  ■  1  and 

(8.65)  on  taking  g^  -  (5jk  gk},  j  fixed. 


» 
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On  applying  Corollary  8.14  to  Che  direcc  sum  space  (8.1)  and  che 

operator  V  it  is  seen  that  (8.13)  of  Corollary  8.2  is  equivalent  to 

(8.5)  of  Theorem  8.1.  Similarly,  Corollary  8.15  implies  that  the 

relations  of  Corollary  8.3  are  equivalent  to  (8.6)  of  Theorem  8.1. 

Proof  of  Corollary  8.2.  It  was  shown  that  (8.5)  follows  from 

Theorem  7.3.  Thus  (8.13)  is  valid  by  Corollary  8.14. 

Proof  of  Corollary  8.3.  The  spectral  mapping 

»***)  for  \  is  unitary  for  a11  U  >  0  by  Theorem 

5.6.  It  follows  from  Corollary  8.15  that  the  analogue  of  Corollary  8.3 

holds  for  4f.  ;  see  (5.38).  It  will  be  shown  that  Corollary  8.3  follows 

U  ♦ 
from  these  relations.  For  brevity  only  the  relation  *P0  ■  1  will  be 

proved.  The  remaining  relations  can  be  proved  by  the  same  method. 

For  the  proof  of  'Pj  *  1  define 

(8.68)  ?0  -  ¥,  r1  :  X  ■*  L2(80)  • 

Then 

t 

(8.69)  V,  V*  -  ?* 

and  it  will  suffice  to  prove  that  ?(  ■  1.  The  following  lemma  will 

be  used. 

Lemma  8.16.  For  all  g  €  V(ii0)  one  has 

(8.70)  (?*  g) (p»y)  -  f  <l>0 (y,  |p|  A)  g(p,X)dX  -  [f?  ,  g(P,«)](y). 

jA0(|p|)  lPl 

Proof  of  Lemma  8.16.  Let  f  €  X'  -  F"1  P(R3).  Then  ?  €  P(R3) 


and  by  Lemma  7.10  one  has 
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($*  g,f)  -  (g.$,f)  -  (g.V>  =  <S.f0> 

(8.71) 

-  [  g(p,X)  [  ^0(y. |p| ,X)  f (p,y)c-2 (y)p“l (y)dpdX. 

JQ0  Jr 

Now  ^0(y,|p|,X)  is  defined  for  (y,p,X)  e  R  x  Q0 .  Extend  the  definition 
to  all  (y,p,X)  6  R  x  RJ  by 

(8.72)  'My»jp|  ,X)  -  0  for  y  e  R,  (p,X)  €  R3  -  , 


and  apply  Fubini's  theorem  to  the  integral  in  (8.71).  This  gives 


(8.73) 


My*  IpI  *x)g(p,X)  dx 

R 


f (p,y)c-2 (y)p-1 (y)dpdy, 


which  implies  (8.70)  because  P(R3)  is  dense  in  X  and  supp  i|»o  (y.  |p  |  . *) 

C  A0 ( | P | ) ,  by  (8.72). 

Proof  of  Corollary  8.3  (completed).  The  relation  ?0  f*  -  1  is 
equivalent  to 


(8.74) 


(?0  ?*f,g)  -  ($*f,?*g)  -  (f,g) 


for  all  f,g  €  L2(Q0).  Moreover,  since  f*  is  bounded  it  will  suffice  to 
verify  (8.74)  for  f  and  g  in  the  dense  set  P(fl„).  Now  Lemma  8.16  implies 
that  for  all  f,g  €  P(fi0)  one  has 


(?*f,$*g) 


I  3  t1p|of(p’-»  ['•,|p|,8(P.-)](y)dpdy 


(8.75) 


1 
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But  by  (5.38),  ^|p|0  ^*p | 0  ■  1  for  all  |p|  >0  and  hence  the  last 
Integral  equals 


(8.76) 


l(f(p.-),g(p,-))A#(|p,)dp  -  (f.8)L2(o  ) 


which  proves  (8.74).  This  completes  the  proof. 

Proof  of  Theorem  8.1.  It  was  shown  that  (8.5)  follows  from 
Theorem  7.3.  Relation  (8.6)  follows  from  Corollary  8.3,  by  Corollary 
8.15. 

Proof  of  Corollaries  8.4  and  8.5.  Corollary  8.4  follows 
immediately  from  the  definitions  (8.16)  and  Corollaries  8.2  and  8.3. 
Corollary  8.5  follows  from  Corollary  8.4. 

Proof  of  Corollary  8.6.  This  is  Just  a  restatement  of  Theorem 
7.3,  equation  (7.12). 

Proof  of  Corollary  8.7.  Equation  (8.27)  for  4^  will  be  verified. 
The  proofs  of  the  remaining  equation  are  similar.  It  is  clearly 
sufficient  to  verify  (8.27)  for  functions  g+  €  Lz (ft)  with  compact 
supp  g+  C  tf*.  If  is  such  a  function  and  f  €  JCc0,n  then 


« 


« 


♦ 


(f.Tjg)  -  <V’S+>  -  <*+.*♦> 


(8.77) 


^(x.y.p.A)  f(x,y)  c  (y)p  (y)dxdy 


!  f(x,y) 

'  R3 

'J'+(x,y,p,A)g+(p,X)dpdA 

. 

g+(p,X)dpdA 


c~2(y)p~1 (y)dxdy . 


This  relation  implies  (8.27)  for  g+  because  JCco°  is  dense  in  3C. 
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Proof  of  Theorem  8.8.  The  representation  (8.30),  (8.31)  follows 
from  Corollaries  8.2  and  8.7.  The  converse  follows  from  the  unitarity 
of  H',  Theorem  8.1. 

Proof  of  Theorem  8.9.  Only  equation  (8.35)  will  be  proved.  The 
other  equations  can  be  proved  by  the  same  method.  To  prove  (8.35)  let 
f  €  D(A),  g  €  P(ftQ)  and  note  that 

(8.78)  (V  f.g)  -  (A  f,V$g)  =  (A  f,¥*g). 

Now  if  h  *  4^[g  then  h  €  3C,  by  Corollary  8.3  and,  by  Lemma  8.16, 


(8.79) 


h(p,y) 


I  >P0(y.  |p|.A)  g(p,A)dA,  (p,y)  e  R3. 

jA0(|p|) 


A 

A  distribution-theoretic  calculation  of  A  h  gives 


(8.80)  A  h(p,y)  *  f  |p|  ,A)  A  g(p,A)dA. 

A0(  I P I ) 

A 

In  particular,  since  X  g(p,X)  €  L2(ft0),  one  has  A  h  £  3f  by  Corollary 
8.3  and  hence  h  e  D(A) .  Combining  this  with  (8.78)  gives 


<¥,A  f.g)  *  (A  f,h)  -  (f.A  h) 


(8.81) 


f  f(p»y) 

• 

j  *,(y, | P [ ,A)  X  g ( p , X ) dX 

JR3 

ja0(1pI) 

.  j 

c-2(y)p-1 (y)dpdy . 


Writing  k(p,A)  *  X  g(p,A)  €  P(fl0)  and  applying  Lemma  8.16  again  gives 


(8.82) 


<V  f,g)  =•  (f,4f’k)  =  (f.'Qc) 


«  (>Jf0f,k) 


*  f0(p.x> 


f0(p,X)  k(p,X)dpdX 
g(p,X)dpdX . 


This  implies  (8.35)  because  the  functions  g  €  P(£)„)  are  dense  in  L2(G0). 

Proof  of  Corollary  8.10.  It  is  only  necessary  to  verify  relations 
(8.37).  By  [8,  p.  530]  these  relations  are  equivalent  to 

(8.83)  p±n(ii)  =  n(u)p±  ,  Pkn(u)  =  n(y)pfc 

for  0  <  k  <  N0  and  all  y  €  R.  These  equations  are  immediate  consequences 
of  Theorem  7.4.  For  example  (7.17)  and  the  definition  of  P+  implies  that 

« 

(8.84)  (f,p+n(y)g)  -  (f,n(y)P+g)  -  J  H(y-X)  f+(p,X)  g+(F.,X)dpdX 


for  all  f,g  €  3C,  which  implies  P+II(y)  *  H(y)P+.  The  other  relations  are  I 

proved  similarly. 

Proof  of  Theorem  8.11.  Equation  (8.50)  can  be  verified  by 
applying  the  definition  (8.47)  to  h^  and  using  the  convergence  statement 
of  Theorem  7.3.  Relations  (8.51)  follow  from  (8.47)  and  Corollary  8.3 
by  direct  calculation.  Equation  (8.52)  can  be  verified  by  reversing  the 
steps  in  the  calculation  (8.45),  (8.46).  Relation  (8.53)  follows  from 
Theorem  8.9  and  equation  (8.47). 

Proof  of  Corollary  8.12.  This  was  indicated  immediately  after 
the  statement  of  the  Corollary. 
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§9 .  Semi- Inf inite  and  Finite  Layers. 

The  purpose  of  this  section  is  to  present  extensions  of  the 
preceding  analysis  to  the  cases  of  semi-inif inite  and  finite  layers  of 
stratified  fluid.  The  methods  and  results  are  entirely  analogous  to 
those  developed  above.  Therefore  the  presentation  emphasizes  the 
modifications  required  in  these  cases.  Proofs  are  indicated  briefly  or 
omitted. 

Semi-Infinite  Layers.  With  a  suitable  choice  of  coordinates  the 
region  occupied  by  the  fluid  is  described  by  the  domain 

(9.1)  =  { (x,y)  |  y  >  0}. 

The  acoustic  field  is  assumed  to  satisfy  either  the  Dirichlet  or  the 
Neumann  boundary  condition  on  the  boundary  of  R+,  Physically,  these 
conditions  correspond  to  the  cases  where  the  boundary  plane  is  free  and 
rigid,  respectively.  The  functions  p(y)  and  c(y)  are  assumed  to  be 
Lebesgue  measurable  and  satisfy 

(9.2)  0  <  Pm  <  p(y)  <  PM  <  <*>,  0  <  cm  <  c(y)  <  <  °° 

for  all  y  >  0  and 


(9.3) 


|p(y)  -  p(°°)  i dy  < 


|  c(y)  -  c(°°)  |  dy  <  °°. 


The  acoustic  propagator  A  defined  by  (1.2)  determines  self adjoint 
operators  A0  and  A1  in 


(9.4) 


X+  =  L2  (R^.,c-2  (y)p-1  (y)dxdy) 


137 


138 


corresponding  to  the  two  boundary  conditions.  The  domains  of  A0  and  A1 
are  subsets  of 

(9.5)  l£(A,r|)  =>  L^(Rp  n  {u  |  7*(p_1Vu)  <=  L2(R^)}. 

The  operator  A0  corresponding  to  the  Dirichlet  condition  is  defined  by 

(9.6)  D(A°)  =  L* (A,R^)  n  (u  |  u(x,0+)  =»  0  in  L2(R2)}. 

Sobolev's  embedding  theorem  for  L2(R^)  [1]  implies  that  the  boundary 
values  u(x,0+)  are  defined  in  L2(R2)- 

The  Neumann  condition  will  be  interpreted  in  the  generalized 
sense  employed  in  [ 15 ] ;  i . e . , 


(9.7) 


(V*(p-17u)v  +  p  *Vu 


Vv}dxdy  »  0 


for  all  v  €  L2(R^).  Thus  the  operator  A1  corresponding  to  the  Neumann 
condition  is  defined  by 


(9.8)  D(AX)  -  L2 (A,R_^)  n  (u  |  (9.7)  holds  for  all  v  6  l£(r|)}. 
The  operators  are  defined  by  A^u  =  Au  for  all  u  e  D(A^)  and  one  has 

(9.9)  Aj  -  Aj*  >  0,  j  =  0,  1. 


This  is  most  easily  proved  by  introducing  the  corresponding  sesquilinear 
forms,  as  in  §1,  and  using  Kato's  representation  theorem  [8,  p.  322]. 

The  spectral  analysis  of  A°  and  A1  may  be  based  on  the 
corresponding  reduced  propagators  A0  and  A1  in 


(9.10) 


3C(R+)  -  L2(R+,c"2(y)P_1(y)dy),  R+  -  {y  |  y  >  0). 


They  are  defined  by 


(9.11)  D(Aj)  =  L*(R+)  n  {tj;  |  (p"V)'  €  3C(R+)  and  1JKO+)  -  0} 


(9.12)  D(Aj)  -  L‘(R+)  n  |  (p"V)'  €  3C(R+)  and  (pip')(0+)  -  0}, 


(9.13)  aJ<|>  =  A^  for  <p  e  D(aJ),  j  -  0,1, 


and  one  has 


(9.14) 


4  j  a  aJ  A  3  n  1 

\  Ay  -  cm  y  »  3  u,i- 


The  results  of  §3  can  be  extended  to  these  operators.  Thus 


(9.15) 


Oc(A^)  *  Se(Ay)  *  [ c2 (°°) la2 ,°°)  , 


(9.16) 


a(A^)  n  (-oo>c2(oo)v2)  c  a0(A^), 


and  the  eigenvalues  in  this  interval  are  all  simple.  They  will  be 
denoted  by  Ajj.(y) ,  1  <  k  <  Nj(y)  <  +». 

Eigenfunctions  of  A^.  These  functions  will  be  denoted  by 
^(y,U),  1  <  k  <  N^(y).  They  are  uniquely  defined  by  the  conditions 
^(*,U)  €  D(A^),  8^(’’U),3f(R+)  = 

(9.17)  (Ay  -  Aj[(y))  4^(y,y)  -  0  for  y  e  R+,  and 


(9.18) 


♦J(0+,u)  "  0,  (p">*')(0+,y)  -  0. 


The  asymptotic  behavior  of  t^(y,y)  for  y  -*■  +“  is  given  by 


^(y,y)  ~  c^(y)  exp  {-y  q' (y,Aj|(y))} ,  y 


1 


(9.19) 


140 

where 

(9.20)  q’(y,A)  -  (y2  -  A  c-2(“))1/2  >  0. 

The  function  ^(y»U)  has  precisely  k  -  1  zeros. 

Generalized  Eigenfunctions  of  A^.  For  A  >  c2(°°)y2,  A^  has  a 
single  family  of  generalized  eigenfunctions  {^(‘.y.A)  [  A  >  c2(°°)y2}. 
They  are  determined  up  to  normalization  by  the  conditions 

(9.21)  (A^  -  A)^(y*y,A)  *  0  for  y  e  R+>  and 

(9.22)  iJ,0(O+,y,A)  -  0,  (p“V  ')  (0+,y,A)  «  0. 

Their  asymptotic  behavior  for  y  -*•  00  is  given  by 

(9.23)  ^(y.y.A)  -  cj(y,A)[e"iyq(U,X)  +  Rj(y,A)  ei7q(u»A)]t  y  ^ 
where 

(9.24)  q(y,A)  =  (A  c_2(»)  -  y2)1/2  >  0. 

The  normalizing  factors  cJ(y,A)  are  calculated  below. 

Generalized  Eigenfunctions  of  A~^ .  These  are  defined  by 

(9.25)  ^(x.y.p.A)  -  (2it)"1  eip‘x  (y,  | p|  ,A) ,  (p,A)  €  0, 

(9.26)  >4(x,y,p)  *  (21T)'1  eip‘x  «|£(y,  |p|),  k  >  1 

where 

(9.27) 


♦ 


t 


t 


« 


I 


(9.28) 


n  -  { (p*  A)  I  A  >  C 2  («»)  I P I  2 } 
\  -  (p  i  I P I  e  0k},  k  >  l 
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and  0^  ■  {y  |  N^(y)  >  k  +  l),  as  before.  The  wave- theoretic  interpreta¬ 
tions  of  the  eigenfunctions  (9.25),  (9.26)  may  be  derived  from  the 
asymptotic  forms  (9.19),  (9.23)  as  in  §1.  With  these  definitions  the 
following  analogue  of  Theorem  7.4  holds. 

Theorem  9.1.  The  spectral  families  {11^ (y)}  of  A-'  satisfy 


(f,IIj(y)f)  -  [  H(y-X)  |fj(p,X)|2  dpdX 

(9.29) 

n3-1  r  .  . 

+  l  H(y  -  X^(|p|))  |f^(p)|2dp 

k-1 


for  all  f  €  3f  where  N»  =  sup  (y)  and 

y>0 


(9.30)  ?j(p,X)  -  L2(fi)-lim  ^(x,y,p,X)  f (x,y)c~2(y)p_l (y)dxdy, 

M+08  J  o  1  |x|<M 


(9 


4 

^(x,y,p)  f (x,y)c-2(y)p-1(y)dxdy. 

1  x  |  <M 


Relation  to  the  Infinite  Layer  Problem.  Theorem  9.1  can  be 
derived  by  the  method  employed  for  the  infinite  layer  problem  in  the 
preceding  sections.  However,  it  can  also  be  deduced  directly  from 
Theorem  7.4.  To  this  end  one  extends  p(y),  c(y)  to  all  y  6  R  as  even 
functions : 


(9.32)  p(-y)  -  p(y)»  c(-y)  ■  c(y) ,  y  >  0. 

Then  it  follows  from  (9.2),  (9.3)  that  the  extended  functions  satisfy 
(1.3),  (1.4)  with 


m 
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(9.33) 


p(-°°)  -  p(“),  c (-“>)  -  c(») 


The  corresponding  operator  in  X  will  be  denoted  by  A,  as  before. 
Property  (9.33)  implies  that 


(9.34) 


q+(y»A)  -  q(y,A),  q|(y»A)  *  q’(y,A) 


where  the  latter  are  defined  by  (9.20),  (9.24).  Moreover,  the  special 


solutions  $j(y,y,X)  of  §3  satisfy 


(9.35) 


<M-y.y.X)  *  4>„(y,y,A) 


<J>2(-y,y,X)  -  <|>3(y,y,A) 


and  it  is  easy  to  verify  that  the  eigenfunction  tk  (y,y)  of  A  is  even 

k  y 

(resp.,  odd)  when  k  is  odd  (reap.,  even).  It  follows  that  the  eigen¬ 


functions  of  A^  can  be  calculated  from  those  of  A^  by  the  rule 


(9.36) 


^(yty)  ■  /2  ^2k(y,y),  y  >  0,  k  =  1,2,- 


^(y.y)  *  &  i|>2k_l(y,y),  y  >  0,  k  -  1,2,< 


The  factor  Jl  is  to  renormalize  ^  from  R  to  R+. 

Concerning  the  generalized  eigenfunctions,  note  that  there  is  no 
ip0(y,y,X)  for  A  because  c(-<*>)  *  c(»)  and  (9.35)  implies  that 


(9.37) 


iM-y.y,A)  -  iMy.y.A) , 


It  follows  that  the  coefficients  R+,  T^  in  (1.33),  (1.34)  satisfy 


(9.38) 


R+(y,A)  -  R_(y,A),  T+(y,X)  -  T_(y,X) 
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and  Che  generalized  elgenfunc cions  of  A^  can  be  calculaced  from  chose  of 


by  Che  rule 


^°(y,y,A)  •  ’I'+Cy.u.A)  -  ^(-y.u.A) ,  y  >  o, 


(9.39) 


^(y.p.*)  -  ^(y.y.A)  +  ^(-y.u.A).  y  >  o. 

The  faccors  c^(y,X)  and  R^(y,A)  of  (9.23)  are  given  by 

(9.40)  c°(u,A)  -  c^U.A)  -  (p(»)/47rq(u,X))l/2  , 

(9.41)  R°  -  R+  -  T+ ,  R1  -  R+  +  T+. 

Theorem  9.1  can  be  deduced  from  Theorem  7.4  by  introducing  the 
operators 


(9.42) 


Jj  :  K(R^)  -  JC,  j  -  0,1. 


defined  by 

(9.43) 


J„  u(x,y) 


u(x,y),  (x,y)  €  R|, 

-u(x,-y),  (x,-y)  €  R^, 


and 


(9.44) 


Jj  u(x,y) 


u(x,y),  (x,y)  6  R|, 

u(x.-y).  (x,-y)  €  R+. 


J0  and  are  bounded  linear  operators  and  using  Che  fact  chat  A  has 
coefficients  satisfying  (9.32)  one  can  show  that  the  resolvents  of  A^  and 
A  are  related  by 


I 


(9.45) 


R(Aj,i;)  -  i  J*  R(A,£)  J y  j  -  0.1. 


From  chis  and  Stone's  theorem  (7.64)  It  follows  that 


(9.46) 


nj(y)  -  -f  J*  IKu)  Jjt  j  -  0,1. 


Theorem  9.1  follows  directly  from  these  relations  and  Theorem  7.4. 

Finite  Layers.  In  this  case,  with  a  suitable  choice  of 
coordinates  the  region  occupied  by  the  fluid  is  described  by  the  domain 


(9.47) 


R^  -  R3  n  { (x,y)  |  0  <  y  <  h} 


where  h  >  0.  The  case  of  a  fluid  layer  with  a  free  surface  at  y  **  0  and 
a  rigid  bottom  at  y  ■  h  will  be  discussed. 

The  acoustic  propagator  A  and  boundary  conditions  determine  a 
self ad joint  operator  Ah  in 


(9.48) 


^  =  L2(R^,c~2(y)p_1(y)dxdy), 


To  define  the  domain  of  A  let 


(9.49) 


Lj,0(R£)  -  Lj(R*)  O  (u  |  u(x,0+)  -  0  in  L2(R*)}. 


The  Dirichlet  condition  at  y  ■  0  will  be  enforced  by  requiring 

D(Ah)  C  L2 ’ 9 (R^) ■  The  Neumann  condition  at  y  *  h  will  be  interpreted  in 

the  generalized  sense  that 


(9.50)  |  }  (7,(p"’17u)v  +  p“l7u»7v}dxdy  ■  0  for  all  v  e  Lj’°(R^). 

iRh 


*  / 


Thus 
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(9.51)  D(Ah)  -  L‘*°(R^)  n  L*(A,R*)  n  (u  I  (9.50)  holds} 

and  A^u  ■  Au  for  all  u  €  D(A^).  As  before 

h  h* 

(9.52)  An  ■  An  >  0. 

The  corresponding  reduced  propagator  A^  in  JC(R^) 

-  L2(Rh,c-2(y)p-1 (y)dy)  (R^  -  {y  |  0  <  y  <  h})  is  defined  by 

(9.53)  D(aJ)  -  L‘0^)  n  {ip  |  (p_1Tp* ) '  6  «(R+),ip(0+)  «  0,(p~V)(h+)  -  0}, 
and  A^ip  -  A^p  for  all  ip  €  D(aJ)  .  One  has 

(9.54)  Aj  -  A**  >  c2p2 

as  before.  In  the  present  case  A^  is  a  regular  Sturm-Liouville  operator . 
Hence 

(9.55)  O(AJ)  -  o#(Aj) 

and  if  a£(p) ,  1  <  k  <  00 ,  denotes  the  eigenvalues  then  l£(p)  **•  ®  when 
k  -*■  00 .  Note  that  in  this  case  0^  *  R+  and  £2^  *  R2.  If  ipjj(y.P)  denotes 
the  corresponding  eigenfunctions  and  \p£(x,y,p)  «  (2tt)  !eip  X4»^(y.|p|) 

L 

then  the  spectral  family  (IIn(p)}  for  An  satisfies 

00  M 

(9.56)  (f,nh(p)f)  -  l  H(p  -  Xj(|p|))  |?J*(p)|2  dp 

k-1  'R2  k  k 

where 

(9.57)  fj(p)  -  L2(R2)-lim  fH  f  i|»J(*,y,p)  f (x,y)c~2 (y)p“l  (y)dxdy . 

M-~>  '  o  J  |  x|  <M 

These  results  can  be  proved  by  the  method  developed  in  §§2-7. 


tucmam  p*oi  autauwot  nu&s> 


§10.  Concluding  Remarks. 

The  spectral  analysis  of  the  acoustic  propagator  A  has  been 
developed  to  provide  a  foundation  for  the  study  of  acoustic  wave  propa¬ 
gation  In  stratified  fluids.  The  result  developed  in  this  report  can  be 
used  to  analyze  transient  sound  fields  in  stratified  fluids  following 
the  method  developed  for  the  Pekerls  model  in  [18].  It  can  also  be 
used  to  establish  a  limiting  absorption  theorem  and  corresponding  theory 
of  steady-state  sound  fields  in  stratified  fluids.  A  third  application 
is  to  the  analysis  of  the  scattering  and  trapping  of  acoustic  waves  by 
the  acoustic  ducts  produced  by  minima  of  c(y).  These  applications  will 
be  presented  in  separate  reports. 

A  number  of  extensions  of  the  theory  are  possible.  Thus,  other 
classes  of  density  and  sound  speed  profiles  could  be  studied.  Examples 
include  cases  where  p(y)  and/or  c(y)  tend  to  zero  or  infinity  at  y  *  ±°° 
or  at  finite  boundary  points  or  interior  points.  The  Weyl-Kodaira  theory 
is  applicable  to  all  such  operators  that  lead  to  selfadjoint  realizations 
of  A^.  Some  of  these  cases  will  be  of  interest  for  applications. 

Another  extension  of  the  theory  that  holds  great  interest  for 
applications  is  to  the  analysis  of  the  scattering  of  sound  by  obstacles 
immersed  in  stratified  fluids.  Mathematically,  this  problem  leads  to 
the  spectral  analysis  of  the  acoustic  propagators  for  stratified  fluids 
in  domains  exterior  to  bounded  sets.  The  analysis  presented  above  is  a 
necessary  preparation  for  such  a  study. 
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Appendix.  The  Weyl-Kodaira  Theory. 

The  general  Sturm-Liouville  operator  may  be  written 

(A.l)  L  <f>( y)  -  w"1  (y){-(p_1  (y)<J>*) '  +  q(y)<t>L 

The  basic  spectral  theory  of  such  operators  was  established  by  H.  Weyl 
[14]  and  K.  Kodalra  [10].  The  purpose  of  this  Appendix  Is  to  present  a 
version  of  the  Weyl-Kodaira  theory  that  Is  applicable  to  the  operator  A^ 
of  this  report. 

It  is  true  that  expositions  of  the  Weyl-Kodaira  theory  are 
available  in  [2,  3,  11]  and  a  number  of  other  textbooks  and  monographs. 
However,  in  these  and  most  of  the  book  and  periodical  literature, 
hypotheses  are  made  concerning  the  form  of  the  operator,  or  the  continuity 
or  differentiability  of  the  coefficients,  that  limit  the  applicability 
of  the  theory.  Thus  most  authors  assume  that  w(y)  =  1  and  many  take 
p(y)  i  1  as  well.  Moreover,  it  is  usual  to  assume  that  the  coefficients 
are  smooth  functions  or  at  least  continuous.  It  is  known  that  if  the 
coefficients  are  sufficiently  regular  then  L  can  be  reduced  to  the 
Schrodinger  form  L  <p  *  -<J>"  +  q(y)$  by  changes  of  the  independent  and 
dependent  variables  [11,  p.  2].  However,  this  technique  is  not 
applicable  to  operators  with  singular  coefficients.  Here  a  version  of 
the  Weyl-Kodaira  theory  is  presented  that  is  applicable  to  operators 
(A.l)  with  locally  integrable  coefficients.  The  concepts  needed  for 
this  extension  of  the  theory  are  available  in  the  classic  book  of 
Coddington  and  Levinson  [2]. 

The  operator  (A.l)  will  be  studied  on  an  arbitrary  interval 
I  ■  (y  |  — 09  <  a  <  y  <  b  <  +<*>}.  The  coefficients  will  be  assumed  to  have 
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the  properties 


(A.  2) 


p(y),  q(y),  w(y)  are  defined  and  real  valued  for 
almost  all  y  6  I, 


(A. 3) 


p(y)  >  0  and  w(y)  >  0  for  almost  all  y  €  I, 


(A. 4) 


p(y).  q(y),  w(y)  are  In  (I), 


where  C(I)  -  { f (y)  |  f  €  Lx (K)  for  every  compact  K  c  I}.  it  is 
natural  to  study  L  in  the  Hilbert  space  3f(I,w)  with  scalar  product 


(A. 5) 


(u,v) 


|  u(y)  v(y)  w 
>  I 


(y)dy. 


In  the  general  theory  of  singular  Sturm-Liouville  operators  two 
operators  in  3C(I,w)  are  associated  with  l.  The  first  is  the  maximal 
operator  Lj  defined  by 

[  D(Lj )  «  3C(I,w)  n  AC (I)  n  {u  I  p^u'  e  AC(I),lu  e  Jf(I,w)}, 


(A. 6) 


L,,u  *  lu  for  all  u  €  D(Lj). 


The  second  is  the  minimal  operator  L0  defined  by 


D(L0)  ■  D(Lj)  n  (u  |  (LjU.v)  -  (u.LjV)  for  all  v  €  D(Lj)}, 


(A. 7) 


L0u  *  lu  for  all  u  6  D(L„). 


It  can  be  shown  that  L.  is  densely  defined  and  closed  and  satisfies 


h  C  L;  -  L, . 


(A.  8) 
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It  follows  that  every  selfadjoint  realization  of  L  in  3f(I,w)  must  satisfy 
(A. 9)  L,CLCL,. 

•ff 

If  LQ  *  »  Lj  then  L  is  said  to  be  essentially  selfadjoint.  The 

classification  of  the  selfadjoint  realizations  of  L  by  means  of  boundary 
conditions  at  a  and  b  will  not  be  reviewed  here.  For  essentially 
selfadjoint  operators  no  boundary  conditions  are  needed  (Weyl's  limit 
point  case) .  The  operator  A^  of  §1  is  essentially  selfadjoint  since  its 
maximal  operator  is  selfadjoint  (cf.  (1.18),  (1.19)). 

The  Weyl-Kodaira  theory  provides  spectral  representations  of  the 
selfadjoint  realizations  of  singular  Sturm-Liouville  operators.  Each 
representation  is  derived  from  a  basis  of  solutions  of  L  Tp  “  Xip  and  a 
corresponding  2x2  positive  matrix  measure  m(X)  *  (m^k(X))  [3,  p. 

1337f f ] .  The  representation  spaces  are  the  Lebesgue  spaces  L2(A,m) 
associated  with  m,  with  norm  defined  by 


(A. 10) 


i  l  F  AX)  F.  (X)  m  (dX). 
'A  j,k-l  J  K  JK 


The  following  version  of  the  Weyl-Kodaira  theory  is  adapted  from 
[3,  pp.  1351-6]. 

Theorem  (Weyl-Kodaira) .  Let  L  be  a  selfadjoint  realization  of  L 
in  JC(I,w)  with  spectral  family  (II^(X)}.  Let  A  ■  (Xj  ,X2)  C  R  and  let 
(y,X)  (j  ■  1,2)  be  a  pair  of  functions  with  the  properties 

(A. 11)  *j(y,X)  e  C(I  x  A),  j  -  1,2, 

(A. 12)  The  pair  t^(y,X)  (j  *  1,2)  is  a  solution  basis  for 

L  ip  m  X\p  on  I  for  each  X  €  A. 
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Then  there  exists  a  unique  2  *  2  positive  matrix  measure  m  =  (m.^)  on  A 
with  the  following  properties. 


(A. 13) 


For  all  f  €  30(I,w)  there  exists  the  limit 


f(A)  ■  (f  l  (X)  ,f2(X))  -  L,  (A,m)-lim  f(y)  (ip  (y,A),i^2(y,A))w(y)dy. 

a^-*-a,b'-*-b  'a' 


(A. 14) 


(A. 15) 


(A. 16) 


The  mapping  U  :  JC(I,w)  L2  (A,m)  defined  by  Uf  ■  f  is  a 
partial  isometry  with  initial  set  n^(A)  JC(I,w)  and  final 
set  L2 (A,m) . 

The  inverse  isomorphism  of  L2(A,m)  onto  (A)  3f(I,w) 
is  given  by 


rti2  2 

(U*F)(y)  -  3f(I,w)-lim  £  <My, A)F.  (A)m..  (dA) 

Ml^A1,p2^A2  Jyl  j , k=l  J  J 


For  all  Borel  functions  ’l'(A)  on  R  with  supp  V  C  A, 
one  has 


U  D(Y(L))  -  L.  (A,m)  n  {f  |  'f(A)f  (A)  €  L2(A,m)},  and 


(U  *P(L)f)(A)  -  Y(A)f(A). 


Discussion  of  the  Proof.  The  theorem  is  proved  in  [3]  under  the 

hypotheses  w(y)  =  1,  p(y),  q(y)  S  C°°(I)  and  p(y)  >  0.  To  prove  it  under 

hypotheses  (A. 2),  (A. 3),  (A. 4)  one  may  first  prove  it  for  the  special 

(k-i )  k 

case  of  the  basis  <J>j(y»A)  that  satisfies  <J>j  7(c,A)  ■  6j  where  a  <  c  <  b. 

The  functions  <^(y,A)  are  entire  functions  of  A  and  the  theorem  can  be 
proved  by  the  classical  limit-point,  limit-circle  method  of  Weyl  as 
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presented  in  [2].  The  general  case  can  then  be  obtained  by  a  change  of 
basis  from  to  In  fact,  this  was  the  procedure  used  by 

Kodaira  in  his  original  paper  [10].  The  first  uniqueness  results  for  m 
are  due  to  E.  A.  Coddington  and  V.  A.  Mardenko  (see  [3]).  The  uniqueness 
proof  given  in  [3]  can  be  extended  to  the  case  treated  here. 

As  emphasized  by  Dunford  and  Schwartz,  the  utility  of  the  Weyl- 


Kodaira  theorem  is  due  to  the  possibility  of  using  different  bases 


for  different  portions  of  the  spectrum  of  L.  When  a  basis  has  been 


chosen  one  need  only  calculate  the  measure  m.  A  general  procedure  for 
doing  this,  due  to  E.  C.  Titchmarsh  [3,  p.  1364]  is  known  for  cases  in 
which  the  i/^(y,A)  have  analytic  continuations  to  a  neighborhood  of  A  in 
the  complex  plane.  However,  such  continuations  are  not  always 
available.  A  procedure  that  is  applicable  when  the  ^(y,A)  have  a 
one-sided  continuation  into  the  complex  plane  is  illustrated  in  §5  above. 
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